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ABSTRACT 


AN  INVESTIGATION  OF  THE  STABILITY  OF  A  SYSTEM  OF 
TWO  SHIPS  EMPLOYING  AUTOMATIC  CONTROL  WHILE  ON 
PARALLEL  COURSES  IN  CLOSE  PROXIMITY 

Charles  Natale  Calvano 

Submitted  to  the  Department  of  Naval  Architecture  and  Marine 

Engineering  on  21  May  1970  in  partial  fulfillment  of  the 

requirements  for  the  degrees  of  Master  of  Science  in  Ocean 
Engineering  and  Naval  Engineer. 

An  introduction  outlining  the  development  of  ship 
operations  in  close  proximity  for  underway  replenishment  is 
followed  by  a  survey  of  the  literature  dealing  with  the 
measurement  of  interaction  effects  between  ships  under  such 
circumstances.   A  brief  outline  of  the  theoretical  treatment 
of  ship  motions  in  open  water  is  given  and  followed  by  the 
modifications  to  the  theory  required  for  treatment  of  the  two 
ship  system.   The  applicable  equations  are  written  for  calm 
water  and  suitably  linearized.   The  techniques  of  treating 
the  differential  equations  as  algebraic  equations,  and  of" 
solving  the  equations  by  determinant  evaluation  are  described, 

Extension  of  available  interaction  and  open  water  data 
for  two  different  ships  to  a  full  data  set  for  a  single, 
modified  ship  is  carried  out,  and  the  numerical  values  of 
resulting  hydrodynamic  derivatives  are  given.   The  numerical 
solution  of  the  equations  and  the  extraction  of  stability 
roots  are  then  described  and  results  for  91  typical  trial 
data  sets  are  listed.   A  discussion  of  the  results  obtained 
is  followed  by  an  evoluation  of  the  correctness  of  the 
solution  method  employed. 

The  computer  programs  used  throughout  are  described  in 
detail  in  four  appendices. 

Thesis  Supervisor:   Martin  A.  Abkowitz 
Title:   Professor  of  Naval  Architecture 
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Chapter  I 
INTRODUCTION 

I.   Backgrou nd-Re plenishment  at  Se a 

In  their  efforts  to  achieve  a  degree  of  independence  from 
land  bases,  the  navies  of  the  world  have,  for  decades,  prac- 
ticed various  methods  of  underway  refueling  and  replenishment 
of  their  ships .   As  techniques  were  improved  and  experience 
was  g^lined,  these  operations  progressed  from  the  status  of 
infrequent,  makeshift  affairs,  to  that  of  routine  evolutions. 
These  replenishment  methods  have  ranged  from  actually  lashing 
together  the  supply  ship  and  receiving  ship  while  passing 
coal  and  provisions  from  one  to  the  other,  (one  assumes  this 
was  attempted  in  only  the  very  calmest  of  seas,  or  in  shel- 
tered roadsteads),  to  the  transfer  of  cargo  via  "vertical 
replenishment"  (VERTREP)  helicopters.   Other  methods  employed 
over  the  years  have  included  the  astern  refueling  method  and 
the  alongside  method  for  transferring  fuel  and  cargo. 

In  the  former,  the  receiving  ship  recovers  a  fuel  hose 
trailed  in  the  water  by  the  tanker  and  takes  on  fuel  while 
steaming  in  the  tanker's  wake.   In  the  alongside  method,  the 
receiving  ship  overtakes  the  supply  ship  on  a  parallel  course, 
laterally  separated  by  from  4  0  to  .18  0  feet.   The  transfer  is 
effected  while  the  ships  remain  abreast  at  this  distance. 

The  alongside  method  is  the  standard  of  the  U.S.  Navy, 
(augmented  by  VERTREP  for  cargo  transfer,  though  not,  of 
course,  for  refueling),  and  such  replenishment  operations 


have  become  almost  a  way  of  life  for  U.S.  combatant  ships. 
Unless,  and  until,  all  U.S.  combatant  ships  become  nuclear 
fueled,  this  method  will  probably  remain  in  wide  use. 

Of  course  there  are  other  situations  than  replenishment 
when  two  ships  might  wish  to  operate  in  close  alongside  prox- 
imity while  underway.   To  generalize  the  treatment  herein, 
the  forces  and  moments  arising  from  the  tension  in  the 
replenishment  rig  lines  is  not  included.   Thus  the  situation 
occurring  whenever  two  ships  must  steam  alongside  each  other 
at  close  quarters  is  that  which  is  treated  in  this  thesis. 

II .   Some  Problem s 

Although  alongside  replenishment  has  become  commonplace, 
this  type  of  close-in  steaming  has  not  ceased  to  be  a  dangerous 
undertaking.   There  are,  of  course,  dangers  arising  from  the 
lack  of  freedom  to  maneuver  to  avoid  other  shipping,  or,  in 
time  of  war,  the  actions  of  an  enemy.   But  the  ever-present 
source  of  danger  to  two  ships  steaming  close  together  on 
parallel  courses  is  this  very  proximity  of  the  ships,  causing 
each  to  constitute  a  collision  hazard  for  the  other. 

While  alongside,  the  ships  must  make  almost  constant 
rudder  and  engine  speed  adjustments  to  maintain  position. 
Because  the  ships  are  close,  the  presence  of  each  affects  the 
forces  and  moments  on  the  other  through  the  influence  of  each 
on  the  pressure  field  and  flow  patterns  around  the  other. 
(See  Figure  1-1) .   Thus,  neither  ship  will  respond  as  it  does 
in  open  water  where  these  influences  are  not  present.   And 
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Pressure  Field  about  a  Ship  in  Open  Water 
(From  Reference  A-l) 
FIGURE  1-1 

as  the  ships  alter  their  relative  positions,  the  direction  and 
magnitude  of  these  influences  are  continuously  altered. 

The  repeated  corrections  required  over  a  period  of  up  to 
several  hours,  introduce  hundreds  of  opportunities  for  human 
error  or  momentary  lack  of  attention  to  give  rise  to  damaging 
consequences  in  the  form  of  a  collision.   Any  method  or  device 
which  could  serve  to  reduce  the  probability  of  such  errors 
occurring  would  serve  to  lessen  a  primary  hazard  for  two  ships 
operating  in  this  way. 


11 


'. 


III.   Automatic  Control 

The  introduction  of  automatic  rudder  angle  and  engine 
speed  control  to  the  alongside  operation  is  an  attractive 
prospect.   Such  a  control  system  would  employ  sensors  to 
measure  longitudinal  and  lateral  separation  of  reference 
points  on  the  two  ships,  (and  their  rates  of  change),  and 
compare  these  values  to  desired  equilibrium  values.   It  would 
then  generate  rudder  and  engine  orders  to  reduce  the  deviations 
to  acceptable  magnitudes .   Such  a  system  could  be  a  boon  to 
conning  officers  and  helmsmen  by  relieving  them  of  much  of  the 
tedium  of  maintaining  position  alongside. 

A  theoretical  investigation  of  such  a  system  is  the 
subject  of  this  thesis.   The  problem  will  be  approached  from 
a  mathematical  standpoint  as  a  two-body  system  of  ship  motions. 
The  hardware  required  to  make  it  operable  will  not  be  examined 
or  described.   It  will  be  assumed  that  a  sensor,  (doppler  sonar, 
for  instance) ,  is  available  to  continually  measure  lateral  and 
longitudinal  ship  separations  and  their  rates  of  change. 

Once  the  deviation  from  desired  ship  separation  and  heading 
is  known,  an  automatic  control  system  can  order  rudder  and 
engine  speed  changes  to  reduce  the  deviations.   The  various 
constants  associated  with  a  theoretical  control  system  meeting 
this  description  will  be  a  part  of  the  mathematical  model  of 
the  two  ship  system.   The  stability  of  the  system  will  be 
investigated,  to  determine  if  it  appears  feasible  to  employ 
such  an  automatic  control  system  to  keep  the  ships  in  position 
alongside.   (By  stability  is  meant  the  tendency  of  the  two 
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ship  system  to  return  to  an  equilibrium  condition  after 
receiving  a  small  disturbance  from  that  equilibrium  condition. 
If  the  system  "returns,  or  tends  to  return,  to  the  equilibrium 
condition  when  the  disturbance  is  removed,  then  it  is  stable". 
[A-3]*) 

IV.   Historical  Aspects---The  Literature 

A.   Taylor. 

As  far  as  the  author  has  been  able  to  determine,  the  first 
systematic  attempt  to  measure  interaction  phenomena  for  two 
ships  on  parallel  courses  in  open  water  was  made  by  David  W. 
Taylor,  the  ubiquitous  Naval  Constructor.   He  conducted  a 
series  of  tests  utilizing  four  models,  each  of  3,000  pounds 
displacement,  at  the  Model  Basin  in  Washington,  in  1909.  [A--2] 

Taylor's  tests  were  conducted  by  towing  the  models  in 
pairs,  either  abreast  of  one  another,  or  with  one  model  a 
fixed  distance  ahead  of  the  other.   The  models  were  not  self- 
propelled,  and  did  not  have  propellers  affixed.   They  were 
towed  under  towing  carriages,  with  their  bows  and  sterns 
prevented  from  translating  laterally,  but  free  to  move  vertic- 
ally, changing  draft  and  trim  while  being  towed.   Water  depths 
were  "many  times  the  draught  of  the  models".   Attraction  and 
repulsion  forces  were  measured  at.  the  attachment  points  near 
the  bow  and  stern  of  each  model.   Table  1-1  lists  the  charac- 
teristics of  the  models  employed  and  Figure  1-2  is  a  typical 
plot  of  the  results  Taylor  obtained. 


*Letter-number  citations  in  brackets  refer  to  the  references 
listed  at  the  end  of  the  body  of  this  thesis. 
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Model 

Beam  (B) 

Draft  (H) 

i1. 

Prismatic 

Midships 

Block 

# 

ft. 

ft. 

H 

Coeff  C 

Coeff  C 
m 

Coeff  C. 
b 

834 

3.692 

1.263 

.560 

.560 

.900 

.504 

838 

3.503 

1.198 

2.92 

.560 

1.00 

.560 

858 

3.586 

.957 

3.75 

.74  0 

.926 

.685 

866 

2.778 

1.235 

2.25 

.740 

.926 

.685 

TABLE  1-1 
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Figure  1-2 
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In  conducting  his  tests,  Taylor  found  that  he  h 
difficulty  towing  the  models  on  a  straight  course.   In  view  of 
the  fact  that  towing  carriages  were  used,  this  is  hfird  to 
understand,  but  was  explained  no  further.   As  a  result  of  this 
difficulty,  he  felt  that  the  forces  measured  could  not  be 
regarded  as  highly  accurate,  but  felt  that  they  showed  the 
tendencies  and  general  nature  of  the  interaction,  or,  to  use 
his  term,  "suction" ,  phenomena. 

The  sequence  of  events,  when  vessel  B  overtakes  vessel  A 
on  a  parallel  course,  seemed  to  be  about  as  follows: 

1.  As  B's  bow  overlaps  A's  stern,  a  small  repulsion 
force  is  felt  at  both  B's  bow  and  stern. 

2.  As  B  continues  to  overtake  A,  B's  bow  is  attracted, 
while  its  stern  is  repelled,  causing  B  to  tend  to 
turn  into  A. 

3.  As  B  draws  abreast  of  A,  both  B's  bow  and  stern 
are  drawn  in  toward  A. 

4.  As  B  pulls  up  farther,  the  force  on  its  bow  switches 
to  repulsion,  while  its  stern  is  still  attracted 

to  A. 

5.  As  B  draws  farther  forward,  to  where  its  stern  is 
approximately  abreast  of  A's  midlength,  both  B's 
bow  and  stern  are  repelled. 

Taylor  noted  further,  that  the  occurrence  of  these  events 
was  somewhat  more  consistent  with  the  "finer"  models,  while 
the  results  with  the  "full"  models  tended  to  be  e~   tic.   The 
effect  of  the  "suction"  forces  on  the  models'  resistance  vras 
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not  measured,  but  Taylor  did  note  that  the  side  forces  "varied 
with  speed  as  the  resistance  of  the  model"  ,"  that  is  the  inter- 
action forces  increased  as  the  square  of  the  velocity. 

It  is  evident  from  the  discussions  accompanying  Taylor's 
paper  that  there  was  disagreement  within  the  Naval  Architect 
community  concerning  the  significance  and,  in  fact,  even  the 
existence,  of  interaction  forces  for  full  scale  ships. 

B.  Gibson  and  Thompson. 

In  1913,  Gibson  and  Thompson  conducted  a  series  of  full 
scale  tests  with  small  steam  launches  in  the  Thames  River 
near  London,  to  verify  the  existence,  and  measure  the  effect 
of  interaction  phenomena.   [B-2]   They  did  note  effects 
similar  in  type  to  those  described  by  Taylor,  but  due  to 
instrumentation  difficulties,  their  report  contains  no 
quantitative  information  of  a  useful  nature  in  considering 
interaction  effects  for  ships. 

C .  Robb . 


A  disastrous  collision  between  the  liner  QUEEN  MARY  and 
the  destroyer  H.M.S.  CURACOA  in  19  4  2  led  to  a  series  of 
experiments  reported  by  A.M.  Robb  in  1949.   [A-8]   In  these 
tests,  the  models  used  were  of  these  two  ships  and  made  to  a 
1:56  scale.   This  resulted  in  the  QUEEN  MARY  model  being 
"rather  more  than  18  feet"  in  length,  and  that  of  the  CURACOA 
"rather  more  than  8  feet".   In  these  experiments,  both  models 
were  self --propelled,  with  first  one,  then  the  other,  being 
constrained  to  run  on  a  straight  course  while  the  other  passed 
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on  a  nearly  parallel  course.   Robb  again  confirmed  the  exist- 
ence, of  significant  interaction  effects  drawing  the  smaller 
ship  into  the  larger,  in  some  circumstances/  and  his  article 
is  accompanied  by  several  sketches  showing  the  trajectories  of 
the  ships  for  various  speeds  and  lateral  separations.   The 
results  presented  are  largely  qualitative  in  nature,  and  do 
not  disagree  with  those  obtained  by  Taylor  40  years  earlier. 
But  again,  in  the  discussions  of  Robb ' s  paper,  the  fact 
that  the  very  existence  of  interaction  effects  was  doubted 
in  some  quarters  was  brought  out.   [B-5]   Professor  Prohaska, 
in  his  discussion,  stated,  "The  question  of  interaction  between 
passing  vessels  has  always  given  rise  to  interesting  exchanges 

of  opinion a  great  many  experts,  amongst  whom  are  numerous 

seafaring  people  with  great  experience,  deny  the  existence  of 
interaction,  or  at  least  think  it  without  importance  in 
unrestricted  water."   (Prohaska,  however,  did  not  share  this 
view. ) 

D.   Restricted  Water  Studies. 

Over  the  years,  there  have  been  numerous  other  experiments 
and  theoretical  developments  concerning  the  interaction  between 
a  ship  and  a  stationary,  constant-geometry  object,  such  as  a 
canal  wall  or  bank.   References  A~6,  A--10,  B--6,  B-10,  B-ll, 
B-12,  B-13  and  B-18  all  deal  with  this  phenomenon.   A  ship, 
being  a  large  object,  can  appear,  in  its  effects,  as  a  wall 
would  appear  to  another  ship  passing  near.   But  of  course, 
the  geometry  of  the  ship-object  is  different  than  when  the  object 
is  a  wall.   And,  also,  neither  ship  is  immovable,  as  is  the  wall. 
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Most  of  the  investigations  concerning  interaction  in 
canals  or  channels  include  shallow  water  conditions.   The 
shallow  water  tends  to  emphasize  the  interaction  effects  caused 
by  the  wall.   This  investigation,-  though,  is  concerned  with  the 
interaction  between  ships  in  deep  water. 

E.   Newton. 


The  single  study  of  interaction  effects  which  proved  to 
be  of  greatest  use  to  the  author  in  attempting  to  obtain 
quantitative  information  concerning  the  interaction  forces  and 
moments  was  that  listed  as  reference  A--1,  by  R.N.  Newton  and 
reported  at  the  First  Symposium  on  Ship  Maneuverability 
conducted  at  the  David  Taylor  Model  Basin  in  19  60. 

Newton  reported  on  experiments  conducted  at  the  Admiralty 
Experiment  Works  (AEW)  in  England,  between  1945  and  1948  and 
which  were  inspired  by  the  growth  of  alongside  replenishment 
experienced  during  the  war.   In  his  paper,  Newton  acknowledges 
the  pioneering  work  of  Taylor  in  this  field  and  laments  the 
near-total  lack  of  work  in  this  area  in  the  years  between 
Taylor's  tests  and  the  AEW  efforts,  although  he  noted  the 
wealth  of  information  available  to  him  concerning  interaction 
effects  in  restricted  channels  and  shallow  water. 

Newton's  experiments  were  with  models  of  a  battleship, 
H.M.S.  KING  GEORGE  V  and  a  replenishment  supply  ship,  R.F.A. 
OLNA.   (See  Table  1-2  for  ship  characteristics.)   The  tendencies 
exhibited  by  the  models  were  confirmed  by  full  scale  trials  at 
sea. 
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K.G.V 

OLNA 

Length  on  W.L. ,    ft. 

740 

.  547 

Beam,  ft. 

103 

7  0 

Draft,  ft. 

29.3 

30 

Displacement,  tons 

39,890 

23,570 

Block  Coefficient 
Corresponding"  water  depth 

.611 

.714 

75 

fathoms 

TABLE  1-2 

Tests  were  conducted  with  both  towed  and  self-propelled 
versions  of  the  models.   The  models  were  allowed  complete 
vertical  freedom,  but  constrained  from  lateral  motion.   Forces 
and  moments  were  measured  at  longitudinal  separations  of  the 
midships  sections  of  the  ships  from  600  ft.  aft  of  the  abreast 
position  to  600  ft.   forward  of  this  position,  at  two  lateral 
separations;  50  and  100  ft.  side  to  side  distance.   With  the 
ships  in  the  abreast  position,  the  side  to  side  distance  was 
varied  from  25  to  14  0  ft.   (All  distances  are  the  corresponding 
full  scale  values.) 

Figures  1-3  and  1-4  show  the  longitudinal  variation  of 
the  forces  and  moments  experienced  by  OLNA  at  the  two  lateral 
separations  and  Figure  1-5  shows  the  lateral  variation  with 
the  ships  abreeist.   These  data  formed,  the  basis  for  quantitative 
interaction  forces  and  moments  used  in  this  investigation,  as 
OLNA  is  not  radically  different  in  hull  characteristics  from 
the  MARINER  class  vessel  used  for  this  study.   In  the  figures, 
OLNA  is  the  ship  being  overtaken .   The  agreement  with  the 
results  obtained  by  Taylor  is  apparent. 
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F. Strom-Tejsen  and  Chislett. 

In  making  a  theoretical  investigation  such  as  this  one, 
it  is  necessary  to  have  realistic  values  for  the  hydrodynamic 
derivatives  employed  in  the  mathematical  model.   (The  use  of 
these  derivatives,  and  the  mathematical  model  are  discussed  in 
Chapter  II,  THEORY.)   For  two  ships  steaming  in  proximity, 
values  of  both  the  open  water  and  interaction  derivatives  are 
needed,  as  forces  and  moments  v/hich  are  always  present  for  a 
ship  in  open  water  will  be  present  along  with  those  caused  by 
the  interaction  phenomena. 

Thus,  in  this  investigation  it  was  necessary  to  find 
sources  for  typical  values  of  both  kinds  of  hydrodynamic 
derivatives.   As  mentioned  earlier,  Newton's  work  with  KING 
GEORGE  V  and  OLNA  was  the  only  available  source  of  quantitative 
interaction  information.   It  was,  therefore,  necessary  to  find 
values  of  the  open  water  derivatives  for  a  ship  similar  to 
either  KING  GEORGE  V  or  OLNA.   Strom-Tejsen  and  Chislett  in 
reference  A~5  determined  most  of  the  open  water  hydrodynamic 
derivatives  of  interest  for  a  MARINER  class  merchant  ship. 
(See  Table  1-3  for  MARINER  characteristics.)   In  reference  A--4, 
Strom-Tejsen  used  these  derivatives  to  theoretically  predict 
the  open  water  maneuvering  behavior  of  a  MARINER,  and 
subsequent  sea  trials  proved  these  values  to  be  quite  accurate. 
Fortunately,  OLNA  and  the   MARINER  class  hulls  are  sufficiently 
similar  to  expect  the  hydrodynamic  derivatives  of  one  to 
reasonably  resemble  those  of  the  other.   Therefore,  in  the 
development  of  the  results  of  this  thesis,  the  open  water 
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derivatives  of  the  MARINER  were  used,  and,  to  account  for 
interaction  effects,  derivatives  derived  from  the  interaction 
force  and  moment  data  of  the  OLNA  experiment  were  added  to 
these  to  obtain  the  entire  list  of  needed  derivatives.   Thus, 
the  ship  type  investigated  is  a  somewhat  modified  MARINER 
class  ship. 


MARINER  Characteristics 


Length,  ft. 
Beam,  ft. 
Draft,  ft. 
Displacement,  tons 
Block  Coefficient  C 


b 


Prismatic  Coeff.  C 


p 


Midships  Section  Coeff.  C 


m 


528.5 
76.0 
29.7  5 
20,900 

.6125 
.6246 
.9807 


TABLE  1-3 
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Chapter  II 
THEORY 


1^ B  ackground 

The  derivation  of  a  mathematical  model  representing  the 
steering  and  maneuvering  of  a  surface  ship  in  open  water  is 
given  by  Abkowitz  in  reference  A--3.   In  its  general  form,  the 
model  accounts  for  non-- linear,  as  well  as  linear,  effects, 
and  in  his  computerized  prediction  of  open  water,  surface  ship 
maneuvers,  Strom-Tejsen  uses  the  Abkowitz  model  with  up  to 
third  order  terms.  [A- 4] 

This  mathematical  model,  in  linearized  form,  serves  as 
the  basis  for  the  analysis  made  in  this  thesis,  with  appro- 
priate extension  of  the  model  to  include  the  interaction 
effects  which  are  not  present  in  open  water,  but  which  are 
important  for  two  ships  on  parallel  courses  in  close  proximity, 
The  model  used  herein  also  is  modified  to  represent  the  two- 
ship  system,  rather  than  the  motions  of  a  single  ship,  and 
accounts  for  the  fact  that  both  ships  are  free  to  move  in  the 
horizontal  plane.   (For  a  more  complete  description  than  the 
outline  which  follows,  of  the  derivation  of  the  open  water 
mode],  reference  A--3  is  suggested). 

II.   Reference  Axes 

The  reference  axes  used  constitute  an  orthogonal,  right- 
handed  system  which  takes  advantage  of  the  existence  of  a 
vertical  plane  of  symmetry  through  the  center line  of  the  ship. 
The  selection  of  the  axis  system  illustrated  introduces 
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simplifications  in  the  equations  of  motion  by  virtue   of  the 
fact  that  the  plane  of  symmetry  of  each  ship  is  parallel  to 
the  principal  axes  of  inertia  of  the  ship. 


G       o 


■i-Q, 


-{>       +0 


■i-f/j 


STBD. 
> 


L/£v 


±3%, 


+  0 


*+fr 


(See  Table  I1--1  for  definitions  of  all  symbols  used 

in  this  chapter) 

Reference  Axis  System 
FIGURE  II -I 


III.   Equations  of  Motion  for  a  Ship  Moving  in  the  Horizontal 
Plane 

A  general  form  of  the  equations  of  motion  for  a  body 
v/hich  is  allowed  to  move  in  all  the  six  degrees  of  freedom 
is  given  by  Abkowitz. 

Normally,  when  dealing  with  steering  and  maneuvering  of 
surface  ships  in  open  water,  the  primary  motions  are  considered 
to  take  place  in  the  horizontal  plane,  and  vertical  motions  are 
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neglected.   As  this  is  also  the  case  for  the  two-body  system 
of  interest,  where  only  horizontal  motions  are  of  concern, 
the  equation  in  Z,  the  vertical  force  equation,  can  be 
dropped.   Further,  choice  of  the  axis  system  parallel  to  the 
plane  of  symmetry  of  the  ships  eliminates  the  cross  terms  of 
inertia,  and,  if  it  is  assumed  that  the  center  of  gravity  of 
each  ship  is  in  its  centerline  plane,  (as,  indeed,  can  be 
expected),  the  quantity  y   is  zero  for  both  ships.   The 
usual  assumption  that  rolling  and  pitching  motions,  if  not 
extreme,  have  little  influence  on  steering  and  maneuvering  is 
also  made,  allowing  the  equations  for  K  and  M  to  be  neglected, 
and  making  the  variables  p  and  q  zero.   (See  Table  II-l  for 
definitions) .   The  equations  of  motion  for  a  single  surface 
ship  moving  in  the  horizontal  plane  thus  become: 

2 

a)  X  =  m[u-rv-x  r  ] 

g 

b)  Y  =  m  [v-l-ru+x  r] 

g 

c)  N  =  I  f+mx  (v+ru)  (2-1) 

z     g 

(where  the  notation  v  =  3v/3t,  is  used) 

The  left-hand  sides  of  these  equations  represent  the  sum 
of  all  forces  in  the  x  and  y  directions  and  all  moments  about 
the  z  axis,  respectively.   The  right-hand  sides  are  the 
dynamic  response  terms  associated  with  the  ship  and  its  motions 

IV.   Taylor  Expansion  of  Forces  and  Moments 

The  forces  and  moments  on  the  left-hand  sides  of  the 
equations  of  motion  can  be  expressed  as  functions  of  the 
properties  of  the  body,  properties  of  the  fluid  and  properties 
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TABLE  II-l 


Definitions  of  Symbols 


x  ,v  ,z        Reference  coordinates  fixed  with  respect  to 
o  2  o      o 

earth . 

x,y,z  Reference  coordinates  fixed  in  ship. 

u,v,w  Velocities  in  direction  of  x,  y  and  z  axes, 

respectively . 

0/.O,\O  Roll,  pitch  and  yaw  angles  about  x,  y  and  z 

axes,  respectively. 

p,q,r  Angular  velocities  about  x,  y  and  z  axes, 

•  •  • 

respectively.       (p=0 ,    q~8  ,    r~\p)  . 

Rudder  deflection  in  radians  (port  rudder  is 
positive) . 

n  Propeller  speed  in  revolutions  per  second. 

X,Y,Z  Sums  of  forces  acting  on  ship  in  x,  y,  z 

directions,  respectively. 

K,M,N  Sums  of  moments  acting  on  ship  about  x,  y,  z 

axes,  respectively. 

a  Longitudinal  separation  of  ships  measured 

from  Ship  A  midships  section  to  Ship  B 
midships  section. 

Lateral  separation  of  ships  measured  from 
Ship  A  centerline  to  Ship  B  centerline. 

Y  Typical  of  notation  used  for  hydrodynamic 

R 

derivatives.   In  this  case:   9Y/3u  for  Ship  B, 

measured  when  all  variables  but  u  are  at 
equilibrium  values. 

At  ,  At„        Control  system  time  lags  associated  with  rudder 

and  propeller  response,  respectively.   (Subscripts 
A  or  B  indicate  applicable  ship,  as  At..  ). 

A 
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k.  ,  k. 

Control  system  constants,  i=l  to  6,  for  Ship  A 

XA   XB 

or  B.   Constants  associated  with  control  system 

response  to  deviations  from  equilibrium,  as 

6 A  =  klAA^JA  ' 

x  /  Y  f  z 

Distance  of  center  of  gravity  of  ship,  (G), 

g   g   g 

from  the  axis  origin,  (0),  in  x,  y,  z  directions 

respectively. 

mA'  mB 

2 

Mass  of  Ship  A  or  B  in  units  of  .  (lb  sec  /ft) 

<5 

The  operator  d/dt,    as  JB  u  =   3u/8t 
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of  the  motion.   When  considering  specific  hull  forms,  and 
using  generally  accepted  scaling  laws ,  the  forces  and  moments 
may  be  considered  functions  of  only  the  motion  and  orientation 
parameters  of  the  bodies.   (When  dealing  with  automatic 
control  of  rudder  angle  and  propeller  speed,  they  are  also 
functions  of  rudder  deflection,  5,  and  propeller  speed,  n.) 

The  reference  system  for  the  two-ship  system  at 
equilibrium  and  an  illustration  of  the  variables  involved 
are  shown  in  Figure  II~2. 

When  considering  motion  in  an  unrestricted  horizontal 

plane,  it  is  clear  that  no  forces  or  moments  arise  from  a 

change  in  ship  location  with  respect  to  a  fixed  geographical 

reference  point,  allowing  x   and  y   to  be  discarded  as 
1  J   o     Jo 

variables  of  interest.   Similarly,  the  equilibrium  yaw  angles, 

((<//, )Q  and  O/'b^o^'  or  dr^^'t  angles,  can  be  dropped  from  the 

list  of  independent  variables  because  forces  and  moments 

arising  from  a  drift  angle  are  accounted  for  by  the  terms  in 

u  and  v.   (If  there  were  no  drift  angle,  u  would  coincide  with 

velocity  made  good  in  the  desired  direction ,"  and  v  would  be 

zero.   Thus  the  forces  and  moments  that  (ijO   and  (ihn)       cause 

4A  o       YB   o 

to  be  introduced  are  taken  up  in  the  terms  in  u  and  v.   The 
inter-relation  of  drift  angle  and  the  independent  variables, 
u  and  v,  can  be  seen  from  tan (drift  angle)  =  v/u). 

It  is. also  possible  to  drop  a,  a,  3  and  £  as  independent 
variables  because  of  the  simple  trigonometric  relations 
between  them  and  v,  v,  u  and  u.   The  usual  assumption  that 
forces  and  moments  produced  by  6 ,  6  and  n  are  negligible,  in 
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terms  of  ship  control,  is  also  made.   (These  forces  and 
moments  are  not  necessarily  negligible  when  considering  forces 
on  the  rudder  stock,  propeller  blades,  etc.   But  they  do  have 
a  negligible  effect  on  ship  motion.   (See  reference  A-4, 
pg.  7).) 

As  Abkowitz  has  shown,  the  left-hand  sides  of  the 
equations  of  motion,  comprising  the  functions  describing  the 
forces  and  moments  acting  on  the  ship,  can  be  developed  into 
useful  form  for  analysis  by  use  of  the  Taylor  Series  expansion 
for  a  function  of  several  variables.   The  force  or  moment 


caused  by  a  change  in  any  variable  is  expressed  as  the  product 
of  the  derivative  of  the  force  with  respect  to  that  variable, 
(with  all  other  variables  at  equilibrium  values) ,  and  the 
change  in  the  variable.   In  the  linearized  form  of  the 
equations,  only  the  linear  terms  are  retained,  though  the 
expressions  can  be  written  to  any  desired  degree  of  accuracy, 
by  retaining  sufficient  terms  in  the  expansion.   (For  example, 
the  Y  force  caused  by  a  change  in  u  is  written  as  (3y/3u)x(Au) 
where  (3Y/3u)  is  the  value  of  (3Y/3u)  with  all  other  variables 
at  their  equilibrium  values.   The  derivatives,  such  as  (3Y/3u) 
are  the  hydrodynamic  derivatives  and  can  be  experimentally 
determined. ) 

Employing  standard  notation,  where  Y   represents  (3Y/3u) 
at  equilibrium,  the  Y  equations,  for  example,  for  the  two 
ships,  retaining  only  the  linear  terms  of  the  Taylor  expansion, 
become : 
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Ship   A: 

a)    YaflAa      +y6AAB+YUAAVYuAAClA+yvAAV 

••"+Y6„ASA+Yn  AnA  ==  mA [ VrAuA+xg  *A] 

A       A  ^A 


b)  Y   AcH-Y„  A3  +Y   Au  -l-Y.  Au_+Y   Av„+. 
aB     3B      UB   B   UB   B   VB   B 


+Y.    A<5    +Y^    Anp    =  m    [v   +r   u    +x      r    ]  (2-2) 

6  B       nR       B  BBBB      g       B 


The  left-hand  sides  of  the  X  and  N  equations  are  similarly 
expanded  for  both  ships. 

V.   Linearizing  the  Dynamic  Response  Terms 

Since  only  linear  terms  have  been  retained  in  the 
expanded  force  and  moment  expressions  on  the  left-hand  sides 
of  the  equations,  a  similar  linearization  of  the  dynamic 
response  terms  on  the  right-hand  sides  must  be  made. 
For  all  variables,  the  following  holds  true: 

variable  =  (variable)   +  Avariable  (2-3) 

where  the  subscript  "o"  refers  to  the  equilibrium  value. 

For  all  time  derivatives  of  the  variables,  (v,  ip ,  etc.), 
the  equilibrium  value  is  zero,  as  the  equilibrium  condition 
is  unchanging  in  time.   Thus  for  these: 

variable  =  (variable)   +  Avariable  =  Avariable       (2-4  ) 
(e.g.   u  =  (u)   -!-   Au  =  Au) 

Making  use  of  (2-3)  and  (2-4),  and  by  neglecting  second 
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V 


order  terms  consisting  of  products  of  small  quantities,  such 
as  (Ar)«.(Au)  ,  the  linearized  right-hand  sides  of  the  Y  equations 
become : 


Y  =  m[v  +  ur-l-xr] 


(2-5) 


Similar  developments  are  applicable  to  the  X  and  N 
equations  and  are  carried  out  in  detail  in  [A- 3] . 

VI.   Linearized  Equations  of  Motion 

Taking  into  account  all  developments  outlined  thus  far, 
the  linearized  equations  of  motion  become: 

Ship  A: 

a)  X   Aa+X0  A3+X   Au  +X.  u  -l-X   Av,+X.  v„ 
aA     3A     UA   A   UA  A   VA   A   VA  A 


+X-    iJu+Xv    $.+Xj5    A6    -l-X      AnA    =   m.u,. 
^A  ^A   A       6A  A      A  A   A 


b)  Y   Aa+Y„  A6+  Y   Au-hY.  u  +Y   Av„+Y.  v, 
aA     ^A      UA   A   UA  A   VA   A   VA  A 


+Y*  i>.+Y.,    iJL+Y..  A6+Y   An,  =  m,  [v^+U   ^  -hx   $  1 
^A  A   VA   6A   A   nA   A     A   A   °A  A   gA  A 


c)  N   Aa+N„  A8+N   Au  -l-N.  u„+N   Av  +N.  v„ 
aA     eA     uA   A   uA  A   vA   A   vA  A 


+N'  ik-i-N.,  i|L-i-Nr  A6,+N   An,  =  I   ffi  +m,  x   (v„+u   iL) 
V  A   V A   6A   A   nA   A     2A  A   A  ^A   A   °A  A 


d) 
e) 
f) 


Corresponding  Equations  for  Ship  B 


(2-6) 
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VII .   Automatic  Control 

The  automatic  control  system  used  to  maintain  each  ship 

in  position  alongside  is  assumed  capable  of  deflecting  the 

rudder  and  changing  propeller  speed  in  response  to  measurement 

of  the  ship's  heading  angle  error  from  equilibrium,  A  if; ,  and 
V 

the  separation  parameter  errors  from  equilibrium,  Aa  and  A3, 
as  well  as  the  rates  of  change  of  these  variables.   The  rudder 
control  is  sensitive  to  if> ,  i|> ,    $  and  3;  the  propeller  control 
to  a   and  a . 

a)  6  =   k  Ai|M-k J>+   k-Aa+k.g 

b)  n  =  kcAa+k,a  (2-7) 

(where  k.  are  the  control  system  constants) 
There  will  be  time  lags,  (At),  in  both  control  systems, 
representing  the  lapse  between  ordering  a  control  response 
and  attainment  of  the  response.   Thus  rudder  deflection  at 
time  t  is  proportional  to  \p    at  time  t,  ,  (where  t  =  t.+At). 
Therefore : 

6(t)  =  kip(t-At)+k  i];(t-At)+k  3(t-At)+k  .3  (t-At)       (2-8) 

In  the  linearized  equations  of  motion  ifi  and  3  are 
functions  of  t.   Thus  these  variables  in  (2-8)  must  be 
obtained  as  functions  of  t,  which  is  done  through  another 
Taylor  expansion,  carried  out  in  detail  in  [A~3] .   If  the 
time  lags  are  assumed  small,  only  linear  terms  need  be  kept, 
and  the  automatic  control  system  equations  for  6  and  n 
become : 


35 


a)    6    -    [k1+(k2-k1At1)J9    -k^t^  2]  Aij> 
+  [k3+(k4-k3At1)^  --k^t^O  2]A(3 


b)    n   -    [kc+(kr-kcAt_)c9   ~k^At„«02]Aa  (2-9) 

(Where  <^  is  the  operator  d/dt;    e.g.  JO    ip  =   82^/9t2) 

Here  At,  is  the  time  lag  associated  with  rudder  response 
and  At„  that  associated  with  propeller  response. 

These  expressions  for  6  and  n  are  substituted  in  the 
equations  of  motion. 

VIII.   Equations  Relating  a  and  3  to  u ,  v  and  ty 

Examination  of  Figure  II--2  shows  the  following 
relationships  to  hold: 

a)    a   =   uBcos^B-uAcosipAH-vAsim|;A-vBsini|)B 


b)    3    =   u^sin^-u,  sinib-.+v^cosib-v^  cosily  (2-10) 

B  yB      A  yA       B  rB       A  YA 


Since    ty      and      ip_    are    small,    cosip      is    1.0    and    cosip      ~    1.0 
A  B  A  B 

And    since   v,  <<u,    and   v<<u^,    and    sinil)      and    sinili      are    small, 
A         A,  B         B  yA  yB 

the   products   v    sinif;      and   v    sinijj      are    second   order.      Thus 


(2-10) ,  a) ,  becomes: 

a  =  uR-uA  =  ((^B)0+AuB)"((UA)0+AUA) 


where  (UB)Q  =  u0  =  ^ua^0'  and 


a  =  AuB-AuA 
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LOa^    jc  w 


Also,    since    \\>      and    i>      are    small,    siivp      z    ty      and 


sinij^    ~i/-'    .       Thus,     (2-10),    b)  ,    becomes 


o  s  VbWY"» 


=  ♦bNV'*a,^io'iVVVVaV1"s 

2nd    order 


Therefore,  the  equations,  (2-10),  become 


a)    a    ~Au   -Au 


b)  3  VVWaYAva 


(2-11) 


^ 


IX.   The  Complete  Set  of  Linear  Equations  of  Motion  f or  the 

Two  Ship  System 

Incorporating  all  the  foregoing  discussion,  the  linear 
equations  of  motion  for  the  tv;o-ship  system  emerge: 

Ship  A:  (2-12) 

a)    X      Acx+X.    AB+X      Au_+X.    u_+X      Av^+X.    v+X*    iL 
aA  3A  uA      A      uA  A      vA     A      vA  A      ^A 


+Xv   $A+x6    [kj    +(k2   -k1  Atx    )jQ     -k2   Atx  cQ2]Aij 


A      A    A     A    A    A 


A    A 


A 


+X   [k  +(kA    -k   At   )«0  -k   At  £  2]A3 
°A    A    *A    A    A       A   XA 


+X   [k_  +(k.  -k_  At   )J3  -k,  At_  c0  2]Aa-m,u,  =  0 
nA   5A    6A   5A   2A      6A   2A  A  A 
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b)    Y      Aa-frY,,    A3+Y      Au    H-Y.    u+Y      Av    +Y.    V^+Y'    ij> 
SA  uA      A      uA   A      vA      A      vA  A      */A 


A 


+Yv  $A+Y   Ik   +  (k   -k   At..  )J5  -k   At  £  2]Ai^ 
^A  A   6A   1A    2A   XA   ^      2A   XA       A 


+Y   [k3  +(k   -k   At   )J3  -k   At  cO  2]A3 
A    A    *a    A    A       A    A 


+Y   [k,  +  (kr  ~kc  At0  )J3  -kc  At„  cO  2]Aa 
nA   5A    6A   5A   2A      6A   2A 


■mA[VUoVXg    ^A]=0 
yA 


f- 


c)    N      Aa+N,    A3+N      Au„+N.    u,+N      Av^+N.    ^,+Nr    iK 
aA  ^A    7      UA      A      UA   A      VA      A      VA   A       V  A 


+NiPAV'N6 Jki    +  <k2    "kl    Atl    ^~k2    *t    <02]A* 
VA  A       ■la  A         A       XA  A         A 


+N       [k      +(k       -k       At,     )£f  -k.    Atn    <#2]A3 
°A       JA  '\       3A       lA  4A       XA 


+N       [k      +(k       -k_    At0    )«#-kc    Atot#2]Aa 
nA       5A  6A       5A       2A  6A       2A 


"^/a-Y./VVa^0 


d) 
e) 
f) 


Precisely  Similar  Equations  for  Ship  B. 


g)    a  =  AVAu?\  =  «&' 


a 


h)    3   -  uoAi|-B-uoAiPA+AvB-AvA  =    #3 
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X.   Identifying  Zero-Value  Derivatives 

Arguing  from  symmetry  of  the  ship  about  its  centerline 

plane,  Abkowitz  demonstrates  that  X  ,  the  change  in 

x-direction  forces  of  the  ship  caused  by  a  change  in  v, 

with  all  other  variables  fixed  at  their  equilibrium  values, 

must  have  a  value  of  zero.   Similarly,  it  is  demonstrated 

that  X.  =  X'    =   Xv  =  0.   ([A-3],  pp.  I~28  and  1-30).   Due  to 
v     if;     \p  ^^ 

the  port  and  starboard  symmetry  of  the  rudder,  X.  is  also 

zero.   ([A-3],  pg .  1-67). 

Further  arguments  from  symmetry  demonstrate  that  in  the 

open  water  case,  Y  ,  Y. ,  N   and  N.  are  all  zero.   But  for 
r  u    u    u      u 

the  two-ship  system  of  interest,  the  interaction  side  force 

and  yaw  moment  are  affected  by  changes  in  speed.   Thus  Y   and 

N   cannot  be  assumed  to  be  zero.   The  effect  on  this  force 
u  — 

and  moment,  of  a  small  acceleration,  is,  however,  still  assumed 

negligible,  and  Y.  and  N.  are  assumed  to  be  zero. 
^   3  ■   '       u      u 

On  examining  the  hydrodynamic  derivatives  due  to  the 
orientation  of  the  ships  with  respect  to  each  other,  further 
identification  of  negligible  derivatives  can  be  made.   The 
derivatives  X   and  XR,  representing,  respectively,  the 
change  in  x-direction  force  caused  by  small  changes  in 
longitudinal  and  lateral  separations  of  the  ships,  will,  at 
most,  be  small.   Experiments  have  shown  that  interaction 
phenomena,  while  affecting  Y  forces  and  N  moments,  have  no 
measurable  effect  on  the  forward  motion  of  the  ships,  or  on 
their  resistance  forces.   ([A-l],  [A~2])   Thus,  X   and  XD  are 

Ci  p 

set  equal  to  zero. 
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Summarizing : 

X      =   X.    =   X-    =   X.,   =   X.  =  X      =   XD    -   Y.    =   N.    =    0  (2-13) 

vv  i|j  y  o  a  p  u  u 

With  this  information,  the  equations  of  motion  are  simplified 


to 


Ship  A: 


(2-14) 


a)  X   Au+X.  u^+X   [kc  +  (k   -k_  At_  )£ 

uA   A   uA  A   nA   5A    6A   5A  2^ 


•kc  At„  3  2]Aa-mAuA  -  0 
6A   2A         A  A 


b)  Y   Aa+YD  Ag-i-Y   Au„+Y   Av^+Y.  v\+Y;  iji 
aA     _3A     uA   A   vA   A   vA  A   ^A 


*/     x.  y 


2-,  ■ 


+Y-  $A+Y   [k   +(k   -k   At   )«&  -k   At  S      ]A^ 

^A  A   6A   1A    2A   XA   XA      2A   XA       A 


-l-Y   [k   +(k   -k   At   )J3  -k   At   fl  ^]A3 
A   JA    ^A   °A   XA      ^A   -""A 


+Y^  [k,  +(kr  -k_  At0  )£}  -k.    At   c82]A 
nA   5A    6A   5A   2A      6A   2A 


a 


■VVVA+XgAV     =    ° 
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c)    N      Aa+N       iB+Nu    4UA+N       AvA.!-N      V^-N-    *A 

/\  /i  *~i.  *"*  *x  n 


+NT|    $   +N.     [k      +(k9    -k      At      )JQT  -k      At     X3     ]  Ai[i 

^A  A         A  A         A       XA  A         A 


+N-    [k-   +(k,    -k-    At1    )£   -k.    At.,   c0  2]A3 

6A       ^A  4A       JA       XA  4A       ±A 


+N       [kc    +(k^    -kc    At„    )JD  -k-    At~    <02]Aa 
nA       5A  6A       5A       2A  6A       2A 


d) 

e) 
f) 


Precisely  Similar  Equations  for  Ship  B 


g)    J3  a   =    Aug    --    Au 


h'^=    Uo&*B    "    V*A    +    AVB    "    AVA 


XI.   Rearranging . Equations  of  Motion 

At  this  point,  terms  are  gathered  so  as  to  have  each 
variable  multiplied  by  the  sum  of  all  its  coeffieients  in 
each  equation. 

Remembering  that  cU   (variable)  -  3 (variable) /3t ,  then 

3u 


r=\, .  .    3  (Au)    3  (u-u  )     3u 


o 


(2-15) 


3t 


But  u   is  constant,  so  jj(Au)  =  ^Qu 
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Thus  : 

Ship   A: 

a)     [X      k.    +X       (ka    -k,.    At0    )Q  -X      ka    AtoP02]Aa 
nA    5A       nA       6A       5A       2A  nA    6A       2A 


(2-16) 


+  [X   +  (X.  -m  )c0  ]Au.  =  0 
UA    UA   A       A 


b)  [ (Y   +Y   k_  )+Y   (k,  -kc  At„  )£   -Y  kr    At„  &Z]ka 
aA   nA  5A    nA   6A   5A   2A      nA  6A   2A 


+  [(YR    +Y-    k-    )+¥.     (k.    -k_    At,     )£   -Yn    k.    At,  3    2]A| 

3A       6A    3A  6A       4A       3A       1A  UA    4A       ^ 


+  C»U    ]AuA+[Y      +(Y^A)(e]AvA 


+  [Y      k      +((Y.    -mAuo)+Y       (k2    -k      At      >  )£ 


A      A  rA 


A         AAA 


+  ((Y«    --mAXa J-Y6.k2,Atl    >5     )A,"a   "    ° 


A  JA  A      A         A 


c)      [(N      +N      k       )+N       (k      -k      At       )JQ  -N      k       At      J0^]Aa 

aA         A    DA  A       °A       °A       ^A  nA    °A         A 


+  [(NR    +N.    k-    )+N.     (k.    -k-    At,     )£)   -N.    k.    At,    ,5      ]AS 

3A       6A    3A  6A       4A       3A       1A  6A    4A       1A 


+N      Au+[N      +(N.    -m_x       )J0    ]Av^+[Nx    k,     +  ((Nr     -m_X      u    ) 
UA      A         VA         VA      A   9A  A         6A    1A  ^      A   gA    o 


+N^     (k_    -k,    At.     ))3   +((NTi    -I      )-N.    k_    At,     )o0   2]Ai|;      -    0 
6A       2A       XA      1A  *A       ZA         6A    2A      XA  A 


4  2 


In  writing  the  equations  for  Ship  B  use  is  made  of  the 
relationships  provided  by  equation  g  and  h  of  (2-14),  to 
substitute  for  Au.,  and  Av_.  in  the  Ship  B  equations.   From 

D  D 

these  equations,  Au_.  =  Au   +J9"a  and 

Av^  =  £f  3  +  Av„  +  u  AiIj,,  -  u  Ail;-.   Thus: 
B  A     o  ^A     o  rB 


Ship  B: 


(2-16) 


d)  [X   k   +(Xn  On      -k   At   )+X   )JQ 

nB  5B    nB   6B   5B   2B    UB 


+  ((X   -mB)-X   k   At   )£  ]Aa 


+  [X   +  (X.  -m  )J9  ]Au_  =  0 
uB    uB   B       A 


e)  [(Y   +Y   k_  )  +  (Y   (kr  -k_  At„  )+Y  \&    -Y   kr  At„  c&2]Aa 
aB   nB  5B     nB   6B   5B   2B    UB      DB  6B   2B 


+  [(Y   +Y   k_  )  +  (YR     (k.  -k-  Atn  )+Y„  )Jfr 
3B   6B  3B     6B   4B   3B   1B    VB 


+  ((Y.    -m   )-Y.    k.    At      )e8"2]A3+[Y      ]Au.  +  [Y„   +  (Y .    -m    )  &  ]  Av 

VB      B         6B    4B       1B  UB         A         VB         VB      B  A 


+  [UoV+Uo(Yv    -"B^^^ 


-I-    [(Y      k      -u  Y      )  +  ((Y-    -m_u   )+Y-     (k0    -k_    At.     ) 

6B    1B      o   vB  ^B      B   o        6B      2B      1B      1B 


-u    (Y      -m    ))£  +((Y      -m   x      )  -y      k      At..     )£f^]^      =    0 
0      vB      B  ^      B    gB         6B    2B      1B  *B 
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f)  [  (N   -|-N   k   )  +  (N   (kc  -kc  At0  )+N   )£  -N   kc  At.  d9  ]Act 
aB   nB  5B     nB   6B   5B   2B    UB      RB  6B   2B 


+  [(NR  +N   k_  )  +  (N.  (k.  -k-  At   )+N.  )J9+((N.  -mRx   ) 
6B   6B  3B     6B   4B   3B   3-B    VB         VB   B  gB 


■N   k   At   ))C02]A3+[N   ]Au  +[N   +(N.  -m  x   )J3  ]  Av 
6B  4B   1B  UB    A    VB    VB   B  gB       A 


+  [u  N   +U  (N.  -m_x  )ff  ]  Aip_  +  [(N.  k.  -u  N   ) 
°  VB   °   VB   B  gB       A    6B  XB   °  VB 


+  ((N*  -m„x   u  J+Njf  (k„  -k.  At,  ) -u  (N.  -m_x   )  )<£ 
^B   B  gB  °    6B   2B   XB   XB    °   VB   B  gB 


+  ((N..  -I   )-N   k   At   )J92]A^   =  0 
^B   ZB    6B  2B   XB        B 


These,  then,  are  the  six  equations  of  motion  in  six 
variables,  in  their  final  form  for  solution.   The  equations 
are  of  the  form:  (2-17) 

(AXl)Aa  +0-1-  (AXSJAu,.  +     0      +0      +     0     =0 
(AYl)Aa  +  (AY2)A3  +  (AY3)Au   +  (AY4)Av   +  (AY5)AiK  +     0     =0 

r\  x"i  r\ 

(ANl)Aa    +     (AN2)A3    +     (AN3)Au      +     (AN4)  Av      +     (AN5)Ai|j       +  0  -0 

r\  s\  /\ 

(BXl)Aa   +  0  -:-    (BX3)Au      +0  +0  +  0  =0 

(BYl)Aa    +     (BY2)A$    +     (BY3)Au      +     (BY4)Av_    +     (BY5)A^a    +     (BY6)A^    =0 

r\  A  A  tj 

(BNl)Aa    +     (BN2)A3    +     (BN3)Au,    +     (BN4)Av„    +     (BN5)A^    +     (BN6)A^    =0 

A  A  A  d 

This  is  a  set  of  six  simultaneous,  linear  differential 
equations  in  six  unknowns.   The  term  BX3 ,  for  instance,  is 
used  to  denote:   Ship  B,  X-equation,  3rd  term,  (or  Au   term). 
Each  of  the  coefficients  is,  itself,  a  polynomial  in  £)     ,  as 
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can  be  seen  from  equations  (2-16) ,  for  which  (2-17)  is  a 
shorthand  notation. 

XII.  Solving  the  Equations  of  Motion 

If  the  terras  AX1 ,  BX1 ,  etc.,  were  just  numerical 
coefficients,  straight-forward  algebraic  techniques  could 
be  used  to  solve  this  system  of  equations.   Yet,  while  these 
coefficients  are,  themselves,  polynominals  in  the  operator  i? 
Abkowitz  shows  that  <£)     can  be  treated  as  an  algebraic 
quantity.   ( [A-3] ,  pp.  1-26,  27).   This  is  possible  because 
the  hydrodynamic  derivatives  are  defined  as  the  slope  of  a 
force  or  moment  versus  a  dynamic  variable  taken  at  the 
equilibrium  condition.   Thus  the  terms  in  the  coefficients 
AX1,  etc.,  other  than  cD  ,  are  constants  in  time. 

The  solution  for  any  of  the  variables  is,  then,  the 
determinant  of  the  coefficients,  with  the  column  corresponding 
to  the  variable  being  solved  for  replaced  by  the  right-hand 
side  column  vector,  divided  by  the  determin?mt  of  all  the 
coefficients . 

For  example: 

(see  next  page) 
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Aa  = 


0   0  AX3   0  0    0 

0  AY2  AY 3  AY 4  AY5   0 

0  AN  2  AN  3  AN  4  AN  5   0 

0   0  BX3   0  0    0 

0  BY2  BY3  BY 4  BY5  BY6 

0  BN2  BN3  BN4  BN5  BN6 

AX1  0  AX3   0  0    0 

AY1AY2  AY3  AY4  AY5   0 

AN1AN2  AN3  AN4  AN5   0 

BX1  0  BX3   0  0    0 

BY1  BY2  BY3  BY4  BY5  BY6 

BN1  BN2  BN3  BN4  BN5  BN6 


DET 


(2-18) 


Evaluation  of  DET  leads  to  a  polynominal  in  3u      ,  (made 
up" of  sums  of  products  of  polynominals  in  JJ     ,  such  as  AX1, 
etc.)/  whose  coefficients  are  various  combination  of  the 
hydrodynamic  derivatives  and  control  system  constants.   Thus 
the  solution  for  each  variable  would  be: 


variable 


0 


,8  n  -l-  a.c& 


n-1 


•!-a 


n+1 


(2-19) 


If  DET  is  other  than  zero,  the  solution  for  all  variables 
would  be  identically  zero  for  all  time  after  a  small  disturbance 

from  equilibrium a  physical  impossibility.   Thus,  the 

denominator,  the  determinant  DET,  must  be  set  equal  to  zero. 
Setting  this  polynominal  in  <fi  ,  resulting  from  the  evaluation 
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of  DET,  equal  to  zero,  the  roots  of  the  polynomial  can  be 
obtained,  where  a  ,  o2,....o      are  these  roots.   Then: 


variable  = 


ur7^  (<a  "Vi* (Jb  "_ai)  (2~20) 


Abkowitz  demonstrates  how  the  operation  -t-k— -,  [z] 
is  equivalent  to: 

[z]  =  eat/e~at[z]dt  (2-21) 


(cQ-a) 


mT.     0       at  r  -at,-.  ,.       at 
Thus  ,-,&  _  -x  =  e    /e    (0)dt  =  c,e 


Carrying  out  this  operation  successively  to  evaluate 
equation  (2-20)  yields: 

.  !  t        a,  t  ,     a_t  a  t 

variable  -  c,e  1   +  c„e  2   c  e  n 

12  n 

where  the  c.  are  arbitrary  constants  of  integration. 

XIIl.  Establishing  the  Stability  Criterion  f or  the  System 

Since  each  of  the  variables  is  a  sum  of  exponential 

terms  in  time,  it  is  clear  that  the  only  way  the  variables 

can  all  go  to  zero  as  time  increases,  is  for  all  the  a.  (or 

l 

all  real  parts  where  the  a.  are  complex)  to  be  negative.   The 
criterion  of  stability  of  the  two  ship  system  is  for  all  the 
variables,  representing,  as  they  do,  deviations  from  the 
equilibrium  condition,  to  go  to  zero  as  time  increases.   For 
when  all  the  variables  have  gone  to  zero,  the  system  has 
regained  the  equilibrium  condition,  meeting  the  definition 
of  a  stable  system. 
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Thus,  the  stability  of  the  system  requires  that  the  real 
parts  of  all  <r .  be  negative.   As  the  cr*.  are  the  roots  of  the 
characteristic  equation  of  stability  resulting  from  setting 
the  polynomial  in  dJ   ,  (which  in  turn  results  from  evaluating 
DET) ,  equal  to  zero,  the  stability  requirement  for  the  two-ship 
system  is  that  the  real  parts  of  all  the  roots  of  DET  be  less 
than  zero. 


Chapter  III 
Procedure 

I.  Background  * 

As  stated  at  the  conclusion  of  Chapter  I,  a  modified 
OLNA/MARINER  class  merchant  hull  is  the  ship  used  in  the 
mathematical  model.   As  indicated  there,  this  was  done  because 
open  water  hydrodynamic  derivatives,  in  non-dimensional  form, 
were  available  for  a  MARINER,  (  [A-4]  ,  [A--5]  )  ,  and  because  the 
only  good,  quantitative  interaction  phenomena  data  available 
were  for  a  ship  quite  similar  to  a  MARINER the  R.F.A.  OLNA. 

II.  Open  Water  H ydr o d ynam ic  Derivatives 

Reference  A--5  presents  experimentally  determined  values 

for  most  open  water  derivatives,  measured  with  a  planar  motion 

mechanism,  for  a  MARINER  model.   Reference  A-4  utilized  these 

values,  with  the  exception  of  X  ,  for  which  a  different  value 

i  u 

was  determined.   The  accuracy  of  the  values  used  in  [A-4]  was 
borne  out  with  full  scale  sea  trials  and  these  values  are  used 
in  this  thesis  for  the  open  water  derivatives,  for  ship  A,  and 
are  listed  in  Table  III-l. 

III.  Obtaining  the  In t e raction  Hydrodynamic  Derivatives 

A.   Y  ,  YQ,  N  ,  ND. 
a        3   a    3 

Values  for  hydrodynamic  derivatives,  as  such,  arising 
from  interaction  effects,  could  not  be  obtained.   The  inter- 
action data  for  OLNA  consisted  of  dimensionless  forces  and 
moments  at  side-to-side  ship  separations  of  50  and  100  feet, 
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for  all  longitudinal  positions  of  the  midships  sections  of 
the  ships  between  -600  and  +600  feet.   Also  available  were 
force  and  moment  data  for  OLNA  at  side-to-side  ship  separations 
of  from  30  to  150  feet,  with  the  ships  abreast.   [A-l]   (See 
Figures  1-3,  1-4  and  1-5) 

Using  the  data  illustrated  by  Figures  1-3,  4  and  5,  the 
only  positions  at  which  Y  ,  Y n ,  N   and  N„  could  all  be 
determined  were  at  side-to-side  distances  of  50  and  100  feet, 
with  the  ships  abreast.   (Points  1  and  6,  Figure  III-l) .   This 
limitation  of  the  examination  of  system  stability  to  only  two 
ship  positions  was  felt  to  be  too  restrictive.   As  the  only 
means  available  for  obtaining  the  desired  derivatives  at  other 
positions  of  interest  was  extrapolation  and  interpolation  of 
the  OLNA  data,  this  course  was  followed  for  both  Y  force  and 
N  moment.   As  data  was  available  at  side-to-side  distances 
(STSD,  for  short),  of  50  and  100  feet,  these  values  were  chosen 
as  bounds  to  the  range  of  interest,  just  as  values  oi:  -600  and 
+600  feet  were  chosen  as  bounds  in  the  longitudinal  direction. 
The  way  in  which  Y  force  was  interpolated  in  the  range  of 
interest  will  be  illustrated. 

The  Y  force  was  known  along  three  lines  of  position,  as 
shown  in  Figures  J. -3,  and  1-5,  and  as  in  Figure  III-l.   To 
determine  Y  at  positions  other  than  those  known  in  Figure  III-l, 
it  was  assumed  that  the  ratio  of  change  in  Y  between  1'  and  2' 
to  the  change  between  1'  and  6'  was  the  same  as  the  ratio  of 
the  change  in  Y  between  1  and  2  to  the  change  between  1  and  6. 
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Or: 


(Y2,-Yr) 


(W 


(Y3,,-Y1.) 


(Yg.-Y  ,) 


(Y6-Yx) 


(Yg.-Y^) 


<VY1> 

(Y6-yi> 


etc 


(3-1) 


The  points  1  through  6  in  Figure  III-l  were  chosen  at 
10  foot  intervals,  and  the  interpolation  for  these  values  of 
STSD-60,  70,  80  and  90  feet  was  performed  at  50  foot  intervals 
in  a.   Thus,  values  of  Y  were  obtained  at  all  intersections  of 
a  grid  with  lateral  spacing  of  10  feet  and  longitudinal 
spacing  of  50  feet,  as  shown  below: 


-I 


1  _ 


-100 


-50' 


I 
I 

o' 


I 


I 

I 
I 

so' 


STSD-  50 
srso  =  c-o' 

ST5D-7o' 

srsor  so' 


loo- 


Figure  III-2 

The  results  of  such  interpolations  for  both  Y  and  N  are 
shown  as  figures  III-3  and  III-4,  where  the  curves  for  STSD  = 
50  and  100  feet  are  from  [A--1]  and  the  others  are  interpolated 
The  interpolated  curves  are  felt  to  give  a  reasonable  picture 
of  the  way  the  interaction  force  and  moment  would  be  expected 
to  vary  as  STSD  is  changed.   And,  of  course,  the  points  on  the 
line  a=0  are  the  same  points  as  those  plotted  in  Figure  1-5. 
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NokJ-  Dwensional    Y  Force     vs.  ex 
at    Side- to- Side:    Distances 

OF    50,  GO,  70,80,   90    AND    100     FT. 


FORCE      OKI      SHIP    A       A5 
SHIP    B     PASSES 


a 


Dimensionless  Y  Force  vs.  Longitudinal  Separation 

FIGURE  II 1-3 
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Dimensionless  N  Moment  vs.  Longitudinal  Separation 

FIGURE  I II -4 
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Obtaining  the  interaction  forces  and  moments  over  this 

range  of  50'  <  STSD  <  100'  and  -600'  <    a    <    +600'  is,  though, 

just  a  step  along  the  way  to  obtaining  the  derivatives  Y  ,  Yft, 

N   and  N„  ,  vzhich  are  what  is  really  desired, 
a       3 

The  longitudinal  separation  parameter,  a,  is  represented 

directly  on  the  curves  of  Y  and  N.   The  lateral  separation 

parameter,  3/  is  related  to  the  side-to-side  distances  shown 

on  the  figures  by  a  constant,  for  a  given  pair  of  ships. 

Since  3  is  the  lateral  separation  of  the  centerlines  of  the 

ships,  it  is  apparent  that: 

3  =  STSD  +  (BEAM-,  .   _  +  BEAM_,  .   „)/2 

Ship  A        Ship  B  ' 

=  STSD  +  Constant  (3-2) 

Thus,  for  two  ships  of  50'  beam,  an  STSD  of  50'  is 

equivalent  to  3=100'.   Therefore,  Figure  III-3,  for  example, 

represents   a  surface  in  space,  (as  sketched  in  Figure  III-l) , 

whose  horizontal  coordinates  are  a  and  3  and  whose  vertical 

coordinate  is  the  value  of  Y  everywhere  on  the  surface.   Then, 

to  get  3Y/8a,  or  Y  ,  at  any  point,  the  slope  of  the  Y  surface 

in  the  a  direction,  at  the  point,  must  be  obtained.   Similarly, 

YD    is  the  slope  in  the  3  direction.   For  example: 
p 

a)  Y(a,3)  =  Y(100,60) 

-  (  (Y (150, 60) ~Y (50,60)  )/ (150-50) 

b)  Y(a,3)  =  Y(100,60) 

-  ((y  (100,70)-Y(100,50) )/(70-50)  (3-3) 


In  this  way,  Y   everywhere  in  the  range  (50 ' <STSD<100 ' ; 
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-550' <  a  <  550'),  V7as  obtained,  as  was  Y   everywhere  in  the 

range  (60  '  <STSD_<90  '  ;  -600 '<  a  <    600').   (The  values  cf  Y   and 

YD  at  the  extremities  of  the  range  of  known  Y  values,  in  the 
p 

direction  of  differentiation,  could  not,  of  course,  be  obtained 
by  this  simple  procedure) . 

The  values  of  N   and  N0  over  the  same  ranges  were  obtained 

a      3 

in  the  same  way. 

A  computer  program,  program  INTCOF ,  was  written  to  take 
the  values  of  Y  and  N  at  the  locations  known  from  Figures  1-3, 
4  and  5,  to  interpolate  to  get  Y  and  N  throughout  the  range, 

and  then  determine  Y  ,  Yn ,  N   and  N„  in  this  range.   The  curves 

a    3    a      3 

in  Figures  III-3  and  III-4  for  side-to -side  distances  of  60, 
70,  80  and  90  feet  are  plots  of  the  Y  and  N  outputs  of  the 
INTCOF  interpolations. 

Similarly,  the  curves  in  Figures  III-5  and  III-6  for 
STSD  60,  70,  80  and  90  feet  are  INTCOF  outputs.   And  in 
Figures  III-7  and  III--8,  the  curves  are  all  INTCOF  outputs. 

Program  INTCOF  is  described  in  detail  in  Appendix  II, 
and  the  results  of  the  program  for  the  modified  OLNA/MARINUR 
ship  are  included  therein. 


B.   Y  ,  N  . 
u    u 


As  noted  in  Chapter  I,  experimental  investigators  in  the 
field  of  ship  interactions  have  noted  that  interaction 
phenomena  vary  with  speed  in  the  same  way  as  the  resistance  of 

the  ship or  that  interaction  effects  vary  as  the  square  of 

the  speed.   The  OLNA  data  used  for  interaction  effects  was 
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obtained  at  a  speed  of  15  knots,  which  was  also  used  as  the 

equilibrium  speed  for  the  results  presented  in  this  thesis. 

(However,  all  force,  moment  and  derivative  values  are 

presented  in  non-dimensional  form  and  are  valid  at  other 

speeds. ) 

Remembering  that  Y  ,  for  example,  is  3Y/3u  with  all  other 

variables  at  equilibrium  values,  and  knowing  Y  for  a  given 

position  at  15  knots,  it  is  a  simple  matter  to  get  Y  at  14  and 

16  knots,  and  then  Y   at  15  knots. 

u 

a)  Y(16)  =  Y(15)  x  (16/15) 2 

b)  Y(14)  =  Y(15)  x  (14/15) 2 

c)  Yu(15)  =  (Y(16)-Y(14))/((16-14)x(1.6878))  (3-4) 

N   is  obtained  by  a  similar  process. 

C.   X  ,  Y  ,  N  . 
n  _  n    n 

The  effects  on  forces  and  moments  of  a  change  in  piropeller 

speed  are  also  required,  and  are  represented  by  X  ,  Y  ,  and  N  . 

xi  xr-  J   n   n       n 

The  derivative  X   is  determined  by  using  the  known  values 
of  Effective  Horsepower,  (EHP) ,  and  RPM  for  a  MARINER  at 
speeds  around  15  knots. 

a)  X(14kts)  =  (EKP(14)x550)/(15xl.6878) 

b)  X(16)  =  (EHP(16)x550)/(16xl.6878) 

c)  X  (15)  =  (X(16)-X(14)  )/(((RPM(16)-RPM(14))/60)      (3-5) 

The  value  of  X   so  obtained  is  in  dimensional  form,  and 
n 

is  then  non-dimensionalized  to  permit  broader  applicability. 
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The  single  screw  MARINER  requires  a  1.2°  port  rudder 

angle  to  maintain  a  steady  course  in  open  water.   [A-4]   Since 

the  value  of  Yr    for  the  MARINER  is  known  from  reference  A-4, 
o 

it  is  possible  to  determine  what  side  force  is  being  generated 

by  this  1.2°  rudder  angle,  which  is  denoted  by  (6  )   .   When 
2  3  J         o  ow 

maintaining  a  straight  course  in  open  water ,  this  rudder  side 

force  is  offsetting  the  side  force  caused  by  the  asymmetry  of 

a  single  screw.   Thus,  on  a  steady  course,  Y   ,,   =-Y 

J  •*  rudder    prop 

The  dimensionless  Y~  given  in  [A-4]  can  be  converted  to 

o 

its  dimensional  form,  which  is  in  units  of  lbs/radian.   The 

actual  Y  force  exerted  can  then  be  calculated  by  multiplying 

the  dimensional  Yr  by  (6  )   /57.3°,  giving  a  side  force  in 

6  •*         o  ow        ^     3 

pounds  which  is  the  negative  of  the  force  being  caused  by  the 
propeller.   The  force  exerted  by  the  propeller  can  be  expected 
to  vary  as  the  propeller  speed  squared,  and,  thus,  knowing  the 
propeller  side  force  at  a  given  speed,  it  is  possible  to 
determine  how  this  force  varies  as  propeller  speed,  n,  varies. 

a)  Y   ..,    =  Y~xNDF(@15  kts)  x  (6  )   /(l  radian) 

rudder     6  o  ow 

b)  Y     (15  kts)  =  -Y   ,-.   (15) 

prop  rudder 

c)  Y     (14)  =  Y     (15)  x  (n(14)/n(15) )2 

prop        prop  ' '. 

d)  Y     (16)  -  Y     (15)  x  (n(16)/n(15) )2 

prop        prop  - ' 

e)  Y  (15)  -  (Y     (16)-Y     (14) ) /(n (16) -n (14) ) 

n         Prop      prop     ' 

(Where  NDF  is  the  non-dimensionalizing  factor  for  Y~) 

The  value  of  N   is  similarly  obtained,  and  the  Y   and 

n  J  n 

N   so  gotten  are  suitably  non--dimensionalized. 
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A  computer  program,  program  DETERM,  was  written  to 

determine  the  values  of  Y  ,  N  ,  X  ,  Y   and  N  .   The  inputs 

u    u    n    n       n  L 

for  the  program  are  the  values  of  Y  and  N  at  the  positions  of 

interest,  gotten  from  program  INTCOF,  and  the  EHP ,  RPM,  speed 

and  (6  )    values  for  the  ship  being  examined.   The  values  of 
o  ow  c 

Y~  and  N~  in  open  water  are  also  inputs.   Program  DETERM  is 
described  in  detail  in  Appendix  III. 

IV.   Calculation  of  Interaction  Derivatives  for  Ship  3. 
All  discussion  concerning  interaction  hydrodynamic 
derivatives,  so  far,  has  been  for  the  OLN A/MARINER  ship  which 
is  assumed  to  be  Ship  A  of  the  two  ship  system.   Ship  A  serves 
as  the  reference  for  measuring  both  longitudinal  and  lateral 
separations,  a    and  3,  and  is  assumed  to  be  the  overtaken  ship, 
located  to  the  left  of  Ship  B,  the  overtaking  ship. 


CJZ> 


This  coincides  with  the  case  for  the  OLNA  experiments  and  all 
plots  of  the  interaction  variables  vs.  a  are  for  the  overtaken 
ship  on  the  left. 

But  hydrodynamic  derivatives  for  both  ships,  at  any 
relative  position  of  interest,  are  needed  to  write  the  equa- 
tions of  motion  for  the  system.   Since  the  modified  OLNA/ 
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MARINER  used  as  Ship  A  is  the  only  one  for  which  it  was  possible 
to  obtain  all  required  derivatives,  Ship  B  must  also  be 
considered  to  be  this  type  of  ship,  in  order  to  be  able  to 
obtain  derivatives  for  it.   If  the  two  ships  are  identical, 
the  force  or  moment  felt  by  Ship  A  at  a  given  a  and  3  are  the 
negatives  of  the  force  or  moment  felt  by  Ship  £  at  the  same 
p  a±  ol'-o..   To  illustrate: 

1.  2. 


v* 


Y  on  Ship  B  at  2  equals  -Y  on  Ship  A  at  1 


>hip  A  is 


attracted  to  B  in  position  1,  then  B  will  be  attracted  to  A 
in  the  same  way  at  2.   But  because  of  the  reference  system, 
(as  illustrated  in  Figure  II-2) ,  an  attraction  for  Ship  A  is 
a  positive  Y  force,  while  for  Ship  B  an  attraction  is  a 
negative  Y  force. 

Thus,  in  program  DETERM,  the  Y  and  N  values  entered  for 
Ship  B  are  the  negatives  of  those  for  Ship  A  when  a~-a.   As 

an  example,  with  STSD  =  70'  and  a  =  300',  the  non-dimensional 

-3 

Y  for  Ship  A  is  .1652  x  10   .   The  Y  value  entered  for  Ship  B 

-3 
in  this  case  is  --.16  89  x  10   ,  which  would  be  the  value  for 

Ship  A  when  a  =  -300'.   (These  values  are  from  the  INTCOF 
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output  data  of  Appendix  II.)   With  the  inputs  so  adjusted, 

program  DETKRM  will  calculate  correct  values  of  Y   and  N   for 
^  u      u 

both  ships.   (The  calculations  for  X  ,  Y  .  and  N   are  indepen- 

n   n       n 

dent  of  those  considerations) .   The  output  of  program  DETERM 

is,  therefore,  X  ,  Y  ,  and  N   for  Ship  A,  (which  are  assumed 

n    n       n        c 

to  be  the  Name  for  Ship  B) ,  Y   and  N   for  Ship  A  and  Y   and  N 

r      u      u        l  u      u 

for  Ship  A,  if  it  were  in  Ship  B's  position. 

To  bypass  the  necessity  for  considering  both  ships 
identical,  provision  was  made  to  let  Ship  B  be  of  a  length 
different  from  Ship  A,  but  of  the  same  form.   All  form 
coefficients,  plus  the  length/beam  and  beam/draft  ratios,  are 
assumed  tho  same  for  both  ships.   This  assumption  implies 
that  the  non-dimensional  values  of  the  hydrodynamic  derivatives 
for  Ship  B  are  the  same  as  those  for  Ship  A,  (with  some  sign 

difference:.,  as  in  Y„  and  N„),  when  all  derivatives  are  non- 

p       p 

dimensional  i  zed  with  respect  to  the  length  of  the  ship  to 
which  they  belong.   But  for  the  two  ship  system,  the  charact- 
eristic length  used  in  the  non-dimensionalizing  process  must 
be  the  samo  for  aH  derivatives.   It  was  decided  to  use  the 
length  between  perpendiculars,  LBP ,  of  Ship  A.   Thus,  if  Ship 
B  is  of  a  length  different  than  that  of  Ship  A,  the  dimension- 
less  derivatives  of  Ship  A  must  be  multiplied  by  (LBPB/L3PA) 
to  get  the  viimensionless  derivative  for  B  with  I'BP   replacing 
LBP   in  the  non-dirnensionalizing  factor. 

This  modification  of  the  Ship  A  derivatives  to  fit  Ship  B 
is  done  in  the  early  part  of  a  third  computer  program,  SOLVE. 
(This  program  is  used  to  actually  solve  the  eauafions  of 


motion  for  the  stability  roots,  as  outlined  below,  but  modifies 
the  derivatives,  as  stated,  before  solving  the  equations) . 
SOLVE  is  described  in  Appendix  iv. 

V.   The  Complete  List  of  the  ) I ydrodynamic  Derivatives 

Table  III-l  lists  the  dei  ivatives  from  reference  A-4  for 

the  MARINER  and  these  are  the  values  used  for  the  open  water 

derivatives  for  Ship  A,  and  wh  i.eh  are  modified  to  apply  to 

Ship  B  in  program  SOLVE. 

Table  III--2  lists  the  dei  ivatives  associated  with  propeller 

speed  which  were  calculated  by  program  DETERM  for  Ship  A.   In 

program  SOLVE,  they  are  used  f,,r  Ship  A  and  modified  for  Ship  B. 

The  values  of  Y  ,  Y  ,  N  ,  m  ,  Y   and  N   are  position- 
er   p    ex    3    u      u     ^ 

dependent,  since  the  values  of    y  and  N  acting  as  a  result  of 
interaction  effects  vary  as  the  ships  change  positions  relative 
to  each  other.   In  addition,  these  derivatives,  even  if  ships 
A  and  B  are  identical,  will  noi  be  the  same  for  the  two  ships. 
Table  III-3  gives  values  of  those  derivatives  for  the  two 
ships,  assuming  them  identical,  for  two  typical  positions  of 
the  ships  relative  to  each  othor. 
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TABLE    III-l 


Derivative           Non-Dimension- 

Non- Dimensional 

.  .    '.".  alizing  Factor 

.  .    Value 

(X.-m) 
u 

l/2pL3 

-840. xio"5 

X 
u 

l/2pL2UQ 

-120. xl0~5 

Y 

V 

l/2pL2u 

-1160.4xl0~5 

(Y.-m) 

V 

l/2pL3 

-1546. xl0~5 

(YfmuQ) 

l/2pL3u 

-499. xlO-5 

(YfmXg} 

1/2PL11 

8.6xl0~5 

l/2pL2u2 
o 

277.9xl0"5 

N 

V 

l/2pL3u2 
o 

-263.5xl0~5 

(N.-mx 
v    g 

)                1/2PL1* 

22.7xl0"5 

(N~mx  u 
r      g 

)               l/2pL4u 
o                     o 

-166. xl0~5 

(NfI2) 

l/2pL5 

-82.9xl0~5 

N6 

l/2PL3u2 

-138.8xl0~5 

P 

mass  density  of  sea  water   1 

.9905  lb-sec2/ft't 

L 

LBP  Ship  A,  in  ft 

u    = 
o 

equilibrium  speed  in  ft/sec 

TABLE  I I 1-2 

Derivative 

Non-Dim.  Factor 

Non-Dim.  Value 

X 

n 

Y 
n 

N 
n 

l/2pL3u 

l/2pL3UQ 

l/2pL\i 

4.62xl0"5 

-.52xl0"5 

.26xl0~5 
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TABLE  III-3 


tf=0';  STSD 

=  60' 

Deriv . 

Non-Dim.  Factor 

Ship  A 

Ship  B 

Non-Dim.  Value 

Non-Dim.  Value 

Y 
u 

l/2pL2UQ 

152.5xl0~5 

-152.5xl(T5 

N 

u 

l/2pL3uo 

-67.1xl0~5 

67.1xl0~5 

Y 

a 

l/2pL  u2 

.16xl0"5 

.16xl0~5 

Y8 

1/2 pL  u2 
K    o 

-.87xl0"5 

.87xl0"5 

N 

a 

l/2PL2u2 
o 

.126xl0~5 

.126xl0~5 

\ 

l/2pL2u2 
o 

.42xl0"5 

-.42xl0~5 

«=  -400' ; 

STSD=80' 

Y 
u 

Same 

-47.6xl0"5 

-4.29xl0~5 

N 
u 

17.3xl0~5 

-16.7xl0~5 

a  ) 

as 

.044xl0"5 

-.093xl0~5 

\ 

_^  .447xl0"5 

~.lllxl0~5 

Above  <- 

N 

-.068xl0~5 

-.035xl0~5 

N3 

-.14xl0"5 

.098xl0~5 

VI.   Method  of  Solving  for  Stabi lity  Roots 

In  Chapter  II,  it  was  shown  that  the  stability  charac- 
teristic equation  for  the  system  was  obtained  by  evaluating 
the  determinant  formed  by  the  coefficients  of  the  variables. 
These  coefficients  are  polynominals  in  D  ,  and  it  was  stated 
that  the  operator  could  be  treated  as  an  algebraic  quantity, 
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so  that  the  result  of  the  determinant  evaluation  was  another 
polynominal  infl  ,  whose  roots  were  the  stability  roots. 

The  coefficient  determinant  to  be  evaluated  is  a  6  x  6 
determinant,  given  as  equation  (2-17).   The  determinant  is 
evaluated  by  pivotal  condensation  ((A-lJ,  pg .  155,  problem 
4-21).   An  illustration  of  pivotal  condensation  follows: 


all  al2 


a21  a22 


a  t  a 
nl   n2 


In 


2n 


nn 


all* 


1    a 


11 
11 


a21  a22 


a  ,  a  „ 
nl   n2 


a 


In 

11 


2n 


nn 


:a  x 

11 


1    0 


21   2  2       2n 


a  -,  a  ~ 

nl   n2 


nn 


:alxX 


22 


n2 


a 


2n 


nn 


(3-7) 


Where  a !  .  =  a .  . -  (a,  . /a, , )  (a . , ) 
13     13    li   11    ll 


The  procedure  would  then  be  repeated  with 


a 


22 


n2 


"2n 


nn 


The  matrix  properties  whereby  the  interchanging  of  two 
rows  or  two  columns  changes  the  sign  of  the  determinant  are 
also  utilized.  After  three  applications  of  pivotal  conden- 
sation and  use  of  these  two  matrix  properties,  the  determinant 
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in  (2-17)  becomes: 

(AX1)  (BX31)  (BN6)  * 


AY4   AY5   AY2 
BY4'  BY5'  BY2' 
AN  4   AH  5   AN  2 


(3-8) 


Where 


BX31 

=  BX3  - 

BX1 (AX3/AX1) 

BY4' 

=  BY4  ~ 

BY6 (BN4/BN6) 

BY5' 

=  BY5  - 

BY6 (BN4/BN6) 

BY21 

=  BY2  ~ 

BY6 (BN4/BN6) 

(3-9) 

Carrying  the  evaluation  of  (3-8)  to  its  conclusion,  and 
utilizing  the  relationships  of  (3-9),  the  determinant  becomes: 

DET  =  [ (AX1) (BX3)-(BX1) (AX3) ] 

x{  [  (BN6)  (BY5)~ (BY6)  (BN5) ]X  [  (AY4)  (AN2)- (AN4)  (AY2) ] 
+  [(BN6)  (BY2)-(BY6)  (BN2) ]X  [  (AY5)  (AN4)-(AN5)  (AY4) ] 

+  [  (BN6)  (BY4)- (BY6)  (BN4) ]X  [  (AY2)  (AN5)- (AN2)  (AY5) ]  } 

(3-10) 

Where  the  various  terms  are  defined  in  Table  III-4. 

Program  SOLVE,  after  appropriately  modifying  all  hydro- 
dynamic  derivatives  to  apply  to  Ship  B  as  v/ell  as  Ship  A,  then 
calls  a  subroutine,  DETER,  to  carry  out  the  mathematical 
operations  to  evaluate  (3-10)  at  one  or  more  positions  and 
with  one  or  more  sets  of  control  system  constants  for  the  two 
ships  at  the  various  positions.   Each  of  the  terms  in  (3-10), 
AX2,  BN6,  AY2,  etc.,  is  a  polynomial  in  D   ,  and  the  various 
multiplications,  subtractions  and  additions  are  carried  out  in 
subroutine  DETER  by  using  IBM  provided  Scientific  Subroutine 
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Package  (SSP)  subroutines  for  arithmetic  operations  with  poly- 
nomials.  Program  SOLVE,  (and  all  the  subroutines  it  calls), 
is  described  in  detail  in  Appendix  IV. 


VII.   Control  Constant  Magnitude 

Reference  A--11  discusses  automatic  ship  control  and 
utilizes  some  typical  values  for  the  control  constants  in  a 
system  sensitive  to  ip  and  ty .       Sensitivity  to  an  error 
parameter,  such  as  ty,    is  proportional  control,  and  to  an 
error  parameter  rate,  such  as  ip,  is  rate  control.   The 
equation  for  rudder  control  in  such  a  system,  (not 
considering  signs),  is: 


6  -  k1ip  +k    ty  (3-11) 

In  [A-ll] ,  Schiff  and  Gimprich  use  values  of  k   between  1.0 

and  5.0,  as  typical  values.   The  only  values  used  for  k   are 
0  or  1.0. 


With  k 


,  =  1.0/  a  1°  heading  error  will  give  rise  to  1' 


of  correcting  rudder,  while  with  k,  =  5.0,  a  1°  error  will 
cause  a  5°  rudder  correction.   If  k„  =  1.0,  a  heading  error 
rate  of  l°/sec.  gives  a  rudder  angle  of  1°  for  correction. 
Thus,  k,  was  considered  to  vary  from  1.0  to  5.0,  (with  smaller 
values  investigated  to  determine  some  trends) ,  and  k„  was 
varied  between  0.5  and  4.   (Defining  6  and  ip  to  be  in  radian 
does  not  affect  k   and  k  ) . 

In  determining  ranges  to  be  used  for  k_  and  k.,  which 
determine   sensitivity  of  the  control  system  to  lateral 
separation  error  and  its  rate  of  change,  3  and  3/  no  basis 
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was  found  in  the  literature.   Thus,  these  constants  were 
estimated  from  what  the  author  expected  to  be  reasonable 
corrections  of  the  rudder  angle  for  lateral  separation 
errors.   Rudder  angle  is  sensitive  to  3  as  below,  (again 
signs  are  not  considered) : 

6  =  k33  +  k43  (3-12) 

Based  on  the  author ' s  experience  in  underway  replenishment 
operations,  it  was  assumed  that  an  error  of  10  feet  in  3  would 
prompt  a  rudder  correction  of  about  1° .   And  a  change  in  3  of 
1  ft/sec.  was  felt  to  call  for  3°  of  rudder  correction. 
Therefore : 

i  «'  'l'°     1  radian  ~  1   ~  n^^n 
3   10ft     57.3°     57.3 

1  radian  ~  3    ~nc„  ,_  ..  _. 

x  —■=■= — 5-5 ~-p^ — 5-  -.052  (3-13) 


4   lft/sec     57.3°     57.3 

As  a  result,  k_  was  investigated  between  .0005  and  .01 
and  k.  between  .01  and  .1. 

A  propeller  speed  control,  sensitive  to  longitudinal 

separation  error  and  its  rate,  a  and  a,  was  assumed  in  the 

theory.   For  this  sensitivity,  not  considering  signs: 

n  =  kea  +  k,a  (3-14) 

o      6 

To  get  an  idea  of  the  magnitude  of  k-  and  k,,  estimates 

were  resorted  to.   It  was  assumed  that  a  typical  error  of 

10  feet  in  a  would  prompt  an  n  change  of  1  RPM,  (about  1/4 

knot  for  a  MARINER)  and  a  1  ft/sec.  value  for  a   would  prompt 

a  change  of  4  RPM.   Therefore: 
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1_RPM        1 ~   1 

10  ft    6  0RPM/RPS    600' 


k 


4  RPM 


6    lft/sec 


-  x 


1 ~  1_ 

6  0RPM/RPS    15 


.0016 


-  .066 


(3-15) 


So  k_  was  assumed  to  be  between  .0005  and  .01  in  value, 

and  k..  between  .01  and  .2. 
6 

VIII.   Control  Constant  Signs 

When  either  ship  develops  a  positive  error  in  i(j  or  if/,  a 
positive  rudder  angle  is  required  to  correct  for  the  error. 
Thus  k,  and  k„  are  positive  for  both  ships. 


■<  ( 


+  y  j  +  H* 


V4-  Q 


+  S 


+y>,  +  y 


As  $  or  3  increase  in  a  positive  direction,  it  is  desired 
to  turn  Ship  B  to  the  left  and  Ship  A  to  the  right,  to 
compensate.   Thus  as  3  increases,  a  negative  5  is  required 
on  Ship  A,  but  a  positive  6  is  wanted  for  Ship  B.   Therefore 
k~  and  k.  are  negative  for  Ship  A  and  positive  for  Ship  B. 


6 


h 
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As  a  increases  (or  ~a  decreases  in  magnitude) ,  it  is 
desired  to  increase  Ship  A's  speed  and  decrease  Ship  B's. 
Therefore  k   and  k   are  positive  for  Ship  A  and  negative  for 
Ship  B. 

Both  rudder  and  propeller  response  will  have  time  lags 
associated  with  them.   The  rudder,  however,  can  be  expected 
to  respond  more  quickly  than  propeller  speed.   Thus  rudder 
time  lags  from  0  to  4  seconds,  and  propeller  time  lags  from 
0  to  6  seconds  are  examined.   (The  desire  is  to  determine  the 
way  in  which  time  lags  affect  stability,  and  errors  in  the 
actual  values  used  will  not  affect  trends  developed.) 

To  summarize,  Table  III-5  lists  all  the  control  system 
parameter  ranges  examined. 

TABLE  III-5 


Ship  A 

Ship 

B 

1 

1  kll  5. 

1 

< 

ki  < 

5. 

.5 

<  k2  <  4. 

.5 

< 

k2  1 

4. 

-.01 

<_   k3  <  -.0005 

.0005 

< 

k3  £ 

.01 

-.1 

1  k4  1  --01 

.01 

< 

k4  < 

.1 

.0005 

<  kc  <  .01 

—    D  — 

-.01 

< 

k5  < 

-.0005 

.01 

<  kr    <    .2 
—   b  — 

-.2 

< 

\± 

-.01 

0 

<At,  <   4. 

0 

< 

Atx  < 

4. 

0 

<_At   <   6. 

0 

< 

At2  1 

6. 
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Chapter  IV 
RESULTS,  CONCLUSIONS  AND  RECOMMENDATIONS 

I.  Introduction 

The  intent  of  this  thesis  is  to  develop  and  solve  the 
equations  of  motion  for  the  two  ship  system,  and  to  obtain 
the  stability  roots  for  it.   The  method,  with  the  computer 
programs  used  to  perform  the  calculations,  (which  would  be 
prohibitive  by  hand) ,  permits  determination  of  whether  or 
not  the  system  is  stable  for  a  given  set  of  control  constants 
and  time  lags  at  given  positions.   In  order  to  examine  the 
accuracy  of  the  solutions  and  correctness  of  the  method, 
various  sets  of  control  system  constants,  for  various  ship 
positions,  were  used  as  inputs  to  the  computer  programs  and 
the  effect  on  the  stability  of  the  system  observed. 

It  is  beyond  the  scope  of  this  thesis  to  attempt  an 
optimization  of,  or  even  complete  examination  of,  the  effects 
of  variations  in  the  automatic  control  parameters  upon 
system  stability.   The  purpose  is  to  present  the  method,  with 
only  enough  calculation  to  assess  the  value  of  the  theoretical 
formulation  and  computer  solution.   As  a  resu.lt,  the  various 
combinations  of  controls  and  ship  positions  examined  are  by 
no  means  exhaustive but  are  representative. 

II.  Resu Its 

A  few  sets  of  control  constants  were  evaluated  at  each  of 
four  positions,  for  ships  of  different  size,  with  and  without 
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time  lag,  and  are  presented  in  Appendix  IV  as  typical  output 
from  program  SOLVE . 

To  perform  a  somewhat  more  complete  examination  of  system 
stability,  a  single  position,  (a  =  0',  STSD  =  90'),  was 
chosen,  and  a  total  of  91  combinations  of  control  constants 
and  time  lags  was  run.   The  results  of  these  runs  are 
summarized  in  Tables  IV-1,  IV--2,  and  IV--3.   In  these  tables, 
the  number  of  zero  roots,  number  of  negative  roots  and  number 
of  positive  roots  obtained  for  the  input  control  constants 
are  indicated.   In  addition,  the  largest  root  for  each  condition 
is  listed,  unless  the  largest  root  is  zero.   In  this  case  the 
largest  negative  root  is  listed.   The  effects  of  changes  in 
some  of  the  control  constants  on  the  value  of  the  largest 
root  are  shown  in  Figures  IV- 1  through  IV-- 5. 

III.   Discussion  of  Results 

The  first  part  of  Table  IV-1  summarizes  the  effect  of 
varying  k, ,  with  all  other  control  parameters  zero  and  no 

time  lags.   These  results  are  illustrated  in  Figure  IV-1. 

k 
Varying   1  ,  alone,  from  a  very  low  value  to  its  maximum, 

decreases  instability  for  a  time,  (the  value  of  the  largest 

root  is  positive,  but  decreasing),  but  then  instability 

k  k       k 

increases  as   1,,  increases.   When   1,.  and   1.-.  are  both 

A  A         rs 

increased  through  this  same  range,  the  system  instability  is 
lessened,  with  stability  actually  acheived  for  part  of  the 
range,  but  the  trend  for  instability  to  increase  for  the 
higher  k,  values  is  still  present. 
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The  system  investigated  to  get  the  results  presented  in 
this  chapter  consisted  of  two  identical  ships.   This 
permitted  a  check  of  the  method  to  see  that  the  solution  is 

symmetrical  for  this  symmetrical  system.   As  a  check,  a  run 

k  k 

was  made  with   1„  =  3  and   1_  =  0,  which  yielded  results 

r>  A 

k  k 

identical  to  those  when   1..  =  3  and  ^1^  =  0.   Thus  the  system 

of  equations,  and  their  solution,  appear  to  be  correct  in  at 

least  this  respect. 

The  second  arid  third  parts  of  Table  IV-1  summarize  the 

results  of  varying  k_  and  k   for  the  tv/o  ships,  with  all 

other  constants  zero.   As  can  be  seen  from  the  table,  and 

from  Figure  IV- 2,  neither  k_  nor  k,-,  alone,  even  when  present 

on  both  ships,  is  sufficient  to  bring  about  system  stability. 

As  k_  is  increased,  for  one  or  both  ships,  with  all  other  k's 

zero,  system  instability  is  actually  increased.   And  as  k 

is  increased,  the  value  of  the  largest  root  is  unaffected, 

though  the  values  of  the  other  roots  do  change  slightly.   (The 

symmetry  of  the  solution  was  again  demonstrated  for  these 

constants . ) 

k     k     k       k 
Next,   3,.,   3_,,   5,  and   5D  were  all  given  values  in 

A       ti  A  D 

k       k 
their  expected  ranges,  while   L  and   1   were  varied.   The 

A  D 

results  are  summarized  in  the  fourth  part  of  Table  IV-1  and 

illustrated  in  Figure  IV-3.   In  these  circumstances, 

k        k 
increasing   1_  and   1   reduces  system  instability  until  a 

A         B 

stable  condition  is  attained. 
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TABLE    IV- 1 


Summary 

of  Stability  Roots  Obtained  for  Various 

Combinations 

of  Control  System  Constants  With  No  Time  Lags 

a 

=  9  0'; 

STSD  - 

:  0' ;   V   -  15  kts. 
o 

Ship 

A  LBP 

=  528.5' 

;   Ship  B  LB?  -  528. 5' 

All  unlisted  k.  = 
1 

0 

1. 

\ 

kl 

#  0 
Rts. 

#<0 

Rts. 

#>0 

Rts. 

Largest 
Rt.  1    0 

Stable? 

0 

0 

2 

6 

2 

.1227xl0~2 

no 

.001 

0 

1 

7 

2 

.1201xl0'"2 

no 

.01 

0 

1 

7 

2 

.9343xI0~3 

no 

.1 

0 

1 

7 

2 

.3293xl0~3 

no 

1. 

0 

1 

7 

2 

.6980xl0"3 

no 

2. 

0 

1 

7 

2 

.7186xl0*"3 

no 

3. 

0 

1 

7 

2 

•7255xl0~3 

no 

4  . 

0 

1 

7 

2 

.7289xl0~3 

no 

5. 

0 

1 

7 

2 

.7310xl0"3 

no 

0 

3. 

1 

7 

2 

.7255xl0*'3 

no 

.001 

.001 

1 

7 

2 

.1175xl0~2 

no 

.01 

.01 

1 

7 

2 

.6861xl0"3 

no 

.1 

.  1 

1 

9 

0 

~.2123xl0"2 

'yes 

1. 

1. 

1 

9 

0 

~.1877xl0~3 

yes 

2. 

2. 

1 

9 

0 

-.9270xl0~"4 

yes 

3. 

3. 

1 

9 

0 

-.  611.1  xl0""4 

yes 

4. 

4. 

1 

9 

0 

~.4534xl0~4 

yes 

5. 

J 5 • 

1 

9 

0 

-.3587xl0"4 

yes 
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TABLE    IV- 1 
(con '  t.) 


2. 


k3 
A 

k3 
B 

#  0 
Rts. 

#<0 

Rt  S  . 

#>0 
Rts. 

Largest 
Rt.  f   0 

Stable? 

-.005 

0 

2 

6 

2 

.2809X10"1 

no 

-.01 

0 

2 

6 

2 

.3720x10"* 

no 

-.03 

0 

2 

6 

2 

.5848X.10"1 

no 

-.05 

0 

2 

6 

2 

.7242X10"1 

no 

-.0  7 

0 

2 

6 

2 

.8348X10"1 

no 

0 

.01 

2 

6 

2 

.3720x10"" 1 

no 

-.005 

.005 

2 

6 

2 

.3720xl0_1 

no 

-.008 

.008 

2 

6 

2 

.4509xl0~3 

no 

-.01 

.01 

2 

6 

2 

.4940xl0~]" 

no 

-.0  3 

.03 

2 

6 

2 

.7821xl0_1 

no 

-.05 

.05 

2 

6 

2 

.9741X.10"1 

no 

-.07 

.07 

2 

6 

2 

.1122 

no 
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TABLE    IV- 1 
~ (con't) 


3. 

k5 
A 

k5 

#  0 
Rts. 

#<0 
Rts. 

#>0 

Rts. 

Largest 
Rt.  V  0 

Stable? 

.005 

0 

1 

7 

2 

•1227xl0"2 

no 

.008 

0 

1 

7 

2 

no 

.01 

0 

1 

7 

2 

no 

.03 

0 

1 

7 

2 

no 

.05 

0 

1 

7 

2 

no 

.07 

0 

1 

7 

2 

no 

0 

-.01 

1 

7 

2 

no 

.005 

-.005 

1 

7 

2 

no 

.008 

-.008 

1 

7 

2 

no 

.01 

-.01 

1 

7 

2 

no 

.03 

-.03 

1 

7 

2 

no 

.05 

-.05 

1 

7 

2 

no 

.07 

-.07 

1 

7 

2 

no 

\* 

=  -.005 

=  -k3 
B 

;   5A  = 

.01  -  - 

k5  . 

4. 

kl 
A 

kl 
B 

#  0 
Rts. 

#<0  • 
Rts. 

#>0 

Rts. 

Largest 
Rt.  ^  0 

Stable? 

1. 

1. 

0 

8 

2 

.2982X10"1 

no 

2. 

2. 

0 

8 

2 

.2163X.10"1 

no 

3. 

3. 

0 

8 

2 

.1245X.10"1 

no 

4. 

4. 

0 

8 

2 

.2264xl0~2 

no 

5. 

5. 

0 

10 

0 

-.3403xl0"2 

yes 
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TABLE    IV- 1 
(con't) " 


kl   - 
A 

k2   = 

B 

2.;  k3 

=  -.005  -  -k3 

B;  \  =  -01  " 

-k5 
B 

5. 

k2 
A 

k 
B 

#  0 
Rts. 

#<0 

Rts. 

#>0 

Rts  . 

Largest 
Rt.  f-    0 

Stable? 

.5 

.5 

0 

8 

2 

.2126X10"1 

no 

1. 

1. 

0 

8 

2 

.2090X10"1 

no 

2. 

2. 

0 

8 

2 

.  2018x10"" 1 

no 

3. 

3. 

0 

8 

2 

.1950xl0_1 

no 

4. 

4. 

0 

8 

2 

.1883X10"1 

no 

\ 

k4 
B 

-.01 

.01 

0 

8 

2 

.1996X10"1 

no 

-.04 

.04 

0 

8 

2 

.1595xl0~1 

no 

-.07 

.07 

0 

8 

2 

.1427X10"1 

no 

-.1 

.1 

0 

8 

2 

.1499X10"1 

no 

6A 

kfi 
6B 

.01 

-.01 

0 

8 

2 

.2163x10  X 

no 

.05 

-.05 

0 

8 

2 

it 

no 

.1 

-.1 

0 

8 

2 

it 

no 

.2 

-.2 

.  0 

8 

2 

ii 

no 

83 


TABLE    IV- 1 
(con1 t) 


kl   = 

A 

k2   = 

^A 

kl   = 
B 

2.;  k3A  =  -.005  -  ~k3B;  \   * 
\   -  ^V  ^A  "  ^B 

■  .01  =  -% 

/ 

6. 

k2 

|k4l 

|k6l 

#  0 
Rts. 

#<0 
Rts. 

#>o  ■ 

Rts. 

Largest 
Rt.  ?    0 

Stable? 

.5 
1. 

2. 

4. 

.01 
.04 
.07 
.1 

.01 
.05 
.1 
.2 

0 
0 
0 
0 

8 
8 
8 
8 

2 
2 
2 
2 

.1956X10"1 
.1493X10"1 

.llSlxlO-1 
.9429xl0"2 

no 
no 
no 
no 

Same,  but   1.  =  kl   -  3. 

A       r> 

7. 

.5 
1. 
2. 

.01 
.04 
.05 

.1 
.1 
.1 

0 
0 
0 

8 
8 
8 

2 
2 
2 

.1030xl0_1 
•  6420xl0*"2 
.4783xl0~2 

"  no 
no 
no 

Same,  but  kl   =  kl   =  4. 

A           D 

8. 

1. 

2. 

.04 
.04 

.  1 
.1 

0 
0 

10 
10 

0 
0 

-.1449x10  2 
-.2810x10 

yes 
yes 
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Then,  with  the  proportional  constants  for  both  ships, 
(k,  ,  k,.  and  kr),  set  at  values  within  their  ranges,  the  rate 
constants,  (k^,  k.  and  k,) ,  were  systematically  increased. 
This  step  is  summarized  in  Table  IV-1,  part  5,  and  the  effects 
for  k„  and  k.  are  illustrated  in  Figure  IV-4.   As  can  be  seen, 
increasing  k   decreased  system  instability,  and  increasing  k. 
had  the  same  effect  until  ]k.|  reached  about  .07,  after  which 
instability  increased.   Varying  k,  had  no  effect  on  the 
largest  root,  though  the  other  roots  changed  slowly  as  k_ 
increased. 

In  part  6  of  Table  IV-1,  the  effect  of  increasing  all 
three  rate  constants,  with  the  proportional  constants  fixed, 
is  shown  to  be  a  decrease  in  instability.   And  in  parts  7  and 
8,  the  same  trend  is  shown  with  the  k,  values  being  set  at 
first  3.0,  then  4.0.   This  last  combination  is  seen  to  result 
in'  a  stable  system. 

Some  investigations  with  non-zero  time  lags  are  summarized 
in  Table  IV- 2,  and  in  Figure  IV-5.   As  can  be  seen,  in  all 
cases  investigated  an  increase  in  time  lags  caused  the  system 
to  become  increasingly  unstable. 

In  Table  IV--3,  the  results  of  giving  only  one  ship  an 
automatic  control  system  are  summarized.   For  this  case  all 

i   =  0.   This  table  indicates  that  a  stable  system  can  be 
achieved  with  only  one  ship  automatically  controlled,  with 
control  constants  in  the  expected  ranges.   Physically,  this 
would  appear  to  indicate  that  one  ship  has  been  "captured" 
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TABLE    IV- 2 


Summary  of  Stability  Roots  Obtained  for  Various 
Combinations  of  Control  System  Constants  With  Time  Lags 

a   =  90*;  STSD  =  0';  V   =15  kts. 

o 

Ship  A  LBP  =  528.5";  Ship  B  LBP  =  528.5' 

All  unlisted  k.  =  0 
l 

klA  =  klB  =  2.,-  k3A  =  -.005  =  -k3B;  \   =  .01  =  -k5fi 

At, 

A 

At 
A 

At 
B 

At 

B 

#  0 
Rts. 

#<0 
Rts  . 

#>0 
Rts. 

Largest 
Rt.  ^    0 

Stable? 

1 
2 

4 

2 

4 
6 

1 
2 

4 

2 
4 
6 

0 
0 
0 

8 
6 
6 

2 

4 
4 

.2399xl0_1 
.2637x10 

•3116X10"1 

no 
no 
no 

Same,  plus  k2A  =  k20  =  2. 

1 

2 

4 

2 

4 
6 

1 
2 

4 

2 

4 
6 

0 
0 
0 

8 
6 
6 

2 

4 
4 

.2246X10"1 
. 2477X10"1 

.2951xl0_1 

no 
no 
no 

Same,  plus  k4A  -  --.04  =  -k4- 
J\                                       B 

1 
2 

4 

2 
4 
6 

1 
2 

4 

2 

4 
6 

0 
0 
0 

8 
6 
6 

2 

4 
4 

.2027xl0_1 
.2452X10"1 
.  3278x10*" 1 

no 
no 
no 
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TABLE    IV- 2 
(con't) 


Same, 

plus 

6A  " 

.1  =  - 

GB 

At 
A 

At 
A 

At 
XB 

At 
B 

#  0 

#<0 

#>0 

Largest 

Stable? 

Rts. 

Rts. 

Rts. 

Rt.  /    0 

1 

2 

1 

2 

0 

8 

2 

.2399X10"1 

no 

2 

4 

2 

4 

0 

8 

2 

.2637X10"1 

no 

4 

6 

4 

6 

0 

8 

2 

.3116X10"1 

no 

Same , 

plus 

<2   =  ] 

A 

c2   = 
B 

2.,-  k4fl  =  -. 

04  = 

-  4  ;   6   = 
B'    A 

.1—  6B 

1 

2 

1 

2 

0 

8 

2 

.1780xl0_1 

no 

2 

4 

2 

4 

0 

8 

2 

.2362xl0-1 

no 

4 

6 

4 

6 

0 

8 

2 

.3473xl0_1 

no 
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TABLE  IV- 3 


Summary  of  Stability  Roots  Obtained  for  Various 
Combinations  of  Single  Ship  Control  With  No  Time  Lags 


a  =  9  0';  STSD  -  0';  V   =15  kts . 

o 

Ship  A  LBP  =  528.5';  Ship  B  LBP  =  528. 5 » 


.005;  k5A  -  .01;  All  kiB  =  0 


'A 


.5 
1. 
2. 

A  . 


k 


A 


.01 
.04 
.07 
.1 


k 


A 


.  01 
.05 
.1 
.2 


#  0 

Rts 


#<0 
Rts 


10 


#>0 

Rts 


Largest 
Rt.  ?    0 


6793x10 


-2 


3092x10 


3914x10 


2852x10 


-5 


--3 


Stable? 


no 


no 


no 


yes 


-.005;  k4A  =  -.05;  k5A  =  .01;  k6A  =  .1; 


'A 


A 


lB  =  0 


"A 


3. 

4. 


10 
10 


4373x10 


6003x10 


-3 


-3 


yes 
yes 


9.1 


by  the  other  through  the  interaction  force  and  moment.   This 
was  an  unexpected  result  and  it  was  felt  that  it  might  be 
indicative  of  an  error  in  the  formulation  or  solution  of  the 
equations . 

To  investigate  this  possibility,  some  extra  computer 
runs,  for  special  circumstances,  were  made.   In  these  runs, 
all  the  interaction  force,  moment  and  derivative  inputs  for 
program  SOLVE  were  set  equal  to  zero.   This  is  equivalent  to 
the  ships  being  infinitely  far  apart,  where  no  interaction 
effects  can  be  felt.   In  this  circumstance,  it  should  prove 
possible  to  obtain  a  stable  system  with  reasonable  control 
constant  values  for  each  ship,  and,  this,  in  fact,  was  the 
case. 

A  further  check  was  made,  where   the  ships  were  given 

control  systems  with  different  constants then  run  again 

with  the  same  two  sets  of  control  constants,  but  switched 
from  one  ship  to  the  other both  runs  made  with  the  inter- 
action effects  all  set  at  zero.   These  two  runs  should  result 
in  identical  sets  of  stability  roots  if  the  solution  is 
correctly  done,  and,  did,  in  fact,  do  so. 

The  last  check  run  to  verify  the  accuracy  of  the  solution 
method  was  with  all  control  constants  zero  for  both  ships. 
The  interaction  effects  were  all  left  at  zero,  as  well.   This, 
in  effect,  puts  the  ships  infinitely  far  apart,  with  no 
automatic  control,  and  should  result  in  neutral  directional 
stability.   .  The  result  Of  this  run  was  four  zero  roots  and 
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six  negative  roots,  indicating  a  dynamically  stable  system, 
neutrally  stable,  probably  with  regard  to  heading 

IV.   Conclusions 

The  investigation  reported  herein  has  prompted  the  follow- 
ing conclusions: 

A.  The  methods  of  formulation  and  solution  of  the  equations 
of  motion  appear  to  be  correct. 

B.  It  appears  possible  to  achieve  a  stable  system  of 
two  ships  on  parallel  courses  in  close  proximity, 
using  automatic  control. 

C.  System  stability  is  most  affected  by  control  sensitivity 
to  heading  angle  error,    — -that  is,  k   terms;  and 
least  affected  by  sensitivity  to  longitudinal  pos- 
ition error  and  error  rate,    and   that  is,  k^ 

and  k.  terms. 
6 

D.  It  is  possible  to  achieve  a  stable  system  with  pro- 
portional control  only,  (K2        k.    k    0). 

E.  It  is  possible,  in  isolated  circumstances,  at  least, 
to  achieve  a  stable  system  with  only  one  ship 
employing  automatic  control. 

P.   Under  some  circumstances,  increasing  the  magnitude  of 
a  control  constant,  and  thereby  increasing  system 
sensitivity  to  one  of  the  error  parameters,  can  have 
an  adverse  effect  on  system  stability,  (as  illustrated 
in  Figure  IV-2). 
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V.   Recommendations 

The  fallowing  recommendations  are  made  fc^  future 
investigation  into  the  practicality  of  employing  automatic 
controls  for  ships  engaged  in  operations  alongside  and  close 
aboard,  as  in  underway  replenishment: 

A.  A  detailed  examination  of  the  theoretical  development 
used  in  this  thesic  should  be  conducted  by  someone 
other  than  the  author,  to  remove  any  reservation  as 
to  the  accuracy  of  the  method. 

B.  The  examination  of  system  stability  should  be  made 

for  two  ships  whose  full  set  of  hydrodynamic  derivatives 
are-known.   (-*-his  would  eliminate  to  need  to 
synthesize  the  characteristics  of  two  different  ships, 
as  the  author  did  with  the  OLNA/MARiNER  modification). 

C.  The  results  of  any  such  theorectical  investigation 
should  be  confirmed  by  model  tests. 

D.  The  method  should  be  extended  to  examine  the  effect  of 
including  an  outside  excitation  on  the  system,  such 

as  waves,  wind,  etc. 
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I  * P PGndix  I 
COMPUTER  UTILIZATION 


I.   Discussion 


Three  computer  programs  have  been  written  as  part  of  this 
investigation.   The  three  programs,  described  in  Appendices  II , 
III  and  IV,  are  used  in  sequence  if  the  investigation  of  the 
stability  of  a  two  ship  system,  such  as  that  studied  herein, 
starts  at  the  level  of  known  information  with  which  this 
thesis  was  begun.   If  more  information  concerning  the  inter- 
action of  the  two  ships  being  studied  is  available,  it  may  be 
possible  to  enter  the  sequence  at  a  later  stage.   (This 
sequence  of  the  programs  is  illustrated  in  the  overall  flow 
chart,  Figure  AI-1) .   Precisely  because  any  possible  future 
use  of  these  programs  may  commence  at  a  different  stage,  the 
three  programs  were  not  unified  into  a  single  large  program. 
The  use  of  the  programs  will  be  summarized  in  the  reverse 
order  from  that  in  which  they  were  used  in  this  work . 

Program  SOLVE,  the  third  program  used,  is  the  heart  of 
the  computer  solution  and  actually  solves  for  the  system 
stability  roots.   It  also  requires  the  greatest  amount  of 
known  information.   If  all  the  interaction  and  open  water 
hydrodynamic  derivatives  for  Ship  A,  plus  the  interaction 
derivatives  for  Ship  A  in  Ship  B's  position,  are  known, 
program  SOLVE  may  be  used  directly. 

Program  DETERM  can  be  used  to  get  the  force  and  moment 
derivatives  with  respect  to  propeller  speed,  n,  and  the 
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interaction  derivatives  Y   and  N   for  both  ships,  if  these 

quantities  aie  not  known.   To  do  this,  DETERM  uses  the  values 

of  Y  and  N  at  the  positions  of  interest,  plus  Yr,  Nr  and  6 

^  r  o    o      o 

open  water  for  Ship  A  as  inputs.   EHP  and  RPM  data  for  Ship  A 
are  also  required  inputs.   The  output  of  DETERM  then  forms 
part  of  the  required  input  for  SOLVE. 

Program  INTCOF  is  used  to  determine  the  interaction 
derivatives  of  Y  and  N  with  respect  to  a  and  3/  and  the 
values  of  these  quantities  at  the  positions  being  examined 
can  then  be  used  as  part  of  the  SOLVE  input.   In  the  process 
of  getting  these  interaction  derivatives  throughout  the  range 
of  interest,  program  INTCOF  also  calculates  Y  and  N  everywhere 
in  that  range,  using  the  interpolation  technique  described  in 
Chapter  III  and  illustrated  in  Figure  III-l. 

The  way  in  which  INTCOF  and  DETERM  were  written  was 
largely  dictated  by  what  information  was  available  concerning 
the  interaction  phenomena  and  were  steps  toward  getting  all 
the  derivatives  needed  for  SOLVE.   As  a  result,  INTCOF, 
especially,  would  hardly  be  likely  to  be  useful  to  an 
investigator  starting  with  a  different  level  of  known  inter- 
action information.   (If,  for  instance,  Y  and  N  were  knov/n 
from  model  tests  everywhere  in  the  range  of  interest,  the 
interpolating  portion  of  INTCOF  would  not  be  needed  and  the 
derivatives  with  respect  to  a  and  3  could  be  obtained  more 
directly).   Program  SOLVE,  however ,  should  be  directly  useful 
for  two  ship  system  stability  investigations  for  any  two  ships 
whose  hydrodynamic  derivatives  are  all  known. 
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Appendix  II 
PROGRAM  INTCQF 

I.  Discussion 

II.  Flow  Chart 

III.  Program  Listing 

IV.  Results 
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Appendix  II 
PROGRAM  TNTCOF 

I.   Discussion 

Program  INTCOF  takes  the  values  of  Y  and  N  where  known 

from  reference  A-l,  interpolates  to  get  Y  and  N  everywhere 

in  the  range  of  interest,  then  calculates  Y  ,  Yn ,  N   and  N0 

a   3   a      8 

everywhere  in  this  range,  using  subroutines  DERIV  and  DERIVE. 

Program  INTCOF  is  highly  specialized  to  the  use  of  the 
data  of  reference  A-l.   Any  investigation  using  Y  and  N 
information  from  another  source  to  get  the  above  four  inter- 
action derivatives  would  be  unable  to  use  INTCOF.   And,  of 
course,  an  investigation  using  the  [A-l]  data  could  use  the 
results  herein,  directly.   For  these  reasons,  the  input  format 
required  for  INTCOF  will  not  be  described  in  detail,  and  the 
remainder  of  this  Appendix  will  merely  contain  a  flow  chart, 
program  listing  and  results  for  INTCOF. 
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ft   Flow      Chart Program     Intcof 


INTCOF      pg  1 


5  TART 
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FROM    iwtco; 
P&  1,  #1.2 


*2.1 


INTCOF     PG.  2 


V7to     INTCOF 

P&    1«  1.3 
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JUBnoi'TINt     DERIV      (NXVy  XINiT,  XDELT,  M,Y,    SLOPl;) 


DERlv      PG.  1 
(itrrcor) 


FROM     IWTCOF 
PG.  2.    K  2.3 


XUP  -  X-l-  2.-K  XDELT 
KP=  K+2 

SLOPE  (Kh  O.O 

SLOPE(K)^(y(M,KP) 
-Y(MJ<))/ 
(XUP  -X) 


XMID   =   X+XDELT 
K^    K  +  l 
X  a    XMID 


KP  =NXY        ?     Valsf. 


TRUE 


V 

RETURN 
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DERIV6       PG.    1 

(ifvjTcor) 
SfBROUT\WE      DeRIVB      (iSlPTS,  BINIT,  BJM?M,Y,    SLOPE) 


FROM     IN/TCOF 
PG    ?.,     #  2.4 


* 


K=  1 
B  =  B  I N  IT 


k  p  =  k  +  a 

SLOPE  (K)-  o.o 
SLopc:  (i<)=  (V(kp,m: 
-Y^mIJAbup-b) 


BNMD  -   B-v-BINC 
K  =   K  +  1 

B  -    BMI  D 


K  P    -    N  P  T  S  ? 


TRUL 


V 

RETURN 


FALSE 
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IS    Y     FCJ 

O.N    SHIP    A 

ST  SO     IS    SI! 

)E    TO    SIDE 

D I  S 1  ft  MC !     8  E  "I  W  1  1 

AT     A  LEA    = 

-6CO. 

.00 

STSD= 

50.  C 

Y= 

-0. 1400E-03 

STSC= 

60.0 

Y- 

-0.1263E-0  3 

STSD= 

70. 0 

Y= 

-0.1160F-     ' 

STSO- 

8  0.0 

Y= 

-0.1058E-03 

STSD  = 

90.  0 

Y= 

-0. 9753 E-04 

ST  SC- 

100.0 

Y- 

-0.9000E-04 

AT    ALFA    = 

-550. 

,00 

ST  SD= 

50,  0 

Y  = 

-0. 2400E-03 

STSC= 

60.0 

Y- 

-0. 2112E-03 

STSD= 

70.0 

Y  = 

-0.  1897E--03 

si  so- 

80.  0 

Y= 

-0.  1681E-03 

S  T  S  C  = 

90.0 

Y  = 

-0.  1508E-C3 

STSD= 

100.  0 

Y- 

-0.1350E-03 

AT     ALFA    = 

-500. 

,cc 

STSC- 

50.0 

y= 

-0.3200E-03 

STSD= 

60.  0 

Y= 

-0.2  762E-03 

STSO- 

70.  0 

Y  = 

-Oo 2433E-03 

STSC- 

30.0 

Y= 

-0.2104E-03 

STSO- 

■0.0 

Y= 

-0. 1841E-03 

STSC- 

100.0 

Y  = 

-0.1600E-03 

AT    ALFA    - 

-45  0. 

0  0 

STSO- 

50.  0 

Y- 

-0. 3800E-03 

ST  SC- 

60.0 

Y= 

-0.3225E-03 

SI  SC- 

7 0.0 

Y  = 

-0.2793E-03 

SI  SC- 

8  0 .  0 

Y- 

-0. 2362E-03 

STSO- 

90.0 

Y- 

-0.2016E-03 

ST  SC- 

1 00.  0 

Y- 

-3. 17  00E-03 

AT     £LFA    - 

~4CC. 

,cc 

STSO= 

50.0 

Y  = 

-0.4000E-03 

STSC- 

60.  0 

Y  = 

-6.3337E-03 

STSC- 

.70.0 

Y  = 

-0.2  840E-0  3 

ST  SC- 

BC  .0 

Y  = 

-0.2342E-03 

ST  SC- 

90.C 

Y  = 

-0. 1945E-03 

SI  SC- 

100.0 

Y- 

-0. 1580E-03 

AT    ALFA    - 

-350. 

0  0 

ST  SC- 

50.C 

Y- 

-0. 3500E-03 

SI  SC- 

60.0 

Y- 

-0.2370E-03 

ST  SC- 

7 0.  0 

Y  = 

-0. 2397E-03 

SI  SC- 

80.0 

Y- 

-0. 19  2 5E -0  3 

S  7  S  C  = 

90.0 

Y- 

-0. 1547E-03 

ST  SC- 

1  00 .  0 

Y  = 

-0. 1200E-03 

SIM  '■'  S    A  G  ) 
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AT     ALFA    = 

-3i 

,  0  0 

STSD= 

50.0 

Y^ 

-0,  2  ■         -03 

STSC-= 

60.0 

Y= 

-0.2079E •  C   i 

S  T  S  D^ 

70.0 

Y  = 

-0. 1689E-03 

STSC= 

8  0.0 

Y= 

-0.1299E-03 

ST  SD= 

90.  0 

Y= 

-0.9863E-04 

STSC= 

100.0 

Y  = 

-0.7000E-04 

AT    /5LFA    = 

-2  50, 

.00 

SI  SD= 

r)  0  ,  r 

Y= 

-0.  1  3 

S  T  S  C= 

60.  0 

Y= 

-0.9384E-04 

SI  SC- 

70 .0 

Y  = 

-0.6671E-04 

SI  SC  = 

8  0.0 

Y= 

-0.3959E-04 

STSC= 

90.0 

Y  = 

-0. 1739E-04 

STSD= 

100.  0 

Y  = 

0. 2000E-05 

AT     ALFA    = 

-2  00. 

00 

STSD= 

50.  C 

Y  = 

0.4000E-04 

STSC= 

60.0 

Y= 

0. 5507E-04 

STSC  = 

7  0  ..  0 

Y= 

0.6637E-04 

STSD= 

a  o.o 

Y= 

0. 7767E-04 

STSC  = 

9  0.0 

Y= 

0.  8671E-04 

ST  SD= 

100.0 

Y  = 

0.95C0E-04 

Al     ALFA    - 

-150c 

0  0 

S  T  S  C  = 

50.0 

Y= 

0.27  0  l 

STSD= 

60.  0 

Y= 

0.25O3E-03 

STSC  = 

70.0 

Y= 

0. 2  364E-03 

STSC  = 

80.0 

Y:: 

0.2221E-03 

S  T  S  C  = 

90.0 

Y= 

0.2105E-03 

STS  D 

100  c  0 

Y= 

0.2  000E-03 

AT     ALFA    = 

-  1  0  0 . 

0  0 

S  T  S  C  = 

50.0 

Y  = 

0.5000E-0  3 

STSD= 

60.0 

V= 

0.4507E-03 

STSC= 

70.0 

Y= 

0.4137E-03 

SI  SD  = 

80,0 

Y= 

0. 3 767 E -03 

STsr- 

90.0 

Y= 

0.3471E-03 

S  7  S  C  = 

100.0 

Y- 

0.  3200E~03 

AT    ALFA    = 

-50. 

OC 

STSD  = 

50.0 

Y= 

0.7000E-03 

STSC= 

60.0 

Y= 

0. 6225E-03 

STSC  = 

7  0  . 0 

Y  = 

0.5643E-0  3 

SI  SD= 

30.  C 

y  = 

0. 5062E-03 

STSC= 

9  0.0 

Y= 

0.4596E-03 

ST  SO- 

100  c  C 

Y- 

0.4170E-03 
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AT    *LFA    = 

0. 

00 

S  T  (: 

50.  0 

y= 

0.    1 

S  T  S  C = 

60.0 

Y  = 

0. 7500E-03 

STSD  = 

70.0 

Y  = 

0. 67 50 E -03 

STSD  = 

8  0.0 

Y= 

0.  6000E-03 

STSC= 

90.0 

Y= 

0.5400E-03 

STSD  = 

100.0 

Y= 

0.4850E 

A  I     fl  L  F  a    = 

50. 

0  0 

STSD- 

5  0.  0 

Y= 

0  .9000E-03 

STSD= 

60.  0 

Y= 

0. 7  344E-03 

STSC- 

70.0 

Y  = 

0.6977E-03 

STSD= 

80.0 

Y= 

0.6110E-0  3 

STSD= 

90.  0 

Y= 

i .  5  4  L  6  E  -  0  3 

STSO  = 

100.0 

Y= 

0. 4 7 80 E- 0  8 

AT     AL FA    = 

ICO. 

00 

STSD= 

50.0 

Y^= 

0.8200E- 0  3 

S  T  S  C  = 

60.0 

Y= 

0.7159E-03 

STS0  = 

70.  0 

Y^ 

0.  6378E-03 

STSC- 

80.0 

Y  = 

0.5597E-03 

STSD  = 

90.  C 

Y  = 

0.4973E-03 

ST  SC- 

100.0 

Y= 

0.440CE-03 

AT    ALFA    = 

150. 

00 

STSD= 

50.  0 

Y  = 

<  6250E-03 

S  T  S  D= 

60  .0 

Y  = 

0. 85 79 E- 0  8 

SI  S0  = 

?  0 .  0 

Y  = 

0.5075E-03 

STSC- 

30.  C 

Y= 

0.4572E-03 

STSD= 

90.0 

Y= 

0 .4169E-03 

STSC  = 

1  CO.  0 

Y= 

0.3 8  OOF -0  3 

AT    ALFA    = 

200. 

00 

SI  SD= 

5.0.  0 

Y= 

0.4500E-0  3 

STSC= 

60.0 

Y= 

0.4089E-03 

SI  sc= 

70.0 

Y  = 

C  .  3  7;!  IF -0  3 

STSD= 

80..  0 

Y  = 

0.3473E-03 

ST  SC- 

90.0 

Y= 

0. 322  6E-03 

SI  SD  = 

100.  0 

Y  = 

0.3000E-03 

AT     ALFA    = 

2  50. 

OC 

S  T  S  D  - 

5  0.0 

v- 

0.3000E-03 

STSO  = 

60.0 

Y= 

0.2781E-03 

STSC= 

70.0 

Y= 

0. 2616E-03 

STSC  = 

8  0.0 

Y  — 

0.2452E-03 

STSD= 

90.  0 

Y= 

0. 2  32  IF -0  3 

SI  sc  = 

100.0 

Y= 

0.2200E-03 

11! 


A  f    ALFA    = 

300. 

00 

STSD  = 

50.0 

Y  = 

0.17  K  --0  3 

STSD= 

60.  C 

Y= 

0.167 3 E -  0 3 

S  7  S  C  = 

70.0 

Y= 

0. 165  2E-03 

STSD= 

30.0 

y= 

0. 1632E-03 

STSC= 

90.  C 

Y= 

0. 1615E-03 

STSD= 

100.0 

Y  = 

0.1600F-Q3 

AT     ALFA    = 

350. 

CO 

ST  si: 

50.0 

Y= 

0.6000E-04 

STSD= 

60.0 

v- 

0.6822E-04 

STSC= 

70.0 

Y= 

0.7438E-04 

S  T  S  C  = 

80.0 

Y  = 

•0.805  5  E-04 

STSD  = 

90.  0 

Y= 

0.  r55^0F-04 

s  t  s  f .:  = 

loo  .o 

Y= 

0.90 OOF- 04 

AT     UFA    = 

400. 

STSO= 

50.  0 

Y= 

-0. 2000E-04 

S  T  S  D  = 

60.0 

Y= 

-0. 3562E-05 

ST  SC--" 

70.  0 

Y= 

0. 876 7 E -05 

STSO 

80.0 

Y= 

0. 2110E-04 

S  T  S  C  = 

90  . 0 

v  = 

0.3096E-04 

STSD= 

10  3.0 

Y= 

0.4  0 00E -04 

AT    ALFA    = 

4  50. 

i"  Q 

ST  SD= 

50.0 

Y= 

-0. 5000E-04 

STSD  = 

60.  C 

Y= 

-0. 3493E-04 

STSC  = 

70.0 

Y= 

-0. 2 363E-04 

STSD  = 

80.0 

Y= 

-0. 12  3  3E-04 

S  T  S  C  - 

90.0 

Y= 

-0.3288E-05 

S  T  S  C  = 

100.0 

Y= 

0.5000E-05 

AT    ALFA    = 

5  00. 

00 

STSC= 

50.0 

Y= 

-0.  6000E-04 

STSC= 

60.0 

Y^= 

-0.4904E-04 

S  T  S  D  - 

70.0 

Y= 

-0.4082E-04 

STSC= 

ao.o 

Y= 

-0. 3260E-04 

STSD  = 

90.0 

Y= 

-0. 2603E-04 

STSC= 

100.0 

Y  = 

-0.2000E-04 

AT    ALFA    = 

5  5  0. 

00 

5TSD= 

5  0.0 

Y= 

-0. 500CF-04 

STSL  = 

60  .0 

Y  = 

-0.46  99  "■-  04 

STSD  = 

70.  0 

Y= 

-0.4473E-04 

STSD= 

80.  0 

Y= 

-0. 42 4 7 E- 04 

STSO 

90.0 

Y= 

-0.'+06  5c-0  4 

STSD= 

100.0 

Y= 

-0. 3900E-04 
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AT  A  L  F A  = 

6C0. 

oc 

STSD= 

50.  C 

Y 

STSO= 

60.0 

Y 

STSD- 

70.  0 

Y 

STSC= 

80.0 

Y 

STSD= 

90.0 

Y 

STSD= 

100.  0 

Y 

-0.  io":;r-o/; 

-0.2096E-04 
-0.  ?9i8E-0'4 
-0.3740E-04 
-0.-4397E-04 
-0.5000E-04 
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STS! 

5  0  . 0 

A  L  F  A  = 

-550  .0 

DY/DA= 

-c.  i    : 

STSD= 

50.0 

ALFA- 

-600.0 

DY/     • 

-c .  1 

STSC- 

50.0 

A  L  F  A  - 

-  4 50.0 

DY/DA= 

-0. 3000E-06 

stsf- 

50  .0 

ALFA= 

-400.0 

D  Y  / 1)  / 

0. 3000E-06 

STSD= 

50.0 

A  L  F  A  - 

-  3  5  0  ,  0 

DY/DA= 

0.140  )  E - 0  5 

STSC- 

50.0 

ALFA- 

-'300.0 

DY/  >A  = 

0.2200E-05 

STSD= 

50.0 

ALFA- 

-250  .0 

DY/OA- 

0 .  3       :    -  0  5 

STSU- 

50.0 

%  LF  A  = 

-2  00.0 

D f/DA  = 

0.4(   )0  :-05 

STSD= 

50.0 

A  L  F  A  - 

-150.0 

DY/DA= 

0.  46  1 

STSD= 

__5_C  .  0 

ALFA- 

-100.0^ 

DY/DA= 

0 .  F             -  0  5 

STSD  = 

5  0.  0 

ALf 

-50.0 

DY/DA= 

0. 3  50 )E-05 

STSD= 

5  0.0 

ALFA- 

0.0 

E>Y/DA  = 

0.  2     .'       -05 

STSC= 

50  .0 

ALFA- 

50.0 

DY/DA= 

-0.3    DOE - 

S7SD- 

50  .  0 

ALFA- 

100.0 

DY/DA= 

-0  .27 50 E-05 

STSO= 

5  0.0 

ALFA= 

150.0 

DY/O 

-0. 37  0  E-05 

STSO= 

50.0 

ALFA- 

200.0 

0Y/DA= 

-0.3250E-05 

STSO= 

50,0 

ALFA- 

2  5  0.0 

DY/OA= 

-0.23C DE-05 

STSC  = 

50,0 

ALFA- 

30  0.0 

OY/OA- 

-0. 2400 E-05 

STSD- 

50.  0 

ALFA- 

3  50.0 

DY/OA= 

-  0 .19           -  0 5 

STSO- 

50.0 

ALFA= 

4  00.0 

DY/DA= 

-0.110 

S  "f  SD= 

50..  0 

AL  FA- 

4  5  0 .  0 

DY/DA= 

-o.f  )c  ;f-06 

STSC= 

5  0  •  0 

UFA- 

500.0 

DY/OA= 

0 .  0  0  0  0  E 

S 1  S  C = 

50.0 

A  L  F  \  = 

550  .0 

OY/OA- 

0.5000E-04 

STSO= 

50,0 

A  L  F  A  = 

-550.0 

0Y/QA= 

-0.1499E-05 

STSO= 

60.0 

Al  FA- 

-  5  00  .0 

OY/0  \- 

-C . 1 112E-05 

ST  SD= 

60  .  0 

M  FA= 

-450.0 

DY/DA= 

-0.5753E-06 

SI  SD- 

60  .0 

MFA= 

-400. C 

DY/OA= 

0.  3  54  BE.- 06 

S  "T  S : 

60.  0 

ALT 

-350,0 

OY/DA= 

0.1259 

S  T  S  D 

6  r    .  0 

A!  FA-  ' 

-  3  0  .     ( 

DY/f)/ 

0. 1932E-05 

SI  5'' 

60 

A  L 

250  .0 

OY/D  i 

0  .26  ■ 

s  r  s  n = 

60  -  0 

A, '   F 

•  200  .0 

o;/;. 

0. *447    -   3  ■ 

S  I  '  C 

60,  0 

A-= 

-150 

>,  .  3956 F    05 

;rsc 

-■■'  ,  0 

\         :     A 

-  10  0- o 

DY/OA- 

),   oi  V 

■    1  SD  - 

ATI     A    : 

..;.:.  5      ■ 

'.    ]  ')>\~ 

0  ,29s  i 

.     SI 

6<  I  .0 

!  f  >■■ 

0.  0 

0  .1619! 

S 1  s 

1  ,0 

■ 

60  ,0 

)Y/ 

'•■'  11!      I>6 

"    _    "; 

loo.  : 

"  '  \A- 

-C.  2 26 51 

STSC  = 

50  .0 

ALFA- 

150.0 

DY/DA= 

-0.3070E-     5 

SI  SC- 

60 .0 

ALF  A- 

200.0 

DY/  )A  = 

-0.2793F- 

SI  SD= 

60,  0 

ALF  \- 

2  5  0.  0 

OV/DA- 

-  0  .  2  4 1 6  E 

STSO= 

6  C  .  0 

ALFA- 

300.0 

DY/DA= 

-0. 209OE-05 

5  T  SD= 

60.0 

ALFA= 

350.0 

DY/DA= 

-0.  1  703 E-0  5 

STSO= 

6  0.0 

ALFA- 

4  00.0 

DY/DA= 

-0.1032E- 

5TSD= 

60.0 

Al  FA= 

450.0 

JY/OA- 

-0.  A5'.5c-06 

STSC= 

60.0 

AL  FA= 

5  00.0 

0Y/0\  = 

-0.1205E-06 

SI  S0= 

60,0 

ALFA  = 

550.0 

OY/OA- 

0. 2 3 03 E -06 
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STSO= 

/  ;,  0 

ALFA  = 

-850.0 

DY/DA= 

-0.1273E- 

STSD= 

7  0 «  0 

ALFA- 

-500  c  0 

DY/OA= 

-0.396          06 

STSD  = 

70*0 

ALFA= 

-450.0 

DY/DA= 

-0.4                36 

STSD= 

70.0 

ALFA- 

-400.0 

DY/DA= 

0, 

STSD= 

70,0 

ALFA= 

-3  50.0 

DY/OA= 

0.1151E-05 

STSD= 

70.  0 

ALFA= 

-300.0 

DY/DA= 

0.  17  30E-0  5 

STSC- 

70.0 

A  L  F  A  = 

-2  50.0 

DY/OA= 

0 .  ?  3  5  3  E  -  0  5 

STSO= 

7  0.0 

ALFA- 

-200.0 

DY/DA= 

G. 3032E-05 

STSO= 

70.0 

ALf 

-150.0 

0  Y  / ' ) 

0.34  7       - 

STSD= 

70.0 

ALF4  = 

-100.0 

DY/DA= 

0.  J279E-05 

STSD= 

7  0 .  0 

ALFA- 

-50.0 

DY/DA= 

0  .  1  5  J    i ' 

STSD  = 

7  0.0 

ALFA= 

0.0 

OY/OA- 

0.1334E- 

STSO= 

70.0 

ALFA- 

50.0 

DY/D/ - 

-0.3719E-06 

STSD  = 

70.0 

ALFA= 

100.0 

DY/OA= 

-0. L901E-05 

SI  SD= 

70. 0 

ALFA^ 

15  0.0 

OY/OA- 

-0.259 fE-05 

STSC= 

70.0 

ALFA- 

2  0  0.0 

OY/OA- 

-0.2459E- 

ST  SC- 

70.0 

ALFA- 

2  5  0  .  0 

OY/OA- 

-0.2129E- 0 5 

SI  SD= 

70.0 

ALFA  = 

3  0  0.0 

DY/DA= 

-0.187  3C-05 

STSC  = 

70.0 

ALFA- 

3  50.0 

DY/O  \  = 

-0. L564E-05 

S  T  S  D = 

70.0 

ALFA= 

4-00.0 

0  Y  /  0  \  = 

-0.930  IE-06 

STSC- 

70.0 

ALFA- 

4  50.0 

DY/DA= 

-C. 4^5 ?E-06 

STSD= 

70.0 

A  L  F  A = 

500.0 

OV/OA- 

-0. 2110E-06 

STSD= 

70.0 

ALFA- 

5  5  0.  0 

DY/OA= 

0.1164E-06 

STSD= 

80.0 

ALF'A- 

-  5  5  0  .  0 

OY/OA- 

-0. 1047E-05 

STSO  = 

8  0  .  0 

ALFA= 

-500.0 

DY/DA= 

-0.6803E-06 

STSO= 

8  0.0 

ALFA= 

-450.0 

OY/DA- 

-0.238  4  E-0 6 

S  T  S  D  - 

80.0 

ALFA= 

-400.0 

DY/0A= 

0.4370E-  )6 

STS~D= 

80.0 

ALFA- 

-3  50.0 

OY/OA - 

0.1044E-O5 

STSC= 

8  0  , 0 

ALFA-- 

n  y  /  D  A  - 

0.  I.  5  2  9-- 0  9 

STSD  = 

80.0 

AL  FA- 

-250.0 

DY/OA= 

0.20  7  5G-0  5 

ST  SC- 

8C .0 

ALL  A  = 

-2  00.0 

OY/OA- 

0. 2616E-05 

SI  SO- 

8  0.0 

ALF  A- 

-150.0 

DY/DA= 

0.2990E-  0  5 

STSD= 

80.0 

Alf  A- 

-100. 0 

OY/OA- 

0.2341E-05 

STSC= 

80.0 

AL  FA= 

-50  .0 

OY/OA- 

0.2233E-05 

STSD= 

80.  0 

ALFA- 

0.0 

OY/OA- 

0. 1048 £-0  5 

STSC- 

8  0  . 0 

ALFA- 

5  0.0 

0Y/0\= 

-0.4 J27E-0A 

STSO= 

30.0 

ALFA= 

100  .0 

DY/DA= 

-0.153  8  -T-  -  0  5 

ST  SC- 

30.0 

ALFA- 

1  5  0  <  0 

DY/O  \= 

-0.2125E-05 

SI  SC- 

80.0 

ALFA- 

200.0 

OY/OA- 

-0. 21205-05 

SI  SC- 

30. 0 

ALFA- 

2  5  0.0 

DY/OA= 

-0.1 341E-05 

SI  sc= 

8  0  . 0 

ALFA  = 

300.0 

OY/OA= 

-0. I64  7E-05 

S  I  SC- 

80.0 

ALFA= 

3  5  0  .  0 

DY/OA= 

-0. 1421E-05 

SI  so- 

8  0.0 

A  L  F  A  = 

400.  0 

OY/OA- 

-0.923 3E-06 

STSC= 

9C  .0 

ALFA- 

4  5  0.0 

DY/DA= 

-0.5370E-      - 

SI  SD= 

80.0 

A  L  F  A  = 

500.  0 

OY/OA- 

-0.3014E-06 

STSO- 

30.0 

ALFA= 

5  3  3.0 

OY/OA- 

-C.4795E 
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STSD= 

90. C 

A  L  F  \  = 

-550.0 

u  -  - 

-0. B653E-06 

S7SD= 

90.0 

ALFA- 

-500.0 

DY/DA= 

-0.50825-06 

STSD= 

90.0 

ALF7W 

-450.0 

DY/O  \= 

-0.10  36r 

STSC- 

90.0 

ALF  A- 

-400.0 

OY/MA- 

0. 4699E-06 

STSD= 

9  0 .  0 

ALFA- 

-  3 50.0 

OY/OA- 

0  .  9  5  8  4  f  • 

STSD= 

90.0 

ALFA- 

-300.0 

DY/DA= 

C .  )             -05 

ST  S0= 

90.0 

ALF  \= 

-2  5  0.0 

DY/DA= 

0.13  5  3E-05 

S  T  S  D= 

90  .0 

ALFA- 

-200.0 

DY/OA= 

0.228 

S7SD  = 

90.0 

ALFA- 

-130.0 

DY/DA* 

0.2604E-05 

STSD= 

90.0 

ALF  ^- 

-100.0 

DY/DA= 

0.24  9!  E 

STSC= 

90  .0 

ALFA= 

-50.0 

DY/OA- 

0. L929E-05 

STSD  = 

90.0 

ALFA- 

0.0 

DY/OA= 

0.319  5E-06 

STSD= 

90.0 

ALF  \= 

5  0.0 

DY/i)A= 

-0.4274E- 

STSD- 

90.0 

ALFA= 

1 0  0.0 

DY/DA= 

-0.1247E-05 

STSO= 

90.0 

ALFA- 

15  0.0 

DY/DA- 

-0. 1747E-05 

STSD= 

90.0 

ALFA- 

200.0 

OY/DA- 

-Co  1  349E-05 

STSD= 

90.0 

ALFA- 

?50.0 

DY/OA- 

-0.161 1E- 

S7SD= 

9  0.0 

.    ALFA- 

3  00.0 

DY/D'\  = 

-0. 1466E- 

STSD= 

90.  0 

ALFA= 

3  5  0.0 

DY/DA= 

-C.1305F-05 

ST  SC- 

90.0 

ALFA- 

4  0  0  .  0 

OY/DA- 

-0.  8-377L:-C6 

SI  SO- 

90.0 

ALf  A= 

450  .0 

OY/OA- 

-0. 5699E-06 

STSD- 

90.0 

ALFA- 

5  0  0.0 

DY/OA- 

-0.3737E-06 

STSD= 

90  .0 

AL  FA- 

5  50.  0 

OY/DA- 

-0.1 795E-06 

S'l  SD= 

LOO.  0 

ALFA- 

-550  .0 

DY/DA= 

-0.7000E-06 

STSO= 

IOC  .0 

ALFA- 

-5  00.0 

DY/OA- 

-0.3500E-06 

s  r  s  d  = 

IOC  .0 

ALFA- 

-4  50.0 

0  Y  /  D  4  = 

0.2000E-07 

ST  SO- 

100.0 

A  L  F  A  - 

-4  00.0 

DY/DA= 

0.5000E-06 

STSC- 

100.0 

ALFA- 

-350.0 

[)Y/OA  = 

0. 3300E-06 

STSD= 

100.0 

ALFA= 

-300.0 

DY/DA- 

0.  L220E-05 

STSD-- 

100.0 

ALFA- 

-250.0 

DY/DA= 

0.  L650E-05 

srso= 

IOC  .0 

ALFA- 

-200.0 

DY/OA= 

0.  1930E-05 

STSO= 

100.0 

ALF  A- 

-150  .0 

9Y/0-A,- 

0.2250E-C5 

STSD= 

IOC.  0 

ALFA- 

-100.0 

DY/DA= 

0.2170E-05 

STSC= 

1.00.0 

ALFA= 

-50.0 

OY/OA- 

0.  1650E-05 

STSD= 

100.0 

A  L  F  A  = 

0.0 

DY/DA= 

0.6100E-06 

STSC= 

100.0 

ALFA- 

50.0 

OY/DA- 

-0.4500E-06 

STSD= 

100  .0 

ALFA- 

100.0 

DY/QA= 

-0.91.  )  E  -  0  5 

STSO= 

100.0 

A  L  F  A  - 

150.  0 

0  Y  /  0  A  - 

-0.1  VDOE -05 

STSC- 

100  .0 

ALFA- 

2  00.0 

DY/OA- 

-  0  .  i    :     E  -  0  5 

S7SD  = 

1  oc.o 

ALFA- 

2  50.0 

DY/DA= 

-0 . 1400E-05 

STSC- 

100.0 

ALFA- 

300.0 

DY/i)A  = 

-0.1300E-05 

STSD- 

100.0 

ALFA= 

3  50  .0 

0  Y  /  0  \  = 

-0.  1200E-05 

STSD- 

100.0 

A  L  F  A  = 

400.0 

DY/OA= 

-0.3500E-06 

STSD= 

100.0 

ALFA= 

4  50.0 

DY/DA- 

-0.6000E- 

STSD= 

100.0 

ALFA- 

500.0 

DY/DA= 

-0.4400E-06 

S7SD= 

100.0 

A  L r:  A  - 

5  50.0 

OY/OA- 

-C. 3000E-06 
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ALFA= 

-6  00. C 

STSD  = 

60.0 

DY/I 

0  .  1  I 

-05 

ALFA  = 

-600.0 

STSD= 

7  < 

0.  ]  02  7E- 

-  0  5 

A  L  F  A  = 

-6  00 .0 

ST50= 

B  0  . 0 

DY/D! 

0.9247E- 

-06 

ALFA= 

-600.0 

STSD= 

90.0 

'  DY/D8= 

0.7877  6- 

-16 

ALFA- 

-5  50.0 

STSD= 

6  0 .  0 

DY/DB= 

0.. 

-0  5 

ALFA- 

-5  50  .0 

s  rso- 

70.0 

DY/D3= 

0.21 

A  L  r  A  = 

-550.  0 

STSD= 

80.0 

DY/0 

0.  Ir?4  2f:- 

ALFA- 

-5  50  .C 

STS0- 

9  0 .  0 

DY/U3= 

0.16546- 

i 

ALFA- 

-5  00  .0 

STSO- 

6  0  ,  0 

DY/OB= 

0.38366- 

-05 

ALFA- 

-500.0 

STSD= 

7  0  > 

DY/D8= 

0. 

A  L  F  A  - 

-500  .0 

STSD= 

80.0 

! )  Y  / 1 1 

0.2  9596- 

ALFA- 

-500.0 

STSD- 

9  0  . 0 

DY/D3= 

0.2  5  216- 

-05 

ALFA- 

-4  5  0.0 

STSD= 

60.0 

DY/D8= 

0.5  )  3  4  E  - 

-05 

A  L  F  A  = 

-4  50.0 

STSD= 

70.0 

OY/0  \- 

0.43  L5E- 

■  0  > 

ALFA- 

-4  50.0 

STSO- 

8  0.0 

DY/D8= 

0.3  3  •  E- 

ALFA- 

-4  50  .0 

stso- 

9  0.0 

D  Y  /  D  D  = 

0.33086- 

-05 

ALFA- 

-400. 0 

STS9- 

6  0.0 

DY/D3= 

0.58016- 

-0  5 

ALFA- 

-4  00.0 

stso- 

7  0.0 

OY/O  3  = 

0.49  736- 

-0  5 

A  L  F  A  - 

-4  00.0 

stso- 

80  .0 

DY/OB= 

— >O.A't  75tL- 

-  D5 

ALFA- 

-4  00.0 

STSO= 

9  0 .  0 

DY/08= 

0.38126- 

-0  5 

ALFA= 

-3  50  .0 

STSO= 

60.0 

DY/D3= 

0.5  5  1-'  l- 

-05 

ALFA- 

-350  .0 

STSD= 

7  0  . 0 

DY/DB= 

0.472i  6 

ALFA= 

-3^0 .0 

STS0= 

3  0.  0 

DY/D3= 

0  .425 

-35 

ALFA= 

-350  .0 

STSD= 

90.0 

D  Y /D3  = 

0. 3623E- 

-05 

A  L  F  A  - 

—  3  ■uC:   .  0 

STSD= 

6  0  . 0 

0Y/03- 

0.455  ' 

-05 

A  L  F  A= 

-30C  .0 

STSO= 

70.0 

DY/DR= 

0 .  3   i  - 

-05 

ALFA= 

-300  .0 

S  T  S  D  = 

80.0 

DY/D3= 

0.3  514E- 

ALFA- 

-300.0 

STSD= 

90.0 

DY/06= 

0  .  2  9  0  3 ;:  - 

-J  5 

A  L  F  A  = 

-250  .0 

STSO- 

- — 60.0- 

DY/O0- 

0.31646- 

■  j  5 

ALFA= 

-2  50.C 

STSO- 

7  C  o  0 

DY/0B= 

0.27126- 

-0  5 

ALFA  = 

-250.0 

STSD  = 

3  0.0 

0Y/04- 

0.2  44] 6- 

-05 

ALFA= 

-2  50  .0 

STSO- 

90.0 

DY/0B= 

0.23790- 

-05 

ALCA- 

-2  00  .0 

STSO= 

60.0 

DY/D8= 

0.13186- 

-05 

ALFA= 

-200.  C 

STSO= 

7  0.0 

DY/OB= 

0.  1130  6- 

-05 

ALFA- 

-200.0 

STSD= 

8  0  .  0 

DY/DB= 

0.10176- 

-05 

ALFA- 

-2  00  .0 

STSO= 

90.0 

DY/0  3= 

0. 86646- 

-0  6 
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ALFA= 

-1 50,0 

ST  Si- 

60.0 

DY/D! 

-1 50.0 

Si  so= 

ro. :» 

DY/OG 

-0  .  I  ■   ■ 

ALF  &= 

-150  .0 

S 1  S  D= 

9  0 . 0 

-0.  1  295  ,"-  );» 

A  L  F  A  = 

-150.0 

STS0= 

9 '-  • 

DY/DI 

L  ]      •   -  - 

AL!ZA  = 

-100.0 

STSD= 

6  0.0 

DY/DB= 

-0.431'  "  -  ,   i 

ALPA- 

-ioi  76 

STS0- 

70.0 

-0.  3  5  >',::-0  5 

ALFA= 

-100.0 

STSO- 

80. 0 

DY/DB= 

-0.2             -05 

ALF  4= 

-100  .0 

STSO- 

90.0 

DY/DB= 

-0.  2r;  JA---)3 

ALFA  = 

-5C  .0 

STSD= 

60.0 

DY/DB= 

-0.6734E-05 

A  L  F  A  = 

-50.0 

STSO= 

70.0 

DY/DB= 

-0.5  815E-05 

ALpA= 

-5  0.0 

STSO= 

30.0 

DY/DB= 

-0. 52345-05 

ALFA  = 

-50.0 

STSO= 

90.0 

DY/OB= 

-0 .4458 E-05 

A  L  F  A  = 

0  .0 

STSD= 

60.0 

DY/DB= 

-0.8750E-05 

A  I.  F  A  = 

0  . 0 

STSD= 

70.0 

DY/DB= 

-0.7500E-05 

A  L.  F  A  = 

0  ,  0 

STSD  = 

ao.o 

DY/DB= 

-0.6750E-     > 

A  L  F  A  ^-- 

0  .0 

STSD= 

90.0 

OY/DB  = 

-0.5750E-05 

A  L  F  A  = 

50.0 

STSO- 

6  0.0 

DY/DB= 

-0  .10  120-0'. 

ALFA- 

r>0.0 

STSO= 

70.0 

DY/DB= 

-0.8671E-05 

ft  L  F  A = 

50  .0 

STSO- 

3  j  .  0 

DY/D8= 

-0.7304E-05 

ALFA= 

50.0 

STSD- 

90.0 

DY/DB= 

-0 .6643E-05 

ALFA= 

100  .0 
1  0  j  .  0 

STSO= 
STSO= 

60.0 

DY/DB= 

-0. 9110E-05 

ftl_FA  = 

70.0 

DY/03= 

-0.7308E-05 

ALFA- 

100.0 

STSD  = 

8  0 .  0 

DY/DB= 

-0.7O275-~r> 

ALFA  = 

100  .0 

STSD= 

9  0.0 

DY/0J,- 

-0.5936E-05 

ALFA  = 

1  50  .  0 

STSO= 

60.0  — 

->DY/OB  = 

-0 . 5373E-      i 

ALFA= 

150  ,0 

s  rso= 

70.  C 

DY/DB 

-0.5034E-05 

ALFA= 

,    150  .0 

STSO= 

8  0.0 

DY/O  i= 

-  0 .  4  5 3 1  E-^0.5 

ALFA  = 

1  )  C  .  0 

STSD= 

9  0.0 

DY/0B= 

•    )  .3860E-05 

ALFA  = 

200.0 

STSD=       - 

60.0— 

-*DY/D  J  = 

-0.3596E-C 5 

ALFA- 

200  .0 

STSD= 

7  0 .  0 

DY/D3= 

-0. 3082E-05 

A  L  F  A  = 

2  00.0 

STSO= 

8  0 . 0 

DY/OB= 

-0.2 7 74 E-05 

ALFA= 

200.0 

STSD= 

9  0.0 

DY/D3= 

-0.2  363E-      > 

ALFA= 

2  50  .0 

STSO= 

60.3 

DY/DB^- 

-0. 1918E-0  5 

ALFA  = 

2  50,  0 

STSO- 

7  0.0 

DY/DB= 

-0. 1644E-C    i 

A  L  F  A  = 

2  50  .0 

STSD= 

30.0 

DY/DB= 

-0, 1479E-05 

A  L  P  A  = 

2  50  .0 

STSD= 

90.0 

DY/DB= 

-0.  126CE-05 

ALFA= 

3  00.0 

ST  SO-  — 

___60.5- 

DY/O 3= 

-0.2  3  ;7S-06 

ALFA= 

300.0 

S  1  S  D  - 

70. r) 

DY/OB= 

-0.2055E-06 

ALFA-- 

) .  0 

STSO  = 

8  0.0 

DY/O 3= 

-0  .1349E-  05 

A  L  F  A  - 

300.0 

STSO= 

9  0 .  0 

DY/DB= 

-0.15753-0') 
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A  L  F  ft  = 

:  .0 

STS0  = 

I ..  0  ' 

DY/D  ■  ■ 

0.71 

. 

ALFA= 

3  30,0 

S  T  S  D= 

70.0 

DY/D3= 

.  0.61  < ',  E  - 

ft  L  F  A  = 

350.0 

STSD  = 

80.0 

DY/DB= 

.5  »  '.  3  E  ■ 

-06 

ALFA= 

3  50.0 

STSD  = 

90.0 

DY/03= 

0.4726E- 

-0  6 

ALFA  = 

400.0 

STS0  = 

6 ;  •  .  0— 
7  0 .  0 

— £Y/DB= 
DY/DB= 

0.14: 
0.123 

-  0  5 

ALFA  = 

400.0 

ALFA= 

400  .0 

STSO- 

80.0 

DY/()B  = 

^r  0.1  ]  1  0  E  - 

-05 

ALFA-- 

400.0 

STSD= 

.0 

DY/DB= 

o  .9  a  j  '.:■■ 

)6 

A  L  F  A  = 

450.0 

STSD= 

60.  0 

DY/D 3= 

0.131 

i 

ALFA  = 

450.0 

STSD= 

70.0 

DY/03= 

0.1.1  ' 

-0  5 

ALFA  = 

4  50.0 

STSD= 

8  0.0 

DY/D 3= 

0  . 1 0 1 7  E  ■ 

i 

ALFA= 

ALFA-- 
A  L  f:  ft  = 

4  50.0 

5  00  .0— 
500. 0 

STSD= 
STSD= 

9  0.0 

.0 

7  0  .  0 

DY/D 3= 

0.8664E- 

— — —  '  \       DC  }  O  C  . 

-06 

-00 

UY  /  !J  '  -- 

DY/D 

0  . 8  2  1 9  E  - 

ALFA= 

5  00  .0 

srso= 

80.0 

DY/l)B  = 

0.7397E- 

-06 

A  L  F  ft- 

5  00.0 

STSD  = 

9  0.0 

DY/DB= 

0.6301  E- 

-06 

ALFA- 

5  50  .0 

STSD= 

60.0 

DY/DB= 

0.2637E- 

-06 

ALFA= 

550.0 

STSD= 

70.0 

DY/DB= 

0.2260E- 

-06 

A  L  F  A  = 

5  50  .0 

STSD= 

8  0.0 

DY/D 3= 

0.2034E- 

ALFA= 

5  50.  0 

STSO= 

90.0 

DY/DB= 

0  .  3  7  3  3  E  • 

-0  6 

ALFA= 

6  00  .  0 

STSD= 

60.0 

DY/DB= 

-0. 9589E- 

- ')  ', 

ALFA= 

600.0 

S1S0- 

7  0  . 0 

DY/DB= 

-0.8  219E- 

ALFA= 

600.0 

STSO  = 

80.0 

DY/DB= 

-0.7397E- 

-C6 

ALFA  = 

600  .0 

STSD= 

90  .0 

DY/DB= 

-  0 .  6  3  0 1 E  - 
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Y    V/;H  iBL  E 

IS    N    N 

IOMENT 

ON    SHIP    A 

A  i     A  L  F  A    = 

-600. 

00 

ST  S 'lI- 

50*  0 

Y  = 

0.  6000E-94 

ST  S  D  - 

60  cO 

Y- 

0.60  00F--0/. 

ST  SD  = 

70.  0 

Y- 

0.6000E-04 

S  T  S  D  = 

30.0 

Y- 

0.6000E-04 

STSC  = 

90.  C 

Y  = 

0. 6 000 E -04 

ST  SC- 

100 c  0 

Y  = 

0. 600CE-04 

AT     ALFA    = 

-550, 

00 

ST  SO- 

50.0 

Y  = 

0.  11  OOF -0  3 

STSC- 

60.0 

Y- 

0. 1036E-03 

STSO- 

70,  0 

Y- 

0.9786E-04 

STSD- 

30.0 

Y- 

0.9280E-04 

ST  SD= 

90.0 

Y- 

0. 3890E-04 

ST  SC- 

100.0 

•   Y= 

0. 8500E-C4 

AT     ALFA    = 

-500. 

0  0 

ST  SO- 

50.  0 

Y- 

0.1390E-03 

STSC- 

6  0.0 

Y- 

0. 12  78E-03 

STSC- 

70.0. 

Y- 

Oc  1. 176F-03 

STSD= 

8  0  .  C 

Y- 

0.  1087E-03 

STSC- 

90.0 

Y- 

0.  1019E-03 

ST  SO- 

100.0 

Y  = 

0.c'500E-04 

AT     ALFA    - 

-450. 

0  0 

STSO  = 

5  0.0 

Y- 

0.  1530E-03 

STSC- 

60.0 

Y  = 

0.  13/--E-03 

STSC- 

70.0 

Y- 

0.1 L85E-03 

STSD= 

8  0.  0 

Y.= 

0. 1042E-03 

STSC  = 

90  .0 

Y  = 

0.9308E-04 

ST  SC- 

100.0 

Y- 

0 . 3200 E -04 

AT     ALFA    = 

•    -£00. 

00 

STSC- 

50.0 

Y  = 

0. L400E-03 

STSD= 

60.  0 

Y  = 

0.1197E-03 

STSC- 

70.0 

Y  = 

0.  1012E-03 

STSC- 

8  0.0 

Y  = 

0. B497E-04 

STSD= 

90.0 

Y- 

0.7249E-04 

ST  SC- 

1 00  e  0 

Y  = 

0.6000E-04 

AT     ALFA    = 

-350. 

00 

ST^C- 

50.  0 

Y= 

0. 70C0E-04 

ST  SC- 

60.0 

Y  = 

0.5728E-04 

SI  SD- 

7  0.0 

Y- 

0.4572E-04 

STSC- 

80.  0 

Y- 

Oc  3  56  IE- 04 

STSC- 

90.  0 

Y- 

0.2780E-04 

STSC- 

10Q.C 

Y= 

0.200CE-04 
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a  i    / :  l  f  \  - 

-300. 

CO 

STSD  = 

50.0 

Y^ 

-0.50  ■ 

STSO= 

60.  0 

Y^ 

-0.4618E-04 

STSD= 

70  .0 

Y  = 

-0.4272E-04 

STSD= 

8  0.0 

Y  = 

-0.3968E-04 

STSD= 

90.0 

Y= 

-0. 3734E-04 

STSD= 

100.  0 

Y= 

-0.3500E-04 

AT    ALFA    = 

-250. 
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-C .3591 F -06 

STSD= 

9  0.C 

Al.FA- 

4  5  0 .  0 

DY/DA= 

-0.365  6  E-06 

STSO= 

9u.O 

ALFA- 

500  .0 

DY/OA= 

-C. 2999E-06 

s  rso= 

90.0 

ALFA- 

r>  5  0  .  0 

DY/DA= 

-C.  1 296 E-06 

STSC  = 

100.0 

ALFA= 

-55  0.0 

D  Y  /  D  A  - 

0.   ':■>.  OE-06 

Si  SD  = 

100.0 

ALFA- 

-5  00.0 

DY/DA= 

-C. 3000 E-0 7 

STSD  = 

1  00  . 0 

Al  FA= 

-450.0 

DY/DA= 

-  C .  3  5 

STSD= 

100.0 

AL"FA= 

-400.0 

DY/DA= 

-  C .  6  2 

STSD= 

100.0 

A  L  F  A  - 

-350.0 

DY/DA= 

-C.9500E-06 

ST  SC= 

100.0 

ALFA- 

-300.0 

DY/DA= 

-0. 1110E-05 

s  rso= 

100.0 

A.LFA  = 

-2  50.0 

DY/DA= 

-0.1370^-05 

STSD= 

100.0 

4LFA= 

-2  00.0 

OY/OA- 

-C .9400 E-06 

STSD= 

100.0 

ALFA-- 

-150.0 

DY/OA- 

-0.6300E-06 

ST  SO- 

100.  o 

ALFA- 

-100.0 

DY/OA- 

-C  .2300 E-06 

STSD= 

100.0 

ALFA- 

-50.0 

DY/OA- 

C.9000E-0 r 

STSC  = 

]  00.0 

ALFA- 

0  .0 

.  DY/DA= 

0.4300E-06 

S1SD- 

100.0 

ALFA- 

5  0 .  0 

OY/OAr- 

C .1010E- 35 

STSC= 

1  0  3  . 0 

A  L  F  A  = 

100.0 

DY/DA= 

C. 1 J00E-05 

STSO= 

100.0 

ALFA- 

1  5  0  .  C 

DY/DA= 

C.  1  20  JE-05 

STSO= 

100.0 

ALFA- 

20C.  0 

DY/DA= 

0. 1000E-05 

STSO= 

100.0 

ALF-A- 

2  50  .0 

DY/DA= 

0.6000E-06 

STSD= 

100.0 

ALFA- 

300.0 

DY/D 

0. 2200E-06 

STSD= 

100.0 

ALFA- 

3  50.0 

DY/QA= 

-C. 1500E-06 

STSO- 

1  00.0 

ALFA= 

400.0 

DY/DA= 

-0. 3700E-06 

STSD= 

100.0 

ALF  \  = 

45  0.0 

DY/9A- 

-0.3503 

STSC= 

100.0 

ALFA- 

5  00.0 

DY/DA= 

- 0 .  2 5( 

STSD= 

1  00  .  0 

\  I  F  A  = 

5  50  .0 

DY/DA= 

-C . 15001 
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At  FA- 
ALFA- 
ALT  \  = 
A  L  F  A  = 
ALFA= 
A  L  F  A  = 
ALFA= 
ALFA  = 
AL  FA  = 
ALFA^= 
ALFA= 
A  L  F  A  - 
ALFA= 
ALFA= 
ALFA= 
ALFA= 
ALFA  = 
ALFA  = 
ALFA- 
ALFA  = 
ALFA- 
ALFu- 
ALFA  = 
ALFA= 
ALFA= 
ALFA= 
ALFA- 
A  L  F  A  = 
A  L  F  A  - 
ALFA- 
ALFA= 
ALFA- 
ALFA- 
A  L  F  A  = 
ALFA- 
AL^A- 


6  00 

.0 

ST  SO 

6  00 

.0 

STSO 

6  00 

.  0 

STSD 

600 

.0 

ST  SO 

550 

.0 

STSD 

5  5C 

.0 

STSO 

550 

.0 

STSO 

55  0 

.0 

STSO 

5  00 

.0 

STSO 

500 

.0 

S  T  S  0 

5  00 

,0 

STSO 

500 

.  0 

STSO 

450 

.0 

S  T  S  0 

4  50 

.0 

s  rso 

4  50 

STSO 

'i50 

.0 

STSO 

4  00 

.0 

S  T  S  0 

400 

.0 

STSO 

4  0  0 

.0      . 

STSO 

4  00 

.0 

STSO; 

3  50 

.0 

STSD 

350 

.0 

STSD- 

350 

.0 

STSO 

3  50 

.0 

STSO^ 

300 

.0- 

STSO: 

3  00, 

.0 

STSO 

300 

.0 

STSO 

3  00 

.0 

STSD: 

2  50 

.0 

STSO 

250 

.0 

STSO' 

2  50 

.0 

STSO: 

250 

.  0 

STSO 

2'00 

.0   - 

STSO 

2  00, 

.0 

STSO 

200 

.0 

STSO: 

200 

.0 

STSD: 

60.0 

OY/D  '- 

0.00  )0E 

: 

7  0.0 

D  Y/D '•<>- 

O.OOOOE 

00 

3  o .  o 

DY/Dl 

'  0  . 0  0  0  0  E 

00 

90.0 

DY/DB= 

0  .00  ' 

60.0 

DY/DB= 

-0.60696- 

-06 

70.0 

DY/OB= 

-0.5419E- 

80.0 

DY/DB= 

-0.4', 

-06 

90.0 

DY/OB= 

-0.39  J>:; - 

60.0*" 

n  w  * n  ri   ■ 

-0.1  0< 

-05 

t ;  Y  / 1 J  b  a~ 

70.0 

DY/DB= 

-0.05 

-  '.   i 

8  0.0 

DY/D3= 

-C.7334E- 

-  36 

90.0 

DY/DB= 

-0.6  3676- 

-06 

60.0 

DY/D! 

-0.17  >4E- 

70.0 

DY/DB= 

-0.153 

-05 

8  0.0 

DY/DB= 

-0.12 

-0  5 

9  0.0 

DY/D3= 

-0.  11 

-0  5 

60.  0 

DY/DB= 

'-  0 .  j  9426- 
-0.1734E- 

-05 

70.0 

DY/DB= 

-05 

0  0.0 

DY/OB= 

-0.14 

-0  5 

90.  C 

DY/DB= 

-0.  1  2  'i  v  ~  - 

-0  5 

60.0 

DY/DB= 

-0  .  1 2  1 4  E  - 

-0.5 

70.0 

DY/03- 

-0. 1034E- 

-  0  3 

80.0 

DY/Di 

-0 . B 1 ' 

-OS 

90.0 

DY/03- 

-0  .  7  3  0  3  E  ■ 

60.0 

DY/DB— 

0.3642E- 

-06 

70.0 

DY/DB= 

0.3  251E- 

-  3  •> 

8  0.0 

OY/DB= 

0.2688E- 

90.0 

DY/DB= 

0.  2  J4  1E: 

-06 

60.0 

DY/DB=  ' 

0  . 2  0  3  9  E  • 

-05 

70.0 

DY/03- 

0 .  1 3  2 1 E  - 

-05 

8  0  . 0 

DY/DB= 

0.1505E- 

-  0  5 

9  0 .  0 

DY/DB=. 

0  .  1 3  1 1  E  - 

-  D  5 

60.0  — 

DY/DB=_ 

0.333     :- 

- :  5 

70.0 

OY/Ob= 

0.34256- 

-0  5 

30.  0 

DY/D3= 

0.2331E- 

-  3  5 

90.0 

DY/03- 

0.2  46< 

-05 
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ALFA= 

-150  .0 

STSO- 

A  LF  A  = 

-150.0 

STSO= 

ALFA= 

-150.0 

STSD= 

ALFA  = 

-150  .0 

STSO- 

ALFA- 

-100.0- 

STSO- 

A  L  F  A  = 

-100. 0 

STSO- 

ALFA= 

-100.0 

STSO- 

ALFA= 

-100.0 

srso- 

ALF/ 

-50  .0 

STSD= 

A  L  F  A  - 

-50.0 

STSD= 

ALFA= 

-50  .0 

STSO- 

ft  L  F  A  = 

-50.0 

STSD= 

ALFA= 

CO 

STSD= 

A  L  F  A = 

0  .0 

S TS  r^= 

A  L  F  ft  = 

0.0 

' S  T  S  0- 

ALFA- 

0  .0 

STSO  = 

A  L  F  A  = 

50.0 

S  T  S  0  - 

ALFA= 

50.0 

STSO= 

ALFA= 

50  .0 

STSO- 

ALFA= 

5C  .0 

STSO- 

ALFA= 

100.0 

STSO- 

ALFA- 

100.0 

STSO= 

ALFA- 

100.0 

STSO- 

ALFA- 

100.0 

STSO= 

ALFA- 

150.0 

STSO- 

ALFA= 

150.0 

STSD  = 

ALFA= 

150  .0 

STSO- 

ALFA- 

1  50.0 

siso-- 

ALFA= 

2  00  . 0 

£JSD  = 

ALFA= 

2  00.0 

STSD= 

ALFA- 

2  00.0 

STSO- 

ALFA= 

200  .0 

STSD= 

ftLFA  = 

250.0 

STSO-- 

ALFA  = 

2  50  .0 

STSD= 

ALFA= 

2  50  .0 

S  T  S  0- 

A  L  F  A  - 

2  50.0 

ST5D= 

ALFA= 

300.0 

STSD  = 

ALFA  = 

3  00  .0 

STSD= 

ALFA  = 

3  00.0 

STSO- 

ALFA= 

300.0 

STSO* 

60.0 

DY/DB= 

0.5 

C  5 

70.0 

DY/DB= 

0.4377 E - 

8  0 . 0 

DY/DB= 

0  .  4  0  3  2  E  - 

90.0 

DY/OI 

0.3512E- 

■0  5 

6  0.0 

0Y/':3- 

0  . 5  6 

•05 

70.  0 

DY/DB= 

■05 

80.  C 

DY/OB= 

0  . *  1 57  E- 

■05 

90.0 

DY/DO 

0.362J    E" 

-0  5 

60.0 

DY/0&= 

,    '      , 

■0  5 

70.0 

DY/DB= 

0  .  '•. r:  1 2  E  - 

■05 

30.0 

0Y/0  3- 

0.3978E- 

90.0 

DY/OB= 

0.3465E- 

6  0.0  - 

OY/0  0- 

0.4200E- 

■  D  5 

7  0  .  0 

DY/DB= 

0.3750E- 

•05 

8  0.0 

DY/DB= 

0.3  100 E- 

-05 

90.0 

DY/OB= 

.  0.2700E- 

-J5 

60.0 

DY/D8=— 

0.2671E- 

-  )5 

70.0 

DY/D8= 

0  .  2  3  34  E- 

-05 

8  0 .  0 

DY/DB= 

0 .  1 9  7 1 E  - 

■  0  5 

90.0 

OY/DO- 

0  .17  17  E - 

-  C  5 

6  0 .  0 

DY/OB=  ^ 

-0.1214E- 

-05 

70.0 

DY/DB= 

0.  ]  334E- 

-05 

8  0.0 

DY/DB= 

0.8960E- 

-06 

9  0.0 

DY/DB= 

0.78  33E- 

-0  6 

60.0 

DY/DB^=— 

0.1214E- 

-  1  6 

70.0 

DY/DB= 

0  .  1 0  3  4  E  ■ 

-  3  6 

80.0 

DY/DB= 

0.89 

7 

9  0.0 

DY/D 

.753?" 

-o  r 

60.0 

DY/DB=- 

-0. 5341 t 

-    v   'J 

70.0 

DY/DB= 

-0.47 

-0  6 

3  0.0 

DY/DB= 

-0.39 ^E 

-06 

90.0 

DY/DB= 

-0. 3 A J^E 

-  U  'J 

60.0 

DY/&S- 

-0.1  214E 

-05 

70.0 

DY/DB= 

-0.1 334 E 

-05 

3  0.0 

DY/D 8= 

-0. 6960E 

-    '    . 

90.0 

DY/DB= 

-0.7  JO 3E 

-0  6 

60.0 

DY/D8  =  _^ 

-0.1457E 

-05 

7  0.0 

dy'/db= 

--".13:;  e 

-0  5 

3  0.0 

DY/DB= 

-O.10  75E 

-  0  5 

90.0 

DY/DB= 

-0.^364^ 

-  0  6 
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A  L  F  A  = 

3  5C 

.  I 

ST$0 

Al  r 

350 

.0 

STSf) 

ALFA- 

3  50 

.0 

STSO 

A  L  F  A  - 

3  50 

.0 

STSO 

ALFA- 

400 

.  0 

STSD 

ALFA- 

4  00 

.0 

STSO 

ALFA-: 

4  00 

.0 

STSO 

A  L  F  A  = 

400 

.0 

STSO 

A  |.  F  A  = 

45C 

.0 

STSD 

ALFA- 

4  3C 

.0 

S  T  S  D 

ALFA- 

4  5  0 

.0 

STSO 

ALFA- 

4  5  0 

.0 

STSO 

ALFA^ 

500 

.0 

STSO 

ALFA- 

500 

.  0 

STSO 

Al.F-A- 

500 

.0 

STSO 

ALUA  = 

50C 

.0 

S  T  S  0 

ALFA- 

5r>:: 

.  0 

STSO 

ALFA- 

5  50 

.0 

STSO 

A  L  F  A  - 

5  50 

.  0 

STSO 

ALFA- 

5  50 

.0 

STSO 

A  L  F  A  = 

6  00 

.0 

STSO 

ALFA- 

600 

.0 

STSO 

ALFA- 

600 

.0 

STSD 

A  L  F  A  - 

6  00 

.0 

STSO 

•   .0 

OY, 

70.0 

DY/DG 

8  0  . 0 

DY/DB 

9  0  .  1 1 

DY/08 

6  0.0 

[V(  /■)  \ 

7  0 . 0 

DY/DB 

8  0 . 0 

DY/DB 

9  0  . 0 

DY/DB 

60.0 

DY/DB 

70,0 

D  Y  /  D  \ 

80.0 

DY/DB 

90.0 

DY/DB 

60.0 

nv/00 

70.0 

DY/D6 

80.0 

OY/  >B 

00.0 

DY/DB 

60.0 

DY/DB 

7  0.0 

DY/DB 

8  0.0 

DY/DB 

90.0 

DY/DB 

60.0 

DY/D  ! 

7  0.0 

DY/DB 

8  0  .  0 

DY/DB 

90.0 

DY/DB 

•0 
-0 
■0 

•o 

o 


0  .  1  2 1 4  E  ■ 
0  .  1 C 
0.8  l 
0.78 
0.13  J5F- 
0.11 
98  ■ 

I 

i  !  1 6 :  ■ 

.      .     ,E 
0 . 1 0  9  2  E  - 
.9  754E- 
-0.  8064E- 
.-0.  70236- 
-0.6  069E 
-0.54) 
-0.4430E' 
-0.39026- 
0.1457E- 
0.13016- 
0.10 7  - 
0.9  $64 ' 


-C 


■0  5 
•  0  5 

3  6 
•05 

3  5 
05 
•06 

05 

36 
0  6 

06 
■06 
06 

06 
■06 

06 
■07 
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Appendix  III 
PROGRAM  DETERM 

I  j Pis  cuss  ion 

Program  DETERM  determines  the  hydrodynarnic  derivatives 

with  respect  to  propeller  speed,  n,  which  were  not  available 

in  the  literature  surveyed.   It  also  determines  the  variation 

of  interaction  effects  with  velocity;  u.   Thus  the  outputs 

are  X  ,  Y   and  N   for  Shio  A;  Y   and  N   for  Ship  A;  and  Y 
n    n      n  u       u        "L  u 

and  N   for  Ship  A  in  Ship  B's  position.   These  latter  values, 
u        ^         ^ 

when  used  as  inputs  in  program  SOLVE  are  modified  to  their 

appropriate  values  for  Ship  B  in  Ship  B's  position.  (X  ,  Y 

and  N   are  position-independent.) 
n  L 

The  values  of  X  ,  Y   and  N   are  calculated  by  using  the 

n  .  n      n  J 

known  values  of  X,  Y  and  N  exerted  by  the  propeller  at 
equilibrium  speed  in  open  water  and  assuming  that,  propeller 
forces  and  moments  vary  as  the  square  of  propeller  speed. 
The  derivatives  with  respect  to  u  also  are  calculated  by 
taking  the  known  values  of  Y  and  N  at  equilibrium  and  assuming 
that  they  vary  as  the  square  of  the  velocity,  u.   (The 
procedure  followed  in  calculating  these  quantities  is 
described  in  detail  in  Chapter  III,  PORCEDURES.) 
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1 1 .   Input/Pi  ]  tput  Format 

Program  DETERM's  input  is  ordered  as  follows: 


ITEM 


PROGRAM 

CARD- 

COLUMNS 

FORMAT 

SYMBOL 

ALBP 

1 

1-  10 

F10.2 

ABEAM 

1 

11-20 

F10.2 

V(l) 

2 

1-10 

F10.2 

RPM ( 1 ) 

2 

11-20 

F10.2 

EHP  ( 1 ) 

2 

21-30 

F10.2 

V(2) 

3 

1-10 

F10.2 

RPM ( 2 ) 

3 

11-20 

F10.2 

EHP  (2) 

3 

21-30 

F10.2 

V(3) 

4 

1-10 

F10.  2 

RPM ( 3 ) 

4 

11-20 

F10.  2 

EHP ( 3 ) 

4 

21-30 

F10.2 

NPOSIT 

5 

1-4 

14 

ALFAO ( I ) 

6 

1-10 

F10.2 

LBP  Ship  A 

Beam  Ship  A 

Speed  at  1  kt.  below 
equil .  speed 

RPM  at  this  speed 

EHP  at  this  speed 

Equil.  speed  (lets.) 

RPM  at  equil. 

EHP  at  equil. 

Speed  at  1  kt.  above 
equil.  speed 

RPM  at  this  speed 

EHP  at  this  speed 

#  positions  examined 

Long,  separation  at 
each  position 

Lat.  separation  at 

each  position        STSD(I)       6       11-20         F10.2 

(one  card  for  each  position) 

Y  force  at  each 

position,  Ship  A     AY (M)         7        1-14         E14.4 

N  moment  at  each 

position,  Ship  A     AN (M)         7       15-28         E14.4 

(one    card  for  each  position) 
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ITEM  PROGRAM      CARD     COLUM         FORM 

SYMBOL 


Y  force  on  Ship  A 

in  Ship  B  position  BY (N)  8  1-14         E14.4 

N  moment  on  Ship  A 

in  Ship  B  position  BN(N)  8  15-28         E14.4 

(one  card  for  each  position) 

Y.  Shi-)  A  AYDEL  9  1-14         El 4  .  4 
o 

N.  Ship  A             .  ANDEL  9  15-28         E14 . 4 

6  at  equil,  Ship  A  ADELO  9  29-3-8         F10 .  2 


The  output  format  for  DETERM  is  illustrated  by  the 
sample  results  section  of  this  Appendix  and  should  be  self- 
explanatory.   In  the  output,  X   is  represented  by  DX/DN , 
Y   by  DY/DN,  etc. 
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1        FLOW    CHART 
START 


AZ 


ALBP;   ABEflMj 
(\'(K),  RPMlK), 
EHPIK),    K-  i,3 


NPOSIT, 
(Ai-FAO<I),  ST50(I), 
1=1   N/POSIT) 
(AY  (M),  AWfMIj 

M=  1;  NPOSlT) 
(BY(W),    BN(N), 

N  =  l    WPOSIT 


A^/'O,   AMD, 
ft  DEL  O 


rps(m)-rpm{m)/sq. 

r 


Pi 


ROGRAM 


Dete 


DETERM       PG.    1 


R  M 


X(M)=  EHP(M)¥  zSbo.  / 
(  V  ( M )  *  l.G  8  7  8  ) 

£ 


AXN-  {(xM3)-  mil/ 
(RP5(3)-  RPS(l)))/ 
(0.5  *RHOX(At_0P*X3.) 
#  V(i)^   1.687  8) 


DAVD-   AYD  xc.5s<RH0 
X(ALQP**2.)*(V(t) 
X  1.6878)#X2. 
XADEL  0/57.3 


DAVD'JP-    PAVD  * 
(RPS(3)/RP5(2))m. 


5LAYD=  (OA^OUP- 
DAYD)/(RPS(3>-RPJft)l 

AVK/=  -  SL-AND/(0.5* 
RHO*(AL6P**3.)* 

V(?-)*  1.GG78) 


DAND-    AND  ¥0.5  * 
RHO*  (ALBP«»3J   * 
((VU)  Vj.GC76)XV2.) 
tf  Aaa.0/5-7.3 
D/UJOUP--  OAUO  X 
(RPS(3)/ftPS  (3j)*«2. 
SLAND  -  (OAMOUP-PAWD 
rfRP5(3)-RPS(l» 
AMKJ  =  -  SUMO  /(O.J 
X-RMOX(AL8P**4.) 
5H'(t)i''   l.G  07  8) 


ALBP,    AOEAM 
ADEI-O,  V(2); 
RPM{2) 


AXW,   A"KMJ  A  MM] 


FROIA      DETERM 
PC- .  2,  ^  2.2    . 


±t  1.1 


DAY=  AV(L)  *(0.5* 
RVIO  X  ({V(2)*1.6878J 
*¥Z.)XAL8P**1.) 
DAVUP=   DAy*(V(3)/ 

SI.AV=  (DAVUP-DAY)  / 
((V(3)-V(2))*1.G878) 
MU*   S"LAV/(0.5*Rl->o 
*\fr)*  1.6878* 
(ALBP**Z.U 


V 


TO     DETERM 
PG.   ?.,    **  2.1 
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DETEUK     PC.   2 


FROtA 
PG.  1. 


DETERK 
.1 


*i.: 


2 


#•2.1- 


DAI^=  AN  (L^(0.5*f 
RMOx((VfI)*f  1.GB78) 
**Z.)*  AL8P#*3.) 
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Y.     SAMPLE-    RESULTS 

SHIP    A 
LBP-  528*50  BEAM  =  76,00 

OPEN    LATER    RUDDER    ANGLE=  1.20 

FQL'IL    SPEED    IN    KTS=  15.00  EOUIL    RPM  = 


65°00 


DX/DN= 
DY/DN- 
DN/DN= 


0.4623E-0A 

-0.5236E-05 

0.2615E-05 


AT  SIDE  TO  SIDE  DIST= 
ALFA  ZERO-        OoOO 


60e00 
DY/DIJ-    0. 1525E-02 
DN/DU-   -0.6710E-03 


AT  SIDE  TO  STOP  DIST= 
At.  FA  7ER0=        0.00 


90.00 
DY/DU= 

dn/du= 


Oo 105  8E-02 
•0.4636E-03 


AT    SIDE    TO    SIDE    DIST= 
ALFA    ZER0=  300.00 


70.00 

dy/du= 

DN/DU= 


0.3359E-03 

Oe 2255E-C3 


AT    SIDE    TO    SIDE    DIST= 
ALFA    ZERO=  ~400*00 


BOoOO 
DY/DU=       -0.4762E-03 

DN/DU=         0.1728E-03 
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FOR  SHIP  A  IN  SHIP  R  POSITION 


AT  SIDE  TO  SIDE  DlST= 
ALFA  ZERO=        0.00 


60.00 
DY/OU= 

DN/DU= 


■0* 1525E-02 

0.6710E-03 


AT  SIDE  TO  SIDE  DIST= 
ALFA  ZERO=        o.OC 


90*00 
DY/!)U^ 
DN/DU= 


•0.1098E-02 
0.4636E-03 


AT    SIDE    TO    S IDE    DIST= 
ALFA    ZERO  3C0.0C 


70.00 
DY/DU= 
DN/DU= 


0.34  34E-0  3 
0. 8686E-04 


AT  SIDE  TO  SIDE  DIST= 
ALFA  ZERO-     ~400«00 


80=00 
DY/DU- 

DN/DU"- 


■0.4290E-04 
•0.1671E-03 
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Appendix  IV 
PROGRAM  SOLVE 

I.  Discussion 

II.  Input/Output  Format 
I  ID:.  Flow  Chart 

IV.  Program  Listing 

V.  Sample  Results 
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Appendix  IV 
PROGRAM  SOLVE 

I.   Discussion 

Program  SOLVE,  using  the  hydrodynamic  derivative  inputs 
for  Ship  A,  and  for  Ship  A  in  Ship  B's  position,  calculates 
the  full  set  of  derivatives  for  Ship  B.   The  derivatives  for 
both  ships  are  dimensionalized  in  the  program,  and  the 
applicable  coefficients  for  the  variables  of  the  equations 
of  motion  are  computed.   SOLVE  then  calls  subroutine  DETER 
which  computes  the  values  of  the  coefficients  of  the  polynomial 
in  JJ     which  is  the  result  of  evaluating  the  determinant  of 
the  equations  of  motion  coefficients. 

Subroutine  DETER,  in  the  process  of  performing  necessary 
polynomial  addition,  subtraction  and  multiplication,  calls  on 
library  routines  of  the  IBM  Scientific  Subroutine  Package  (SSP) . 
These  subroutines  are  PADD,  PSUB  and  PMPY  and  are  fully 
described  in  the  IBM  Application  Program  publication  H20-0205-3, 
entitled  "System/360  Scientific  Subroutine  Package,  (360A-CM-03X) 
Version  III".   These  routines  are  included  in  the  program 
listings  in  this  appendix  and  were  used  without  modification. 

Once  DETER  has  evaluated  the  determinant  to  get  the 
characteristic  polynomial  in  S  ,  execution  is  returned  to  the 
main  program,  SOLVE,  to  determine  the  roots  of  this  equation, 
which  are  the  stability  roots  for  the  system.   The  roots  are 
determined  by  a  fourth  SSP  subroutine,  POLRT ,  used  without 
modification  and  listed  in  this  appendix. 
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II.   Input /Output  Format 

Program  SOLVE 's  input  is  ordered  as  follows: 


ITEM 


PROGRAM     CARD     COLUMNS 
SYMBOL 


FORMAT 


Ship  A  Characteristics   ALBP,  ABEAM,  1 

ADRAFT,  ADISPL, 
ACB,  ACP,  ACM 

BLBP  2 

VO  3 


Ship  B  Length 

Equil.  Spd .  in  kts. 

Dimensionless  open 
water  derivatives, 
Ship  A 

(AXU=Xu  ,  etc.) 

Dimensionless 

derivatives  from 
DETERM 


AXU,  AXUD,  AYU, 
AYVD,  AYPD,  AYPOD, 
AYDEL,  ANV,  ANVD , 
ANPD,  ANPOD,  ANDEL 

4-15 

AXN,  AYN,   16-18 
ANN 


1-7  0 

1-10 
1-10 
1-14 


#  Positions  examined    NPOSIT 


Long ' 1  and  Lat . 

separation  at  each 
position 


Dimensionless  Y  , 

a 

Y0 ,  N  ,  NQ,  Ship  A 
6    a    6      ^ 


Dimensionless  Y  , 

u 

N  ,  Ship  A 
u      L 


Dimensionless  Y  , 

u 

N  ,  Ship  A,  Ship  B 


ALFA ( I ) , 
STSD(I) 


19 
20 


1-14 

1-4 
1-20 


Position 


(one  card  for  each  position) 

AYA(I) ,  AYB(I) ,       1-56 

ANA  (I)  ,  ANB(I) 

21 

(one  card  for  each  position) 

AYU (I) ,       22      1-28 

ANU  ( I ) 

(one  card  for  each  position) 

BYUINT(I) ,    23      1-28 
BNUINT(I) 

(one  card  for  each  position) 


7F10.2 

F10. 2 
F10.2 

E14.4 


E14.4 

14 
2F10.2 


4E14.4 


2E14.4 


2E14.4 
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ITEM  PROGRAM      CARD    COLUMNS       FORMAT 

SYMBOL 


Dimensionless  Y  ,  Y  ,    BYAINT(I),    24     1-28        2E14.4 

cu  •   r   ou-   T,a         BYBINT(I) 
Ship  A,  Ship  B 

position 

(one  card  for  each  position) 

Dimensionless  N  ,  N  ,    BNAINT(I),    25      1-28         2E14.4 

e,  •   ,   e.  .    "    P     BNBINT(I) 
Ship  A,  Ship  B 

position 

(one  card  for  each  position) 


The  output  format  for  SOLVE  is  illustrated  by  the  sample 
results  section  of  this  Appendix.   In  the  output,  Y   is 
represented  by  DY/DA,  etc, 
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II.   Solve     plow   chart 

START 


A_ 

/albp, 

ABEAM, 

A  DRAFT, 

AD  15 PL,  ACB, 

AC  P,  ACM 

BLBP 

BBELAM^ 

ABEAM  * 

(BLBP/AL8P) 

BPRAFT- 

-ADRAFT* 

(BLBP  /AL6P) 

BDISPU- 

AC 6  X-  BL8P* 

BBEAMVBDRAFt/35. 

AL3P,  BLBP,   ABEAM, 
SBEAM,   ADRAFT, 
BDRAFi;  AD  IS  PL, 
BDI5PL,  ACB,ACB#ACt3 
f\CP,     ACM,    ACM. 

vo 


■  1   , 

/AXU,  AXUD,  AW, 

AYVD, 

AY  PD; 

AYPDD, 

AY  DEL 

ANV, 

AK'VD 

ANF'D,  ANPDD, 

AN'DEL 

AXK>,  AYM,  fiNK) 


NPOSIT 


ALFA  (I),   ST5DU), 
I«l     WPOSIT 


~i>- 


P 


ROGRAH 


JOL 


50LVE        pG  J 


VF. 


1 


/   AVA(I),   AYB(I), 
AWAIT),  AW6(X), 
AYU(I),  AWUlT), 
BYUIWTCO,   BWUIWT  (I) 
BYAIWT(I),  BVBIWT  IT) 
BWAIWT  (I),   BNBINT  (D 
T=  1.  NPOS  IT 


Dll-tl-  O.5"*RH0* 
ALBP**3 
DlHl=Dlttl*AL8P 

DIM3  ~DIM2  *ALBP 


DIM4-    0.5*RHO*UO* 

AL8P>V<*Z 
DIMS  -    DIM4-*  ALBP 
DIMS  =   DIMS  *  ALBP 


DIM7  =  0.5'*Htto 

*ALBP*UO*tf£ 
DIMS^  DIM?  KALBP 
DIM  9-   DIMS  *ALBP 


/XUD-  AXUD*DlMl 
AYVD-  AWD  *DIM1 
AYPDD^  AYPDD  *  DIM Z 
ANVD^  AIJVD  *  DIMZ 
AlJPDD-  ANPDD  *DIM3 
AXU  =  AXU*  DIK4- 
AYV  ^AYV*  DIM4 
AYPD- AYPD  *DIM5 
ANV^  ANV  KDIM5 
AXU^  AXU  *  D1M5 
AVNf  ^  AYW*    DIMS 
AMPP^  AMPD*  DINKS 
AMU  -  AK)W  *   DIM6 


[> 


V 


AYDEL--  AYOELK-DIMC 
ANDEL^  AW\>EL*blM9 


Cl«  GLBP/AL&P 
Ct^  CI **2 
C3-C1^3 
C4-  Cl  x*Z 
C5-CL**5 


BXUr  AXU*C£ 
BXUO-  AXUDXC3 
BXW-    AXU*  C3 
BYM  =    AYNJ  X  03 
BYDEL-  AYDEL*Ct 
BYV  =    AYV*C2 
BYVD=  A/VD  XC3 
GYPD-  AYPD  v  C'i 
BYPDD-  AYPPD  XC4 
8NW-  ANKI-C1 
BWDEL- AWDEL  XC3 
BNfV^    ANV  *C3 
BWVD^  AWVD  X  C4 
BKJPD^    AKIPD  XC4 
BMPDD^  AUPDD  *C5 


#1.1 


AXO,  6XU,   AXUD,  6XUD,  AW 
BYV,  A.YVDy  8YVD,  AYPD, 
BVPD,  A^PDO,  BYPDD, 

AYD£L,8YD£L,AW\Afl3NV| 
ANVD,   BNVD,  ANPD,8MPI} 
AUPDO,  BUPDD,   AWDEL, 
BNJDF.L,  AXU,  BXU,  AYW, 
BYM,    AMU^BNN 


V 


TO    5oLVF 
PG.  a,  *i-2.i 
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Solve:    pg  'i 


FROM     SOLVE 
PG.       J-, 


t* 


2.1 


NCOWS 


lOOO 

k  «-i,  n  posit 


TRoia      SOLVE 
P&   3,  #  3. 3 


^■Z.4- 


MMK)'-  AYA(K)*  DIM 7 
AveOO"  AVtMK)*  DIM  7 
AYu(|<)-  avu(K)  *  DIM.4 

AWU(k)  =  ANUlK)*0IM5 
AWA  fK^=  ANAfK)<DlM& 
AWG  (K)*  ANB  (K)  X0IM& 


BYAOO"  BYAIWV(K)?- 

C1*DIM7 
BYB(K)=  BYBINT(K)  * 

Cl  ¥   D1M7 
BNA(K)=BNAINT(K)* 

C£  *D1M  8 
BNOUO-BMBINTtK)* 

C£  ¥  DIM  8 
6NJU(k^BYL/IK)T  tK)* 

CZ>'  DIKA 

BNUlK^=  BtvUlvjT  CK>  »• 

C3  *  DIM5 


8ETA(»0  =   STSDlK)  -l 
(ABEAM-*- B BEAM)  It 


ALFA(«)#STSD{WJ 
BETA(K) 


AVAtK),  BVAIK), 
AYB  (K\  BYB(K>, 
ANA  IK),  BNA(K), 
ANB(K>     BK'B(K), 

avu  iK)  ( eyuiKi^ 

AKJU(K);    RNU(K) 


looi 

MJ<—  1  NCOWS 


FROM    Sol.\)E    J 

— >U?.-z 

lfAK(l)  P.K(XJ 


At Tl,  ADT2, 
BDT1,  BDT2 


(M<lM),  M=l;&), 
IBK(M),  M-  i,€), 
A  DTI,  ADTa,BDTi 
BDT2 


V 


V 


AXl(l)  =  AYU  *  BKfe) 
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TT..    SAMPLE      RESULTS 


SI 

■IP    A 

Li}F  = 

5  2  8,50 

BEAM* 

76,  CC 

[  RAI T= 

29.7  5 

r  i  s  p  l  = 

2O9C0.O0 

C  Li- 

0 . 61 2  5 

CF  = 

0.62  46 

0  = 

C .  9  8  C  7 

SHIf 

k  1 5  .  0  0 

68.31 

26.7  A 

1518  2.00 

C  .  6 1 2  5 

C  .  6  2  4  6 

0.9807 


SHIP 


0  X  /  u  u 

-C  .8400 E 

•6 

DX/0UDOT 

-C.  123  4E 

C7 

DY/CV 

-  C  .  8  1 2  3  E 

0  5 

DY/DVDGT 

-C  ,2  2  71c 

C7 

CY/DPDCT 

-C.1846E 

;   ce 

DY/CPCCCT 

0  . 6  6  7  7  E 

07 

DY/DDEL 

C  .  A  8  9  8  E 

06 

DN/DV 

-C.9  748E 

37 

DN/DVCOT 

0 . 1 7  63  E 

0  8 

DN/OPOOT 

-C.3266L 

1 0 

DN/CPCCOT 

-  0 . 3  40  2  E 

11 

DN/COEL 

-  C  .  1  2  r'<  3  E' 

. 

CX/DN 

G.171CE 

CG 

CY/CN 

-C .1937E 

D5 

DN/CN 

Ci5113E 

C7 

SHIP 

■0  .67 
•0.89 
■  C  .  6  5 

-0.1 6 

-0.  1 

0  .42 

0.39 

•G.70 

0.11 

-0.21 

-G.  19 

-0.93 

0.  12 

-G.  14 

0 .  k  5 


B 
8  5  E 

f  : : 
6  1L 

34CE 

57E 
57  E 
77E 

y:  e 

18  E 

c;  5  E 
87  E 
42  E 
C  6  E 
95  E 


CA  . 
06 
C5 
07 

06 
C7 
0  6 
C7 

1" 
11 
C8 
C6 
C5 
n.7 


AT  ALI 


0.0    STSD= 


60.0     Af\C 


THE    POSITICN    DEPENDENT    DERIVATIVES    ARE 


13  2.2 


C  V  /  C  A 

CY/C8 
DN/CA 
CN/C3 

r.  y  /  c  u 

DN/DU 


SHIP     A 
0. 5399E 

-0  .2918E 
0.223CE 
0.  74  )3I 
0.1  06  7E 

-0.  2432E 


CC 
01 

03 

0  3 
05 

0  7 


SHIP  8 
0.  4  85  3E 
C.2623E 
0. 180  IE 
-C..5  9  8CE 
-  0  .  8  6  2  2  E 
C.1PC2E 


CO 
CI 
03 
C3 
CA 
07 
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Sri  I  P    A    CGN  rHCL  SYSTEI  CI         '     NT  S    1    1 

C .  C  0.0  0.0 

C  .  0  0  •  0  C.  0 

SHIP  B  CENT  PEL  SYSTEM  CONSTANTS  1  TC  6 

0.0  O.Q  0.0 

0.0  0.0  0.0 


SI-  IP  A  T  [ME  LAGS  ARE 
SHIP  B  TIME  Lf-GS    ARE 


0.0 

0  .  D 


o.c 

0.0 


TEE    DETERMINANT     IN    ASCENDING    D 

POWERS    IS 

U  •  U 


F  C  L R  T     TEST    V 


CHAR.     EON. 


0.0 

C.  12  55E 

32 

C.510  M 

34 

C. 922 IE 

3  6 

C. 1460E 

3  9 

0.218 6 E 

41 

0. 1949E 

4  3 

C  .  7  6  5  8  E 

4  4 

0  .8479E 

4  5 

C .  2 80 2E 

4  6 

ALUt      IER' 

EAL    RGCT: 

S    A 

C  .0 

-0.1250    - 

-CI 

-  0  ,  -  ■ '  ■  m  6  1  • 

-0  2 

-C.6  807E- 

-02 

-0. L280E 

cc 

-0.7573t  ■ 

-02 

0.173 

-02 

0.  1  73   II 

-32 

-C.  14  24E 


17  8 


SHI  P    A    C  C  NT R C  L     S  V  S  T  E  <v    C  ] N  S  T  /■  NTS     1    T  G 
l.OOi    100  4. CO       00  -C.C1 

-0.100000  0.010000  0.2 


SHIP    B    COM 
1 .000000 

C .  1  c  c  c  c  c 


;L     SYSTt  M    CONST;  N1  S    1     1  C    6 
4  .0000  .     LC   100 

-0.01  0000  -0.2000  )0 


SHIP    A    TIKE     LAGS     ARE 
SHIP    B    TIME     LAGS    ARE 


0  . 0 

0.  0 


o  .0 


THE    DETERMINANT    1  fa    ASCENDING    C 


FOW  ER 


IS 


FCLRT     r  js 


CHAR  .    fcCN 


0.1520E 

34 

0.3168E 

3  6 

0.18  8  3  F 

38 

0.6  835E 

3  0 

0.162 3 r 

0  . 2  6  4  5  E 

42 

0.2997E 

43 

0.2399E 

44 

0.  1652E 

45 

C. 9  746E 

4  5 

C.  2802E 

46 

V^LUE     IER= 

REAL     ROOTS    ARE 

-C.6400E- 

-01 

TC7190E- 

•02 

-0.1149E 

CO 

-0.1745L- 

■01 

-C.  1745E- 

•01 

-0.31 3 1 E - 

-01 

-C.  313]  ,-;- 

-01 

-  0  .  1 2  9 1 E 

cc 

0.3243E- 

-01 

0.3243E- 

-CI 
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Sh  IF    A    CC.M  FCL     SYSTEM    CONST  A  1    I 

l.C(  4.0(  0(  0  0  -o.  OJ 

-0.1C0CCC  0.01  C.2CC00C 

SHIP    B    CONTROL     SYSTEM    CONSTANTS    1     TC    6 
1  .OOOOOC  4,000000  0.C1    C 00 

C.1CC000  -0,01  H  -0.2C 


SUP    A    T  I M  E 
SHIP    B    TIME 


The    DETER! 
POWERS     IS 


FCLR1     TiST 


CHAP.     EGN 


LAGS     APT. 

2  .  C  0 

LAGS    ARE 

2.00 

IN ANT     IN    ASCENDING 

o.is  • 

34 

G.3046E 

36 

C.l£ 

38 

0.549 

39 

0  •  1 2  0  9  E 

41 

C.  102  IE 

4  2 

0.  19  361 

4  3 

0  .  1 4  V  7  E 

4  4 

C.  1084E 

45 

0.71  4 ")  E 

45 

C.2221E 

4  6 

VALUE     IER  = 

0 

REAL    ROOTS    ARi 

-0.6552E- 

-01 

-0.7190E- 

-0  2 

rC.3311E- 

-01 

-0.3  311E- 

-01 

-C.  1124E 

:  d 

-  0  .  1 6  3 1 E  - 

-ci 

-0.1631E- 

-0  1 

-C. 1261 E 

c  0 

C.443  2E- 

■01 

C. 4.432  E- 

-31 

4.CC 
4  .00 


AT     ALFA-  0.0    STSD- 

THE    POSITION    DEPENDENT 


90  .0     A; .-.:    31  TA= 
JERI VAT  IVES    A 


162.2 


D\/DA 
CY/CB 
D  N  /  D  A 
CN/C  > 
CY/DU 
DN/UU 


S  H  I  P     A 

C  •  2  7  3  3  E  0 1'.1 

-0.191 8E  CI 

0 . 1 06 6  E  0  3 

C.475  9E  C3 

0.768  6E  04 

-  c .  1 7 1  o  l  o  7 


SHIP    B 

0.2456E 

30 

0 .  172  41 

CI 

0  .  8  6 1  1 E 

0  2 

-0. 3  844E 

-  C  .  6  2  C  8  E 

C4 

0.124  5E 

07 
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SHIP    A    CC.MBCL     SYSTEM    CONSTANTS     1    TO    6 

CO  0.0  0.0 

0.0  0.0  0.0 


SHIP    3    COI    fRCL    SYSTEM    CONSTANTS    1     TO    6 


0.0 
r    ^ 


0.0 
0.0 


0  . 0 


SHIP  A  T  IME  LAGS  A  ' 
ShIP  B  TIME  LAGS  ARE 


0 .  0 
0.0 


c.c 

C  .0 


THE  DETERMINANT  IN  ASCENDING  D 
POWERS  IS 


0.0 

c.c 

C. 8124E 

31 

0.3  302E 

34 

C.6707E 

3  b 

C  .  1  3  2  9 E 

3  9 

C  .2171E 

h1 

0.  1949 E 

A3 

C.765>!- 

44 

0.8479E 

45 

C.2  302E 

FCLRT  TcST  VALUE  IER- 

CHAR.  FCN.  REAL  ROOTS  ARE 

CO 

C  .0 
--C.1165E-01 
-  0  .  v3  8  3  2  f  -  C  2 
-C.68Q8E-C 2 
-0.1280E  CC 
-C. 7571 E- 02 

C.1296E-02 

C.  1296E-02 
-0.142'tE  OC 
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SHI 
SHIP    B 


1  I  ME 
I  [ME 


LAGS     '■ 
LAGS    ARE 


0  . 0 
0.0 


ri  ,0 
0,0 


THE    DETERMINANT    IN 
POWERS     IS 


ASCENDING  I) 


0,1519 E 

34 

0.3167E 

36 

C  .  1 8  8  3  E 

38 

0,6833 E 

3  9 

C  .  1 6  2  2  E 

41 

0.2645E 

42 

0.2996E 

43 

0  .  2399 E 

4  4 

0.  1652E 

45 

C.9746E 

45 

0.2802E 

4  6 

FCLRT  TEST  VALUE  IER= 


CH/R.  FCN. 


RE/JL 

ROOTS 

,     AR  E 

-0. 

6  3  9  9  E  - 

CI 

"~  u  • 

7190E- 

0  2 

-0  . 

1 1 4  9  E 

0  0 

-c. 

L745E- 

01 

-0. 

1745E- 

01 

3131E- 

•01 

-c. 

313  1E- 

•01 

-0  . 

1  29  .1  E 

0  0 

3243E- 

■01 

c. 

324 3E- 

•01 
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SHI  '"    A    CON  h-  OL     SYSTEM    (  rANIS    1    TC-6 

1  .000000  4 .  -C. D100 

-C.1C0000  0.010000  0.200000 

SHIP    B    CCNTPOL    SYSTEM    CONSTANTS    1    TO    6 
1.  COC   ,:  0  4  .COOOOO  0.010000 

C . 10C0 )0  -0.010000  -C. 2C0C( 


SHI  V 
SHIP 


T]  ME 

T  I M  E 


L  A  G  S 
L  A  G  S 


ARE 

ARE 


2.00 
2.  00 


4.0  0 
4.00 


THE    DE T E R M I N A N T     IN    ASCENDING 
PCKERS     IS 


r-dra    TEST 


0.1519E 

34 

C.  3C4  S£ 

36 

0  .  1 6  3  6  E 

3  8 

C  .  5  4  R  8  E 

3  9 

0.12C8E 

41 

0.  1821 E 

42 

C.  1^3  6E 

A3 

0.1497E 

4  4 

C.l 084E 

45 

C  .  7  1  '.  J  E 

45 

0.2221 E 

46 

VALUE     IEK- 

CMF.  EGN.  REAL  ROOTS  A 
-0. 6552E-C1 
-0. 7190E-C2 
-0.33UE-01 
-  C  .  3  3 1  1 E  -  C 1 
-0.  1124c  00 
-C.1631 E-OI 
-0.  16 3 IE- CI 
-0.1261E  ' 
0  .  4  4  3  1 E  -  C 1 
C.4431E-01 


AT    ALFA=  2C0.C    STSD= 

Tl  E    PCS  IT  I CN    Oc:f;E  MC EN  1 


7 CO    AND    RETA  = 
CERIVA1  IVES    AR L 


1  4  ?  .  ? 


CY/CA 
DY/C  3 
CN/CA 
C  N  /  C  8 

D  Y  /  D  U 

r:  n  /  o  l 


SH  !  P    h 

-0  .6  24 7 E  00 
-C.6854E-C1 

C.28  7  6E  CO 

-0.2293 E  03 

0.23  51E  C4 

C.8342E  05 


SHIP    B 

0 .  5  1  i     1 

CC 

-0.1170E 

01 

-0. 2^clE 

03 

-0.46; 

C2 

0.  1942E 

04 

3.23  1 

C6 
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St- IP    A  CONTRCL    SYSTEM  CONSTANTS     1     1 
CO                            0.0  0.0 

CO  0.0  0.0 

SHIP    B  CCNTRCL    SY ST EN  CONSTANTS    1     1 

0.0  0.0  0.0 

0.0  CO  0.0 


SHIP  A  !"IM£  I  AGS  ART. 
SHIP  E  TIME  LAGS  ARE 


0.0 
0  . 0 


CO 
0  . 0 


THE    DETERMINANT     IN    ASCENCING    C 
PCWERS    IS 


FCLRT  TEST 


CO 

C  .0 

-C. L766E 

31 

-C7C38E 

33 

0.115  7f 

36 

C. 1C44E 

39 

0 . 2 14  IE 

41 

C .  1  9  4  9  b 

43 

C. 7  65  9c 

44 

C.8479E 

4  5 

C. 2E02E 

4  5 

VALUE     lER- 

CHAR.  EGN.  REAL  ROOTS  ARE 


0 


-0 


.  7 577 E- 02 

.6806 E-02 
.1279c:  CC 
.7340E-02 
. 2f34E 
.  1  4  2  4  _ 
.3 93 3 E-02 


C2 

CO 
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SHIP     H    CCNTFCL    SYSTEM    CONSTANTS     1     TO    6 
l.CCCCM  i  4.  COCOOO  -  0.  0] 

-C.1C0C00  0  .01 OCOC  C.2(  DOOC 

SHIP  B  CONTROL  SYSTEM  CONSTANTS  1  TC  6 
i  .COOOCO      4.00C 000      0. C10CCC 

o.  lonoci    -o.oioooo    -c? 


SI-  IP  A  T  IME  LAGS  ARE 
SHI  P  B  1 IME  LAGS  ARE 


0.0 
3.0 


0.0 


Tl-E    LETER.MIN/SNT     IN    ASCENDING    C 
POWERS    IS 


FCLRT    TEST 


CHAR,     EON. 


0  .  1 5 1 9  E 

3  4 

0.316  6 E 

36 

C. 1882E 

36 

0.68 3 1 E 

39 

C,  1622E 

41 

C  . .  2 1 4  4  E 

A  2 

0.2996E 

43 

C. 2399t 

44 

0.1632:; 

45 

0  .97';  i 

4  5 

C.26C2E 

46 

VALUE  I. 

PEAL  root: 

5  ARE 

-0  .  6  39  7  E- 

-01 

-C.7190E- 

-02 

-C.  1149E 

cc 

-0.174 5 E ■ 

-01 

-C.  1745E- 

-CI 

-0.  31 3 IE- 

■01 

-0  .3131  E- 

-01 

-C.  129  IE 

00 

0.3241E- 

-01 

0.3241  E- 

-01 
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SHIR  A  C(  N  I   L  SYS1  ''.■■    CONS!  AN  rS  11  0  6 
1  .CCCOOO      4  .0000  (C     -0.0] 


0.1C00CG 


I 


C. 20C00C 


SHIP  B  CCNTRCl  SYSTEN  CONSTANTS  1  TC 

i .coo:-     4.  x  ooc    c. cio 

0  .100000    -o.cn  0  V    -C.2CC  i 


SUP  A 
ShIP  E 


T  IKE     LAGS    A! 
T  [ME    LiCS     A  i 


2.  00 

2  .  C  0 


':     c    r   C 

i\  .  0  c 


ThE    DETERMINANT     IN    ASCENDING    D 
POWERS    IS 


POLRT     TrST    V 


C  H  /  R .     E  C  N . 


0.151    '£ 

34 

0.304  5E 

36 

0.1635 E 

38 

0. 5486E 

3  9 

C  .  1 2  0  8  E 

41 

o.  ie2  0E 

42 

0.  1936E 

43 

0  ,1497c 

H  4 

C« 1C85E 

4  5 

C.7140E 

4  5 

0.2  22  11: 

4  6 

ALUE     IER  = 

EAl     RCCTS    ARE 

-0. 6550E- 

-  r  i 

-  C  .  7 1  9  C  E  - 

-C2 

-C.3311E- 

-01 

-C.3311E- 

-01 

-C.  U-24E 

CO 

-0.1631E- 

-01 

-C. 1631E- 

-01 

-0.1261E 

00 

C.443CE- 

-CI 

0.44  30 

•01 

AT    ALFA=       -4CC.0    STSD=  80.0     AND    BETA= 

THE     POSITICN    DEPENDENT     DERIVATIVES    ARE 


DY/CA 
D  Y /C  E 
CN/C  ft 
DN/C8 
CY/DU 

i:  I ;  /  c  u 


SHIP     'v 

0  . 1 4  5  7  E 

0.  1492E 

-0.12C8E 

-0. 2523E 

-0. 3333E 

0  .53  93c; 


01 
C3 

03 
C4 
06 


SHI  P    B 

-0.273  4E  CO 

-0.3     :  7 E  C C 

-0.495  8E  C2 

0.1403E  03 

-0.242  6E  C3 

-0  .448  81  C t 
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SI1P    A  CONTROL     SYSTEM  CONSTANTS     L     TC    6 
0.0                            0.0 
CO                           0.0  0.0 

SHIP    B  CCNTRCL    SYSTEN  CCNSTANTS     1    TC    6 
0.0                           CO  CO 

0.0  0.0  CO 


ShIP  A  TIME  LAGS  ARE 
ShlP  B  T  IKE  LACS  ARE 


0.0 


CO 

CO 


TLE    DETERMINANT     IN    ASCENDING    D 
F 0 k E  R  S    IS 

CO 
0  .0 


F  C  L  R  T     TL  ST 


CHAR.     EQN. 


-C.36C9E 

31 

-  C  .  1 4  5  8  E 

34 

C  .8493 E 

34 

C. 985Qh 

38 

C ..  2  1 3  2  E 

41 

0.1948E 

4  3 

C 765  9 E 

44 

0.8479E 

45 

C.2802E 

46 

VALUE    IER  = 

REAL    ROOTS 

\ 

CO 

0  . 0 

-j  C  .  9 1 5  7  E  - 

■02 

-C  l^lWa- 

C2 

-0.6129E- 

■02 

-0. 6129E- 

■02 

-C1279E    00 

0  .3502 E-02 

-C.1424E     30 

-0.6 8 06 E-02 
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SHIP    A    CGNTROL    SYSTEM    CONSTANTS    1     FC    6 
1  .GOOOOO  4,0000  3  -OoJ I 

-C.  100  000  0 ,010000  C.2  30 

SHIP    B    CCNTRGL    SYSTEM    CONSTANTS    1     Tl 
L.CCCCCC  4.00000C  0,010000 

C, 103000  -0.010000  -C.2C000C 


SH]  P 

S  H  I  P 


A     1  I  ?>  E 
B    T  I M  E 


LAGS    ARE 

LAGS     ARl 


0 . 0 
0.0 


0  .0 

G.C 


THE    DETERMINANT     IN    ASCENDING    0 
FCWERS     IS 


0,1.519E 

34 

C  .  3 1 6  6  E 

36 

0 , 1 8 82  E 

38 

C, 6831E 

39 

0.1622E 

41 

0.2644 E 

42 

0.2  996E 

43 

0.23  99E 

44 

C.1652E  4  5 
0.9746 E  4  5 
G.2  802E    46 


PCLRT 

TEST 

VALUE     IEK= 

CH£R. 

ECN. 

REAL    RCCTS     AS 
-C.6397E-01 

_  /"I       11  r'M         <~  O 

-0.1149E    00 
-C.l  745E-01 
-0.  1745 E-01 
-C.3131E-01 
-C. 3131 E- CI 
-0.1291E    00 
C. 3241 E-01 
0.3241E-E1 
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SHIP    A    CCNTRCL    SYSTEN    CONSTANTS    1     T 


1  .CCOOCO 
■  C  .  1 C  0  0  0  0 


4  .  C  C  0  C  0  0 
3.010000 


-  c 


0 . 2 


SHIP    E    CCNTRCL    SYSTEM    CONSTANTS     1    TC    6 
1. CCCOOO  4.000000  0.010     30 

C.  1C00CC  -C. C10CCO  -C.2  00C 


Shi  F 
SHI  P 


TI.V 

T  I  M 


L  d  C  S 

LAGS 


ARE 

■ 


2.00 
2.00 


4  . 0  C 
'«  .  0  C 


THE    DETERI 
POWERS     IS 


FQPT     TEST 


CMP.     ECN. 


IN ANT    IN    i 

\SCENOING 

C.  151 96 

34 

0.3044E 

3  6 

C.163  5E 

38 

C.5486E 

3  9 

c.  120  hi; 

41 

0  . 1 8  2  0  E 

4  2 

0.1S36E 

43 

C. 1497 E 

44 

0.1C85E 

4  5 

C. 7140E 

45 

C  .22?1l: 

46 

VALUE     IE  R - 

0 

REAL    ROOTS     \ 

-  C  .  6  5  5  0  E  - 

-CI 

-C.7190E- 

-02 

-  C  .  3  3  1  1  E  - 

-01 

-0.3  31 1E- 

-01 

-C.  11  241: 

00 

-0.  16  3.1  E- 

-01 

-  C .  1 6  3 1 E  ■ 

-01 

-C. 1261E 

cc 

0.4429E- 

-01 

C.442  9E  - 

-CI 
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